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Symmetriebrechung im Hubbardmodell 
ein bosonisierungsansatz 

zusammenfassung 

Die meisten bekannten Hochtemperatursupraleiter gehoren zur Stoffklasse der 
Kuprate, die sich gut durch das zweidimensionale Hubbardmodell beschreiben 
lassen. Um das Zusammenspiel verschiedener Eigenschaften wie Antiferromag- 
netismus und Supraleitung zu verstehen, berechnet man das Pliasendiagramm 
des Hubbardmodells als Funktion der Ladungsdichte und Temperatur. Fiir diese 
Rechnung eignen sich insbesondere exakte Renormierungsgruppengleichungen, 
die wir im Formalismus der mittleren effektiven Wirkung verwenden. Zu diesem 
Zweck leiten wir eine aquivalente Formulierung des Hubbardmodells her, die die 
Form einer Yukawatheorie besitzt und aus der Informationen iiber langreich- 
weitige Ordnung in verschiedenen Kanalen durch die Berechnung bosonischer 
Erwartungswerte gewonnen werden konnen. Es gelingt uns, die wesentlichen 
Eigenschaften des Phasendiagramms von Hochtemperatursupraleitern zu repro- 
duzieren. AuBerdem zeigt unsere Untersuchung, wie das Mermin- Wagner Theo- 
rem mit der Existenz antiferromagnetischer Ordnung bei nichtverschwindender 
Temperatur zu vereinbaren ist und wie sich die Beriicksichtigung verschiedener 
bosonischer Fluktuationen auf das Phasendiagramm auswirkt. 



Symmetry breaking in the Hubbard model 
A bosonization approach 

Abstract 

Almost all known high temperature superconductors arc cupratcs, which can be 
suitably modelled by the two dimensional Hubbard model. To understand the 
interplay of various long range properties as antiferromagnetism and supercon- 
ductivity, one can calculate the phase diagram of the Hubbard model in the 
charge density-temperature plane. This analysis is conveniently carried out by 
means of exact renormalization group equations that we apply in the formalism of 
the effective average action. For this purpose, we derive an equivalent version of 
the Hubbard model that takes the form of a Yukawa theory. From this modified 
model long range order in various channels can be extracted by simple calculation 
of bosonic expectation values. We are able to reproduce the main features of the 
phase diagram of high temperature superconductors. Furthermore, our analysis 
shows how the Mermin- Wagner theorem can be reconciled with the existence of 
antiferromagnetic long range order at non vanishing temperature and how the 
inclusion of different kinds of bosonic fluctuations affect the phase diagram. 
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Chapter 1 
Introduction 



Tomorrow by the end of the day we shall come to a mountain 
of black stone bight the Magnet Mountain, for thither the currents 
carry us willy-nilly. As soon as we are under its lea, the ship's sides 
will open and every nail in plank will fly out and cleave fast to the 
mountain, for that Almighty Allah hath gifted the loadstone with a 
mysterious virtue and a love for iron, by reason whereof all which is 
iron traveleth toward it. And on this mountain is much iron, how 
much none knoweth save the Most High, from the many vessels which 
have been lost there since the days of yore. 

From: The Arabian Nights, The Third Kalandar's Tale, translated 
by Sir Richard Burton (1850) 



The investigation of electromagnetic properties of condensed matter systems 
is one of the oldest branches of physics. Even the ancients knew about the 
mysterious magnetic force exerted by certain materials, and during the middle 
ages, magnetism was one of the most popular subjects of alchemical specula- 
tion. Over the last two centuries, the fast progress in physics greatly enriched 
our knowledge of possible electromagnetic properties of different materials. The 
long known magnetism is now interpreted as only one possible type of long range 
order, called ferromagnetism. Other long range structures, like antiferromag- 
netic or ferrimagnetic ordering were discovered. On the other hand, materials 
can be classified with respect to their conductivity: Conductors, semiconductors 
and insulators are known. More recently, the discovery of superconductors led 
to completely new developments — both theoretically and experimentally with 
interesting applications. 

Although many of these properties are understood in principle, many unsolved 
problems remain. One of these problems is even the qualitative understanding of 
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Doping 

Figure 1.1: The generic phase diagram for a p-doped cuprate superconduc- 
tor. By AF we denote regions of anti ferromagnetic, by SC regions of d-wave- 
superconducting behavior. The origin corresponds to zero temperature and no 
doping. 

the phase diagram of high temperature superconductors (fig. ( |1.1|) ). This work 
is devoted to develop a formahsm to deal with such systems exhibiting many 
different ways of symmetry breaking and — in the framework of this formalism 
— to shed some light on the origin of the different phases. 
As for all physical problems, we face two kinds of problems: 

1. The question of modeling: The typical high temperature superconduc- 
tor possesses a rather complicated chemical structure. The question is how 
much information about this structure we are allowed to neglect without 
loosing anything significant giving rise to the phase diagram we want to 
explain. By reducing the amount of complexity, we achieve two goals: We 
hopefully end up with a model which can be treated by standard calculation 
methods and that furthermore contains the essential information about the 
actual system in a very condensed form. This should free us from com- 
plications obscuring our view on the true nature of the physical properties 
under investigation. Fortunately, such a model exists, the Hubbard model. 

2. The question of calculation technique: Suppose the Hubbard model 
actually contains enough information to describe high temperature super- 
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conductors. However, it is far from clear how to extract this information, 
given the fact that for nearly forty years the Hubbard model has success- 
fully resisted any attempt to being solved. The main obstacles in doing so 
are 

(a) the different nature of antiferromagnetic and superconducting behav- 
ior, both of which nevertheless should be treated on equal footing in 
our formalism. Besides, we see that for example ferromagnetic or s- 
wave-superconducting behavior is not present in the phase diagram 
of cuprates. Our formalism should not only provide an answer to 
the question why the antiferromagnetic and rf-wavc-supcrconducting 
phases are where they are, but also why other phases are not present. 
We attack this problem by a bosonization procedure. The idea is to 
artificially introduce additional "particles" into the description of the 
model corresponding to physical degrees of freedom like antiferromag- 
netic ordering etc. This allows to discuss antiferromagnetic, supercon- 
ducting and other properties in terms of expectation values of bosonic 
fields — a very convenient method to discuss these phenomena in one 
common language. 

(b) the strong coupling between the electrons. This prevents us from using 
perturbation theory to derive our results. Non-pcrturbative methods 
are needed. Particularly suitable for this kind of problem are renor- 
malization group techniques., which we will apply in the setting of the 
effective average action. 



The focus of this work will be on the bosonization procedure. A lot of recent 
work has been devoted to analyze the Hubbard model in the framework of purely 
fermionic models. We hope to convince the reader that partial bosonization of the 
Hubbard model — giving interesting physical degrees of freedom a particle inter- 
pretation — greatly improves our physical intuition concerning this complicated 
system, which simplifies the motivation of approximation schemes. 

One key ingredient in the formalism to describe antiferromagnetism and su- 
perconductivity in the same formal language will be a certain viewpoint we adopt: 
We interprete both phenomena as spontaneously broken symmetries of the under- 
lying model. Whereas this perspective is quite natural for antiferromagnetism, 
it deserves some explanation in the superconducting case. The usual textbook 
approach to superconductivity is by starting with a microscopic model, motivat- 
ing it by calculating properties of the model and comparing them to experiment. 
To emphasize our point, we will show how the usual properties of a supercon- 
ductor follow merely from the breakdown of C/(l)-symmetry of some underlying 
model that we will not further specify. This topic will be covered in the following 
section of this introduction. The last two sections will discuss high temperature 
superconductors in general and our model for them, the Hubbard model. After 
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a rather formal chapter introducing our starting point, the partition function of 
the Hubbard model, we will describe our bosonization procedure. A mean field 
calculation already reveals the main features of the phase diagram. After that, we 
describe our renormalization group procedure and investigate properties of the 
system beyond mean field — particularly implications of the Mermin- Wagner 
theorem and the influence of charge density and antiferromagnetic fluctuations 
on the superconducting behavior of the model. 

1.1 Antiferromagnetism and superconductivity 
as spontaneously broken symmetries 

Electrons are conveniently described by field operators i'{x), where ip{x) is a 
basis for irreducible linear representations of 

1. SU{2) in the sense that ip{x) = {■iIji^{x),i/ji{x))'^ and the elements of SU{2) 
are represented by U{6) = exp{ia6) and of 

2. U{1) in the sense that ip^x) = ipii^x) + iip2{x) and the elements of U{1) are 
represented by U{9{x)) = exp{i9{x)). 

— * 

Here 9 and 9{x) are used to parameterize the elements of the Lie groups and a 
is the usual set of Pauli matrices. The two dimensional space spanned by ip^ (x) 
and ipi{x) will be called the spinor space. Note that we consider global SU{2)-, 
but local ?7(l)-transformations. The reason will become clear below. 

1.1.1 The S'[/(2)-symmetry 

From rotational invariance in spinor space, we expect the Lagrangian of our the- 
ory to be composed of scalars with respect to SU (2)-transformations. This means 
that the Lagrangian itself is invariant under 5'C/(2)-transformations. The SU{2)- 
symmetry is broken, if for example the (space dependent) expectation value of 
the operator f{x)i/j'^{x)aip{x) does not vanish, where /(x) is an arbitrary non 
vanishing scalar function. In this case, one speaks of the spontaneous symmetry 
breakdown from SU{2) to U{1). 

Ferromagnetic and antiferromagnetic behavior can be inferred from the space 
dependence of f{x). Assume the electrons are strongly located at the sites of a 
quadratic or cubic lattice, so that the above operator expectation value (fiiplaipi) 
is taken at discrete lattice sites i. In the case that {fiiplffipi) is independent of i, 
we have ferromagnetic behavior if /j = fj and antiferromagnetic behavior 

if fi — —fj, where i and j label nearest neighbor lattice sites. 

Since ferromagnetic or antiferromagnetic properties already follow directly 
from the operator expectation value breaking the global SU (2)-symmetry, we do 
not have to bother dealing with the complications of local symmetries. 



1.1. Spontaneously broken symmetries 
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1.1.2 The ?7(l)-symmetry 

In contrast to the S'f/(2)-case we here consider local transformations. This is 
necessary since the defining properties of a superconductor do not follow directly 
from the form of the symmetry breaking operator expectation value. Instead, 
many properties (like the Meissner effect) are connected to the f/(l)-gauge field 
Afj_{x). We therefore assume our Lagrangian to be invariant under the [/(l)-gauge 
transformation 

A^{x) A^{x)+d^Aix), ^^^^ 
ip{x) exp(— ■jeA(x))'?/'(x), 

where we replaced 6{x) by the more common — eA(x) with the electron charge — e. 
A(x) is an arbitrary function, but with A(a;) and A.[x) — 27rn/e, G Z regarded 
as identical. We will now show how the typical properties of a superconductor 
can be derived by assuming spontaneous symmetry breaking of this f/(l)-gauge 
symmetry. This discussion follows [|I|. 

The symmetry is broken if operators like f{x)ip'^{x)ia2'ip{x) develop a non van- 
ishing expectation value (the ia2 between the electron fields is needed since due 
to the anticommutation rules for fermionic fields an operator like f {x)ip'^ {x)ip{x) 
would be equal to zero identically). In this case, the symmetry f/(l) is bro- 
ken down to Z2. Z2 corresponds to the transformations il){x) —>■ ±ijj{x) (with 
A(a;) = or A(x) = — vr/e) which leave f{x)ijj'^{x)ia2ip{x) invariant. According 
to the spatial symmetry of the function f{x) we distinguish s-wave-, p-wave-, 
d-wave, etc. f/(l)-symmetry breaking. Instead of considering the transformation 
properties of ipii^) ^^^^ 'ip2{x) in iIj{x) = ipii^) + i'4'2{x) it is more convenient to 
write 

ip{x) = exp{—ie(j){x)) p{x) . (1.2) 

We define (f){x) to be periodic in vr/e {not in 27r/e as one would expect from this 
notation that resembles the decomposition into an absolute value and a phase 
factor). Then (f){x) serves as a basis for U{l)/Z2 and p{x) for Z2 (if we had taken 
(f){x) to be periodic in 27r/e, (pix) would have corresponded to U{1) and p{x) to 
the trivial group). By this definition, we obtain two fields, one of which belongs 
to the broken U{l)/Z2- and one to the unbroken Z2-symmetry. For U{l)/Z2 we 
have the transformation properties 

A^{x) A^{x) + d^A{x), (P{x) ^(P{x) + A{x), p{x) ^ p{x) (1.3) 

(0(x) -|- A(x) is understood to be taken modulo vr/e) and for Z2 

A^{x) A^{x) + df,A{x), (f){x) ^ (f){x), p{x) -p{x). (1.4) 

Since no mass term tlj''{x)ilj{x) = p^{x)p{x) in the Lagrangian involves 0(x), 4>{x) 
is a massless mode of our theory, the well known Goldstone boson. If we are well 
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within the symmetry breaking regime, the dynamics of the system is dominated 
by this Goldstone boson. In a first approximation, we may neglect the quantum 
fluctuation contributions of the massive modes altogether. Then p{x) plays the 
role of some fixed external field and enters the Lagrangian as a parameter. In 
this case, we end up with a Lagrangian for the electromagnetic and Goldstone 
boson content 

L = -\j d'xF^.F^^^ + L,[A^ - d^ct>i (1.5) 

which is valid in the region of symmetry breaking and not too close to the point 
where the broken symmetry becomes unbroken. The exact form of the functional 
Ls is not known to us; however, the dependence on — d^(j) is dictated by gauge 
invariance. Classically, —Ls can be interpreted as a potential for our theory. We 
will assume that this potential possesses a minimum for vanishing external fields 
An (i.e. the system is stable if external electromagnetic fields are absent) and 
vanishing Goldstone fields, which means that the minimum occurs in ^4^ — 9^0 = 
0. This is all we need to derive the main properties of superconductors. 

We see immediately that if the potential possesses a minimum in A^—d^cl) = 0, 
we have A^ = 9^0, so that the magnetic field vanishes: B = rot A = 0. This is 
the famous Meissner effect Deep within a superconductor we have no magnetic 
field. Closer to the point where the broken symmetry becomes unbroken, i.e. 
closer to the spatial border of the region of superconductivity, A^ — d^cf) no longer 
vanishes. To describe the behavior of the superconductor near the border of the 
superconducting region, we may expand the energy to second order in |A — V0| 
around 1^4 — V0| = 0. The linear term vanishes since we assumed the energy to 
possess a minimum at this point. The quadratic term has the form 

A£;pen = \A-V4>?L^/\\ (1.6) 

where is the volume of the superconductor, A is some length depending on the 
material and in the second line |A — V0p is some average value of [^(x) — V0(x)p 
over the region of integration. Ai?pen describes the energy cost for allowing a 
magnetic field to penetrate the superconductor. Since \A — V0| is of order 5L, 
B being the magnetic field, we have 

A£;pen ~ (1.7) 

On the other hand, the magnetic field carries an energy density of order B^, thus 
the energy cost to expel the magnetic field from the superconductor is 



(1.8) 
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The magnetic field will be expelled from the superconductor, if the energy cost 
to expel the weak magnetic field from the superconductor is much smaller^] than 
the energy cost we have to pay if the magnetic field is to penetrate the supercon- 
ductor: AEcx ^ AEpen or in other words \ <^ L. This means that in materials 
with small A the superconducting region from which the magnetic field is expelled 
is large and vice versa. For this reason, A is called the penetration depth of the 
superconductor. 

Similarly we can see that superconductivity is destroyed, if the magnetic field 
B is larger than some critical magnetic field B^. The fact that some material 
becomes a superconductor means that the superconducting state is energetically 
favored in comparison to the normal state, say by the energy L^A, where A is 
the energy density gap between the superconducting and normal state. As we 
have argued above, the energy cost to expel a magnetic field from the supercon- 
ductor is of order B^L^. If the energy cost to expel the magnetic field is larger 
than what we energetically win by favoring the superconducting state, B > -\/A, 
the material will no longer remain to be a superconductor. The critical mag- 
netic field is then given by Be ~ \/A. However, note that this is only true for 
uniform superconductors. Especially high temperature superconductors are able 
to tolerate much larger magnetic fields than one would expect from these simple 
considerations without losing their superconducting properties. This is due to the 
fact that these materials form magnetic flux vortices, tiny tubes of non vanishing 
magnetic fields that traverse the superconductor. By this mechanism the energy 
cost for expelling the magnetic field is reduced, allowing the material to remain 
superconducting for large magnetic fields ("type II superconductor"). 

We will now come to the most significant property of a superconductor, the 
fact that the resistance equals zero. Imagine a wire made of superconducting 
material with L ^ A, where L is the radial dimension of the wire. Bend the wire 
into a closed ring. Then we know that well inside the wire 1^4 — V0| vanishes. We 
therefore can find a closed curve C (following the linear dimension of the wire) 
along which \A — V0| always vanishes. Now start at some point P on C with 
the fields given at this point by and (pp. Going around the ring following C 
until we reach our starting point Q = P of the closed curve, the fields are Aq 
and (pQ. Since P and Q are equal, we should have Aq = Ap. However, since 0(x) 
is periodic in vr/e, we may have 4'{x)q = (f){x)p + nvr/e, n G all of which are 
equivalent. Therefore the magnetic flux surrounded by our wire is 

/ BhdF = (b Adx= (h V(f>dx = rni/e, (1.9) 
Jf Jc Jc 

where F is the area surrounded by the wire, n is a unit vector perpendicular to 
this area and C is the closed curve connecting P and Q = P. This result tells 
us that the magnetic flux is quantized. A given magnetic flux with tt, 7^ is 

^ "Much" to be on the safe side with all our approximations. 
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maintained by currents flowing in the superconductor. Since there is no way to 
smoothly change the magnetic flux, these currents cannot smoothly decay, which 
means that the resistance of the superconductor is zero. 

The last effect we would like to discuss occurs if two pieces of superconducting 
material 1 and 2 are brought together. Let F be the area of the junction. Then 
the Lagrangian describing the system near the junction is 

Lj = J dxG[A{x),(j)i{x),(P2{x)]. (1.10) 

The integral over x goes over some short line perpendicular to the surface of the 
junction, connecting two points 1 and 2 situated inside the two different materials. 
01 and 02 are the Goldstone modes in the two materials. If we assume that no 
gradients of Goldstone fields and no components of magnetic fields parallel to the 
surface of the junction are present, we may simply write 

L, = FG[l,0i,02], (1.11) 

where we have absorbed the integration over x into G. Gauge invariance tells 
us that G = G'[Ayi0 = dxn{S/(f) — A)], n being a unit vector perpendicular 
to the surface of the junction. The integral is necessary to guarantee the correct 
behavior in the case of vanishing vector potential. In this case 

AA=o0 = y dx^V0 = 02-0i = A0 (1.12) 

so that we end up with a gauge invariant expression as it should be. We want to 
calculate the current fiowing through the junction. The current density is given 

by 

^ = G'(A^0)F^ = -G"(A^0)^ (1.13) 

5A 5A 

and in the case of vanishing vector potential 

J = -G'(A0)^. (1.14) 

The next step is to express A0 by the voltage between the two materials. For 
this purpose, note that the charge density is given by 

^ ' 5A^{x) S(j){xy ^ ' 

so that —J^{x) is the canonical conjugate to 0(a;). In the Hamiltonian formula- 
tion, this yields 



(1.16) 



1.2. High temperature superconductors 
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The voltage V{x) is nothing else than the change of energy density per change of 
charge density, so that 

^{x) = -V{x). (1.17) 

As a side remark, note that this shows that for some superconductor in a station- 
ary state for which = we have V{x) = which is again the zero resistance 
property of a superconductor. If we now assume that our two superconductors 
are kept at a constant voltage and the voltage difference is given by AV, we get 

= -AVt + const. (1.18) 

Using this result in (|1.14| ), we have 

J= -G\-AVt + const)^. (1.19) 

Since A0 is periodic in vr/e, this shows that the current oscillates with frequency 

u = e\AV\/7i. (1.20) 

This is the Josephson effect. It allows high precision measurements of e/h (if we 
had bothered not taking h to be unity), since frequencies and voltages can be 
measured very accurately. 

We would like to recall that all the results we derived in this section were solely 
based on the assumption of a broken f/(l)-symmetry. No explicit dynamical 
model (as the Ginzburg-Landau- or BCS-Lagrangian) was needed to find the 
main properties of a superconductor. This point of view allows us to directly 
identify regions of broken ?7(l)-symmetry in the Hubbard model (which is the 
dynamical model we will use) with regions of superconducting behavior in much 
the same way as we naturally identify regions of broken S't/(2)-symmetry with 
regions of ferromagnetic or antiferromagnetic behavior. 



1.2 High temperature superconductors 

We showed in the last section that superconducting properties can be derived 
by assuming ?7(l)-symmetry breaking of a gauge theory. In this section we re- 
view the history of superconductivity and especially that of the discovery of high 
temperature superconductors. For a recent and more complete overview, see 0. 

In 1911, Heike Kamerlingh Onnes found the first material exhibiting super- 
conductivity by cooling down a mercury wire to 4K. Nowadays we know that 
superconductivity can be observed for many conductors and semiconductors at 
low temperature. However, until 1986 the wide range of materials found to ex- 
hibit superconducting properties had in common that their critical temperature 
Tc (the maximum temperature for which superconductivity occurs) did not ex- 
ceed 20K. In 1986, the first high temperature superconductor (LaBa)2Cu04 with 



10 



Chapter 1. Introduction 



a critical temperature of 35K was found by Bednorz and Miiller . The following 
years witnessed a series of records of critical temperatures for high temperature 
superconductors. The material with the currently highest critical temperature 
known is HgBaaCaaCusOg with = lUK [|. 

Many of the recently found high temperature superconductors are so called 
cuprates, materials with a very anisotropic structure. In contrast to the super- 
conductors known before and to what is usually discussed in dynamical models 
like the BCS-theory, the superconductivity in cuprates has (ia,2_j^2 -wave-symmetry 
0. The cuprates consist of two dimensional Cu02-layers and La-, Sr-, Ba-atoms 
between these layers. For La-interlayer atoms, one effectively finds one electron 
per lattice site of the Cu02-layers. By replacing La bei Sr or Ba, one removes 
electrons from the Cu02-layers, which is called p-doping. Most of the electronic 
dynamics is constrained to the layers. We will exploit this fact by modeling a 
high temperature superconductor by a two dimensional model, neglecting the 
weak coupling between different layers. 

Experimentally, the phase diagram of a cuprate is qualitatively shown in fig. 
( |1.1| ). From inelastic neutron scattering experiments it is known that although 
the antiferromagnetic long range order disappears for strong doping, antiferro- 
magnetic fluctuations are present even in the superconducting domain. There are 
speculations that these fluctuation have an important impact on the supercon- 
ducting order. Furthermore, a whole variety of quantum fluctuations in different 
channels that do not correspond to any long ranged order is under discussion to 
explain the phase diagram. This discussion is additionally fed by the discovery 
of a so called pseudo energy gap below some temperature T*. It is not known 
if this pseudo energy gap is connected to any kind of long range order (as the 
energy gap A discussed in the last section is connected to superconducting long 
range order) , but it is strongly suspected that the key to an explanation of high 
temperature superconductivity lies in the understanding of this pseudo energy 
gap. 

All in all, at the present our understanding of the phase diagram is very 
limited. We hope to convince the reader that the bewildering variety of degrees of 
freedom discussed for the cuprates to explain their properties calls for a formalism 
which is able to include all these degrees of freedom in a transparent and unified 
way and that with our bosonized version of the Hubbard model we are able to 
provide such a formalism. 

1.3 The Hubbard model 

The Hubbard model has been proposed independently by Hubbard, Kanamori 
and Gutzwiller |^ in 1963 as a model for strongly interacting electrons on a 
lattice. There is a wide range of electromagnetic properties of condensed matter 
systems that were or are under investigation by modeling them by the Hubbard 
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model: Ferromagnetism, antiferromagnetism, conductor-insulator transitions and 
— more recently — high temperature superconductivity. The large spectrum of 
physical properties that are tried to be understood by means of this model is 
accompanied by an equally large spectrum of different calculational techniques 
used to approximately solve it. An exact solution of the model is only known 
in one dimension [0]. For two or more dimensions, in general approximations or 
numerical methods have to be used. Unfortunately, the results are not stable 
against choice of the method: A lot of contradicting results have been published 
during the last decades. This is the reason why exact solutions for particular 
values of the parameters of the model play an important role as tests that any 
reliable approximation has to pass. 

The defining features of the Hubbard model are: 

• The electrons are strongly located at the atoms of the lattice. This means 
that the electron field operator is given by ipi, where i label the lattice sites, 
instead of some continuous operator ip[x). 

• The Coulomb interaction between electrons at different lattice sites is ne- 
glected. Any electron interacts only with a possible second electron at 
the same lattice site. Due to the Pauli principle, only two electrons with 
opposite spin at one lattice site are allowed. 

• The electrons have the ability to hop between lattice sites. 

For our purposes, we will additionally make the following assumptions: 

• The lattice is two dimensional and quadratic. This is motivated by the 
actual chemical structure of the cuprates that we want to provide a model 
for. We completely neglect the weak interlayer coupling and the slight 
distortion of the lattice structure away from the ideal quadratic structure. 

• Electron hopping occurs only between nearest neighbor lattice sites. This 
should be the dominating effect, since the electron hopping amplitude be- 
comes smaller with the distance of the lattice sites between which hopping 
may occur. Furthermore, we assume that the hopping amplitude is the 
same for all nearest neighbor pairs. 

With these preliminaries in mind, we can write down the Hamiltonian for the 
Hubbard model 



aj„ and the creation- and annihilation operators for an electron at lattice 

site i with spin a. The first term describes the hopping between different lattice 




(1.21) 



12 



Chapter 1. Introduction 



sites. With our assumptions 



—t , if i and j are nearest neighbors 
, else. 



(1.22) 



The sign in front of t is purely conventional. In particular, we will not assume 
that t > (and indeed there is no simple argument to decide which is the correct 
sign of if:). As we will see, we do not have to bother with this question since all 
our results only depend on t^. The second term of the Hamiltonian describes the 
local Coulomb interaction between electrons at the same lattice site. We take 
f/ > to have a repulsive interaction [U > raises the energy of placing two 
electrons on the same lattice site, which corresponds to a repulsive force). 

The parameters of the model are obviously U and t. If we decide to measure 
all quantities with respect to U, the Hamiltonian may be written as 



Introducing new variables that are dimensionless and measured with respect to 
U, we finally have 



This transcription is unusual in the context of analyzing the Hubbard model in 
this form by means of the renormalization group, as in this case one is mostly 
interested in investigating the flow of the four fermion coupling constants. How- 
ever, in our new approach that we present in this work we will not consider the 
flow of the four fermion coupling, so that this transcription is convenient. 

Another parameter that should be fixed for the model is the number of elec- 
trons on the lattice Ng. If the number of lattice sites is Ng, we have < iVg < 2Ns. 
The case Ne = Ng is called half filling and it is especially interesting, since exact 
results are available for it (at least in the limit of large U). Furthermore, half 
filling corresponds to the undoped cuprate (where because of the chemical struc- 
ture each atom provides one free electron to the system) whereas doping changes 
the number of electrons away from half filling. In contrast to t, the fixed electron 
number does not enter directly into the model as a parameter. We will have to 
include it by specifying a source for the electron charge density^, dynamically 
varying this source to keep the expectation value of the charge density constant. 

The last parameter is temperature. It enters our description when writing 
down the partition function for the Hubbard model as a statistical quantity. We 
will come to the details. 

^This source is nothing else than the chemical potential. 




(1.23) 
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As we mentioned before, rigorous results are important to test approxima- 
tions. For a review, see e.g. |^. One of these exact results is the fact that the 
Hubbard model has an antiferromagnetic ground state at temperature T = 0, 
sufficiently large U and for half filling in agreement with what is experimentally 
found for cuprates (cf. fig. ( |1.1| )). For T > another exact result (the Mermin- 
Wagner theorem) forbids the existence of an antiferromagnetic ground state in 
two dimensions. This is somewhat disturbing, because we would like to predict 
antiferromagnetic order exactly in the region of the phase diagram where it is for- 
bidden by the theorem. A possible explanation would be to argue that although 
the coupling between the layers of a cuprate is weak, it cannot be neglected when 
applying the Mermin- Wagner theorem — for three dimensions, antiferromagnetic 
order is allowed by the theorem. One of the subjects of the last chapter, where 
we analyze the properties of the model using our formalism, is to show that 
it is possible to reconcile the Mermin- Wagner theorem with the occurrence of 
antiferromagnetic long range order for T > even in two dimensions. 

Over the last years, efforts have been made to investigate the properties of 
the Hubbard model numerically by renormalization group techniques [0. In all 
of these approaches the flow of various four fermion interactions was calculated, 
conflrming the main symmetry breaking instabilities of antiferromagnetism and 
d-wave superconductivity in the Hubbard model. These instabilities are inferred 
from the divergence of the four fermion couplings. The divergence of couplings 
at the onset of spontaneous symmetry breaking prevents these approaches from 
following the flow into the broken phase and is one reason for constructing the 
alternative formalism presented in this work, which is more suitable for this task. 
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Chapter 2 

The partition function of the 
Hubbard model 



The starting point of this work will be the partition function of the Hubbard 
model. In this chapter the general method to derive the partition function once 
the Hamiltonian is given in second quantized form is presented. Much of the 
material in this chapter can be found in textbooks covering statistical field theory 
(cf. e.g. 0) and will be known to the experienced reader. However, the last topic 
of this chapter, the formulation of the partition function via coherent states, 
deserves some explanation. We restrict ourselves to fermionic systems. 



2.1 Many particle systems 

The two ingredients for a quantum theory are states and operators. We will 
generalize these concepts from the one particle system to the many particle system 
in this section. First we cover the generalization of state kets. 

Consider a system with identical fermionic particles. The Hilbert space for 
one particle be H.. Then the Hilbert space for the A^-particle system is given by 

Hn = n(^'H(E)...(^n . (2.1) 

N times 

If is an orthonormal basis for the one particle Hilbert space Ti of particle 

i, we can define a basis for Hn 

Oat) = ® . . . IttAf) . (2.2) 

Orthonormality and completeness directly carry over from the one particle basis 
of TC to this basis oiTiN'- 

(ai . . . aTvl a'l . . . a^) = a\) . . . {onI oln) = ^a^a'^ ■ ■ ■ ^a^a'^ 

|ai . . . Oat) (ai . . . oatI = ^ (ail . . . ^ Ioat) (aAr| = 1. ^'^■'^^ 
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For systems with identical fermions, any physical state has to be antisymmetric 
under particle exchange. We therefore define the totally antisymmetric basis 

\ai... aN) = sg"(^) ^^'(1) ■ ■ ■ Q;p(JV)) , (2.4) 



/TV! ^ 

where the sum runs over all permutations of the particles. The scalar product 
now reads 

{ai . . . aN\ a'l ■ ■ ■ a'^) = ^sgn(7') {ai\ 0:^(1)) . . . (aivl ^^'(Ar)) 

V 

1 , if the permutation transferring cti . . . into a[ . . . is even, 
-1 , if the permutation transferring ai . . .aN into a'l - . . a'j^ is odd, 
, else. 

(2.5) 

The completeness relation is 

\ai . . .aN) {oci . . .aN\ — N\ (2-6) 

Ot\ ...OlM 

We now come to the second ingredient of a quantum theory, the operators. 
Suppose we are given a basis of a one particle (labeled with i) Hilbert 

space which consists of eigenstates to some operator t/j 

Ui\Ui)^Ui\Ui), (2.7) 

where Ui is the eigenvalue to C/j. A one particle operator U in the many particle 
system with a general basis {Ictj)} for the ith particle is then defined to be 

JV 

U \ai. . . aN) — Ui \ai . . . aN) , (2.8) 

i=l 

where Ui only acts on the |ai)-part of |q;i . . . a^v)- For example for non interacting 
particles, if we take U to be the energy operator and {Ictj)} to be the energy 
eigenbasis for the ith particle, this means that the energy of the many particle 
system is the sum of the single particle energies. The matrix elements of a one 
particle operator are given by 

N 

{a^... aN\ U \f3i . . . (Jn) ^J^U ^^^l '^*^) ^^^^l ^ l'^^) " (^.9) 

Similarly, we define the two particle operator V by 

V |ai . . .ajv) = ^ 2^ Vij\ai...aN) (2.10) 

l<i,j<N,i^j 
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with the matrix elements 

{a^...a^\V\f3,...M = lY.]l^'^^\ («*«^-| ^ • (2-11) 

2.2 Creation and annihilation operators 

Up to now we considered A^-particle systems, where was some fixed number. 
In quantum field theory however, the number of particles may change. Instead of 
an A^-particle Hilbert space Hn the underlying space is the Fock space JF, which 
is the direct sum of all A^-particle Hilbert spaces 



^=0^^. (2.12) 



N=0 



It is very convenient to introduce creation and annihilation operators on this Fock 
space and to express states and operators by means of these. Since basis kets 
belonging to Hilbert spaces with different N are orthogonal, the completeness 
relation simply reads 



oo ^ 

|0)(0| + 5^— J2 l«i---«iv)(«i...«iv| = 1. (2.13) 

A'"=l ' ai...ojv 

The creation operator is defined by 



a+|Ai...A,v) = |AAi...A^), (2.14) 

transforming a state in Tijv to one in Hn+i- Remember that we only treat the 
fermionic case. For bosonic systems, additional factors appear in this definition 
to guarantee normalization. As a consequence, any state may be written as 

|ai...a^)=a+ ...a+^|0), (2.15) 



where |0) is the vacuum state. Using (2^) and (|2.14|) , we find the anticommuta- 
tion relation 

{a+,a+}=0. (2.16) 
The annihilation operator is defined by 

aa = (a+)^ (2.17) 



From ( p^.l6D we immediately find 

{ao,,ap} = 0. (2.18) 
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By using (^), {^J7^, {^A§ and (gj) we show that 

(—1)*-*^^'' |q;i . . . tti-itti+i . . . On) , if the zth particle is in state 



tta \oii . . . a^) 



, if no particle is in state 

(2.19) 



( |2.14| ) and ( |2.19[ ) yield the last anticommutation relation 

{a^,a+}=6^p. (2.20) 

In the same way we can show that the operator ha = a^cia counts the number of 
particles in the state a: 

N 

ha \ai... Oat) = ^ 6aa, |ai . . • OiN) ■ (2.21) 

i=l 

Of course for fermionic systems, this number is either 1 or 0. The operator 



N = ^na (2.22) 

a 

counts the total number of particles in the system. 

Note that basis changes from a one particle basis {]«)} to another one particle 
basis {!«)} are easily implemented on the creation and annihilation operators. 
From 

\a) = {a\ a) \a) (2.23) 

we find 

a 

One and two particle operators are usually expressed by creation and annihi- 
lation operators. For a one particle operator [/, we have with ( p.4|) and ( |2.9| ) 

N 

(tti ... ttivl f/ ... /?7v) = ^ f/i («!... ajvl /3i .. . (3n) , (2.25) 

1=1 

where we have assumed that is an eigenbasis of f/j. By using the identity 

N 

J2u, = J2Uana, (2.26) 

1=1 a 
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where the sum over a goes over all possible one particle states and Ua is the 
number of particles present in the state a, we conclude that 

U = J2Uaata^. (2.27) 

a 

In a general basis (not necessarily an eigenbasis of Ui) we find by using ( p.24| ) 

^ = 5Z^V«K (2-28) 

with A, /i labeling basis kets of the general basis and 

^AM = $^(A|a)f/a(a|/i). (2.29) 

a 

In much the same way, but somewhat more involved, we can repeat these steps 
to derive the desired form of the two particle operator 

= ^ ^ Vxf.^upO't^^t^p^''- (2.30) 

Xfii/p 

Recalling the Hamiltonian of the Hubbard model ( |1.21| ), we do now under- 
stand the specific form of the terms. The hopping term describes the one particle 
hopping from one lattice site to another, and t is a matrix element giving the 
transition amplitude. The second term is a two particle Coulomb interaction term 
(therefore two creation and annihilation operators) and U describes the strength 
of this interaction. 



2.3 Coherent states 

Up to now, most of our results do not depend on the choice of the basis. We 
will exploit this fact by specifying a very special basis that is useful to derive 
the partition function in the next section. This basis is composed of coherent 
state kets. A coherent state is defined to be an eigenstate of the annihilation 
operator: 

aa 1^) = 1^) . (2.31) 

ipa is the eigenvalue to the annihilation operator aa- Note that since the annihi- 
lation operators for fermions anticommute, the same is true for the eigenvalues. 
This means that the ipa are Gmssmann numbers. If is an eigenket to Oq,, then 
{ip\ is an eigenbra to a+. We call the corresponding eigenvalue t/'*: 



(2.32) 
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Of course, the ip^ are also Grassmann numbers, since the creation operators for 
fermions anticommute. Additionally, we demand the properties 

{V'a,V;} = 0, {iP^,ap}^0. (2.33) 

The full set {■0a, V'a} contains the elements of the Grassmann algebra of all eigen- 
values of coherent states. Note that not only the coherent states no longer 
correspond to states of some definite particle number, but also these states do 
not belong to the Fock space introduced in the last section. Instead, a coherent 
state is a superposition of different kets from this Fock space with Grassmann 
valued coefficients. 

We can now proceed to construct coherent states from the vacuum state, cal- 
culating scalar products, completeness relations and operator expectation values 
using coherent states as we did in the last section for ordinary Fock space states. 

Coherent state kets can be constructed from the vacuum ket by 

m = 11(1 - V^.a^) |0) . (2.34) 

a 

To prove that this is consistent, apply an annihilation operator to both sides: 

13 

= n(l- VJK(1-V'aa:)|0) 

= n(l-^^4)V'a(l-V'aa+)|0) (2-35) 

= lip) . 

In the same way, we can show that we can construct coherent state bras from the 
vacuum bra by 

(V'| = (0|n(l + €««). (2.36) 

a 

It is now straightforward to calculate the scalar product of two coherent states 

(V'|V'')-n(l + C^a)- (2.37) 

a 

The proof of the completeness relation 

I n #:#a n(i - ^'i'p) 1^) (^1 = 1 (2-38) 
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is equally simple, but more lengthy. One proceeds by taking the matrix elements 
of both sides with respect to two A^-particle states in the ordinary Fock space, 
expressing all states by annihilation and creation operators applied to the vacuum 
state. The integral is the usual one for Grassmann numbers 

dijl = j dtp* 1 = 0, jdtpij = jd'tp*tp* = l, (2.39) 

and for multiple integrals the innermost integration is performed first. 

For the partition function we will need the trace of an operator in coherent 
state representation. This trace is given by 

TtA= l'l[drj^^ll{l-r^^p){-^\A\^P). (2.40) 

To prove this, start with the trace in some arbitrary orthonormal basis, insert the 
completeness relation for coherent states and use the completeness of the original 
basis. One also needs that 

di)*dip {a\ ip) {ipl P) = j d^dip (-V-l /?) ("1 ^) ■ (2.41) 

The minus sign comes from the exchange of the integration variables hidden in 
the coherent state kets and bras (cf. ( |2.34| ) and ( p.36| )). 



The last equation we will need in the next section is the expectation value of 
a normal ordered operator A(a^, a^) 

(^1 A(a+, 1^') = + rp%)A{ro., ^«)- (2.42) 

/3 



It follows immediately from the normal ordered form of A and ( p. 37]) 



2.4 The partition function 

The grand canonical partition function is 

Z = Tre-^(^~^^\ (2.43) 

(3 is the inverse temperature /3 = 1/T, is the chemical potential, N is the total 
number operator introduced in ( p.22| ) 

iV = 5^a+a„, (2.44) 

a 

and H is some Hamiltonian expressed by creation and annihilation operators. 
Again, we work on a lattice with sites labeled by i. With a being the spin-3- 
component of the electron (a G {j, |}), we use the collective index a with a = ia. 
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We assume H to be given in normal ordered form. In particular, for the Hubbard 
model we have 



(2.45) 



From ( P?1D| ) we haveQ 



(2.46) 



We cannot apply (|2.42|) directly, since the exponential is not normal ordered. 
To cure this problem, we proceed as usual in the derivation of path integral 
expressions by dividing (3 into M small "time slices" e, so that j3 = Me and write 



AI times 



(2.47) 



Between all of these factor we insert the completeness relation (|2.38|) in the form 



Yldipl, ,.dipa',k e 



li'k) {i'k 



(2.4J 



where k = 1, . . . , M — 1 labels the inserted states. By setting ipM = —i^, i^o = i^, 
iplj = -ip* and = i'*, we have 



a'k 



. k=l 



M 



. A;=l 



IK* 

k=l 



'k\ e 



IV'fe-l) 



M 



M 



nn^Cfc^^".'^ 6-^^^-^'^;^^-'= n^^'^i^"^''"'"^^ i^^^-^) • 



\k=l 



(2.49) 



Since H and N are normal ordered, the same is true for e "^^^ if e is small. 
Then we can use ( p. 42] ) to calculate the expectation value 

(^fc| e"^^^"'^^^ \ipk-i) = e^°'^°''='^"^'='^e"''^^^'^".'=''^"^'="^^"^^'^'^^.'='^°'*-^\ (2.50) 
where 

^(Cfc. ^a,fe-i) = H{a+ ipa,k-i)- (2.51) 



^Recall that for Grassmann numbers 1 — J2j''Pj''p7 = exp(— '0*V'7)- The exponential 
notation is more convenient in the derivation of the partition function. 
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Inserting this resuh in ( |2.49| ) yields 

^ = / ( n n d^kd^p- 

\k=l a 



(2.52) 



Taking the continuum hmit, we get a functional integral expression for the par- 
tition function: 



(2.53) 



Z = / PC(r)P^„(r) 

'V<.(/3)=-^«(0),t/.S(/3)=-^/.S{0) 



exp(- / t/r I >>;(r)(^-/i)C(r) + iJ(C(r),C(r)) 



Particularly for the Hubbard model, we have 

exp dr ijl{T) - /i + ^^(r) (2.54) 



V'<.(/3)=-Vc(0),VS(/3)=-VS(0) 
/3 



2 



where 



Note that the "derivative" ^ is a purely formal transcription of the discrete 
version 

gVar)^^.^C(r)-V.Jr-e) 

since the difference between ifja^k and ipa,k-i (which are Grassmann valued) is not 
"small" in any sense. When we have to actually calculate such a "derivative", 
we will return to the discrete version. 

A remarkable feature of the path integral expression for the partition function 
is the antiperiodic boundary condition ipa{P) = — V'o(O)- If we trace back our 
steps, we see that the anti-periodicity is caused by the minus sign in {—ip\ in eq. 



( |2.40|) that followed from the fact that the ^pa are Grassmann valued. Therefore 



24 



Chapter 2. The partition function of the Hubbard model 



the anti-periodicity is typical for fermionic systems. If we repeat all the steps in 
this chapter for a bosonic system, we would find periodic boundary conditions 
for the functional integral. The anti-periodicity has important implications for 
the properties of a fermionic system and we will discuss it in more detail in the 
last section of this chapter. 

Finally, we clean up our notation by defining spinors 



^.(^) = , ^k^) = (^m^)' m^))- (2-57) 

Using 

E^*-(^)^*(-)(^)^^-(^)^^(-)(^) = -i^KrWr)4{r)A{r), (2.58) 
the partition function becomes 

JVc«(/3)=-V«(0),VS(/3)=-^S(0) 

exp [- dr ^J(t) " + ^) V'i(T) (2.59) 

i 

2.5 Matsubara sums 

The anti-periodicity conditions ■iJjq{0) — —ipa{P) and ■0a(O) — ~'^a{P) ^^^^ ^^^^ 
il^air) and ipai'^) expanded as a series 

oo 

= ^ T^a,nexp(i(2n+ l)7rTr), 

"=-~ (2.60) 
CW= E r<„exp(-z(2n + l)7rTT) 

n=— oo 

with T-independent coefficients, where T — 1/P is the temperature and the factor 
T in front of the expansion coefficients is conventional. These sums are called 
Matsubara sums and one usually introduces the Matsubara frequencies 

cu^ ^{2n + l)7rT (2.61) 

so that 

oo 

^a{r) = E 2"^a,nexp(^a;^r), 

— (2.62) 
= E ^C,nexp(-iw^T). 
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The index F for the Matsubara frequencies indicates that these are fermionic 
frequencies. In the bosonic case, we would have had periodic boundary conditions 
for the functional integral and a Matsubara sum of the same form as for fermions, 
but with o;^ replaced by 



CO, 



^ = 2m:T. (2.63) 



One remarkable difference between fermionic and bosonic system is that for T > 
we have ujf^ > Vn, which is not the case for a;^, since cu^ = if n = 0. We will 
use this positivity property when specifying a regularization scheme for fermionic 
propagators. 



By using 



rl3 

/ dr exp{-i{u^ - u^)r) ^ pSnm (2.64) 
Jo 



and 



Jo 



n,m=— oo 

oo 



n=—oo 

we find for the partition function of the Hubbard model 
Z = [ Vra{r)VMr) 

•/ V-a (/3)=-^a (0) , (/3)=-V'S (0) 

^ ^n=— oo ij 



(2.65) 



(2.66) 



ni...n4=— oo 



One important feature of this expression is the behavior in the high and 
low temperature limit. First note that the Matsubara frequency acts as some 
kind of mass term in the fermionic propagator of this theory. The mass term 
is proportional to temperature. This means that in the high temperature limit 
modes with large n are suppressed. The dynamics of the system is then dominated 
by the Matsubara modes with n = and n = —1, which yield uj^ = ±ttT. 
For our model, it follows that the system is completely two dimensional. As 
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opposed to this, in the low temperature limit the modes tend to form a continuum 
around the zero mode. In this case, the Matsubara sum can be approximately 
replaced by an integral. This renders the system effectively three dimensional. 
The dimensional decrease from T = to T > is called dimensional reduction and 
has an important impact on the qualitative properties of the system. For example, 
the Mermin- Wagner theorem shows that in two dimensions no antiferromagnetic 
long range order can exist. However, in three dimensions this is not true. Because 
of dimensional reduction, this means that for T = antiferromagnetic long range 
order is allowed in our two dimensional model, and indeed there are rigorous 
results proving the existence of an antiferromagnetic phase for zero temperature 
in the two dimensional Hubbard model f^. 



Chapter 3 

Partial bosonization 



The partition function of the Hubbard model describes a purely fermionic model. 
Information about spontaneously broken symmetries is encoded in the renormal- 
ization group flow of the quartic couplings. In principle it is possible to extract 
this information by analyzing the size and momentum structure of these quartic 
couplings. However, the momentum dependence not only reflects interesting de- 
grees of freedom, but also arises from complicated short range fluctuations, which 
we would like to ignore in a simple truncation scheme for solving the renormaliza- 
tion group equations. The aim of this chapter is to explicitly extract interesting 
quartic terms and the relevant momentum structure of their couplings by artifi- 
cially introducing bosonic "particles" corresponding to them. We call this new 
model the "colored Hubbard model". The momentum dependence of the cou- 
plings in this new theory no longer contains essential information about the long 
range behavior of the system, and we will neglect this momentum dependence 
in calculations. We will introduce these approximations later on; in this chapter 
everything is exact, and the "colored Hubbard model" as introduced here is an 
equivalent transcription of the original Hubbard model. 



3.1 The Fermi surface 



To gain a better understanding of the terms appearing in the partition function 
of the Hubbard model, consider the theory with no Coulomb interaction and at 
T -C iQ. In this case, the action of the simplified model can be read off from 



^It is perfectly possible to actually set T = 0. In this case the Matsubara sum in (3.1) takes 
the form of an integral. However, we will never need this transcription explicitly, so that we do 
not bother to write it down here. For our purpose, it suffices to think of T to be infinitesimally 
small. 
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n=— 00 xy 



where we neglected iuj^ which is oc T. We replaced the abstract index i for the 
lattice sites by a two dimensional vector x labeling the lattice sites by two integer 
values for the two spatial directions. By Fourier transforming the fields 



i^xn = J J^^^M exp{txq), 



(3.2) 



we find 



n d q 

n=— 00 """^ ^ ' 

For T <C 1, we know that the chemical potential is equal to the Fermi energy 
yU = Ep and all energy states are filled up to this energy. The second term in 
( |3.3| ) describes the one particle energies E{q) (note that the model has a simple 
one band structure). We conclude that an electron state with energy E{q) is 
occupied, if — 2t(cosgi + cosg2) < /i. A contour plot of E{q) is given in fig. p.l|). 
The equipotential line of quadratic shape corresponds to — 2t(cosgi +cosg2) = 0. 
As is evident from the plot, this is nothing else than the half filling case. We 
therefore identify fi = with half filling, fi ^ describes the doping of the 
system away from the undoped half filling state. 

3.2 The colored Hubbard model 

The starting point of the transcription will be 

Z = [ Vra{r)VMr)exp{-SF-S,oup-Sj) (3.4) 



with 



yn, 



ri.=-oo xy , , 

(3.5) 



^ 00 

Scoup = (/^^ni+n2,n3+n4) V^imV^tcnsV^inaV^i 



2 

711. ..714 = — 00 



xn4 • 
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Figure 3.1: Equipotential lines of the function — cos qi — cos g2- The Unes corre- 
spond to states of the same energy. This plot may be read as a plot of the Fermi 
energies for different values of the chemical potential. 
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Figure 3.2: Labeling of the sites in the colored Hubbard model. The labels of 
the original model are given in parentheses. The coarse lattice sites of the colored 
model are indicated by a ♦. 



We have introduced a "hat" " to indicate fields. Symbols without " will denote 
expectation values of these fields. Additionally, we have introduced a source term 
Sj for fermion fields and fermion bilinears to be able to use Z as a generating 
functional. We will specify 5*-,- later when we need it. Also note that the term 
involving the chemical potential (which has the form of a source term) is now 
included in Sj, serving as a source for the charge density. It is quite natural to 
do so, since we know that the charge density is essentially controlled by doping, 
rendering it a quantity controlled by external conditions on the system. The first 
step in the bosonization procedure is to realize that the most interesting degrees 
of freedom of the Hubbard model have to be implemented non locally. This means 
that if we want to decide whether a system exhibits e.g. antiferromagnetism, we 
have to compare electron spins at different lattice sites. In the same way, we are 
not able to decide whether a system exhibits s- or d-wave superconductivity, if we 
do not take into account the relative sign of electron pair expectation values at 
different lattice site pairs. The idea to deal with this complication is to introduce 
a coarse lattice which consists of plaquettes. Each plaquette contains four sites 
of the original lattice (cf. fig. |3]^). The plaquettes — or equivalently the lattice 
sites of the coarse lattice — are labeled by a two dimensional vector n, which 
takes integer values. The lattice sites belonging to a given plaquette are numbered 
clockwise. We will call these four labels colors. Instead of ipxn, where x labels 
the sites of the original lattice, we have now ipnan, where n labels the sites of the 
coarse lattice, a is the color label with a G {1, 2, 3, 4} and in both cases n is the 
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Matsubara mode. Explicitly, we have 

+l,X2)n) , , 

_ . - _ - (3.6) 

4'nAn = "4^ {xi_,X2-l)ni'^n1,n = '4'{xi+l,X2-l)n 

with n = a;/2, x = (xi,X2). 

The advantage of this transcription is that we can now write down fermion 
bilinears describing antiferromagnetic or superconducting behavior that are local 



on the coarse lattice. Before we do so, we repeat the discussion of sec. |3.1| in this 
new language. 



3.2.1 Fourier transforms 



In a first step, we define Fourier transforms of the spinors as we did in sec. ^TT 
The naive way is to simply define 

^nan= / — — ^V^an (g) exp (mqr) , 

(3.7) 



^Ln = j ^^V'L(g) exp(-mqr) 



in much the same way as before. However, in this case phase factors arise in 
the Fourier transformed expressions in the action, since we neglected the spatial 
dependence encoded in the color index in these transforms. The elegant way to 
perform the Fourier transforms is to define 



with 



V'rian = j ^^^V^an(g) exp(i(n + Za)q), 

d?q - 

i^lan = j j^^lniQ) exp(-i(n + Za)q) 



Zi = (-l/4,l/4) 22 = (1/4, 1/4) 
24 = (-1/4, -1/4)23 = (1/4, -1/4). 



(3.^ 



(3.9) 



Then no phase factors arise in the Fourier transformed expressions of the action. 
Note that the phase factors in this definition are no longer periodic in 27r, which 
means that the same is true for ipaniQ), since the integrand must be periodic in 

ClS db whole. 
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E=/ E-^E I'iB^. (3.10) 



This is a good place to clean up our notation with regard to Fourier trans- 
forms. We define 

Q = {uJn,q), X = (r, n), 
QX = uJnT + nqr, 
"P — — ^ d\ 

727r)2' 

X -"^ n Q n=-oo ^ ' 

6{Q-Q')=f36^,nWS{q-q'), 

5(X-X') = 5(r-r')5„,.'. 

Note that S{q — q') is periodic in 27r and that 5(r) is periodic in (5 for bosons and 
antiperiodic for fermions. The definitions hold for a;„ in the bosonic {ujn = uj^) 
as well as in the fermionic (a;„ = u^) case. We can now write ipnan = "ipaiX), 
'^ani.q) = ^aiQ) and similarly for ^* . 

With these abbreviations, the complete Fourier transforms read 

= ^ exp(z(gX + 

^ . (3-11) 
^l(X) = ^liQ) exp(-z(QX + z,q)). 

Q 

Furthermore, we will often use the notation i/jiQ) for the vector with compo- 
nents ipa{Q)- Note that since the ipaiQ) themselves are two dimensional spinors, 
this means that the objects ipiQ) live in the product space of spin and color and 
have 8 components altogether. 

Another notation we will use concerns infinite sums like Ylx Note that 
(mathematicians hopefully forgive this) 

J2 = J2Y1 ^(^) exp(^QX) = J2 ^(QMQ) = ^(Q = 0)- (3.12) 

X Q X Q 

We define 

V = ^ = <5(Q = 0), (3.13) 

X 

which can be interpreted as the two dimensional volume of the system (that we 
assume to be large) divided by temperature. 

3.2.2 The Fermi surface 



We can now repeat the calculation of sec. p.l| , taking U = 0, T <^ 1 and all 
sources except /i equal to zero, so that we end up with 

S = Sp + Scoup ~l" Sj 

= i^\Q) (-Ai - 2t (cos(gi/2)Ai + cos(g2/2)5i)) ^(Q). (3-14) 
Q 
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Ai and Bi are 4 x 4-matrices in color space and are defined in the appendix [A. 2 



This result should be compared to (|3.3| ). The first difference to be observed is that 
the cosines are no longer periodic in 2tt, but in Air, which is a direct consequence 
of our Fourier transform. This means that all the cosines are positive in the 
interval of integration. 

To see how the Fermi surface emerges in this picture, we temporarily switch to 
fields "^{Q), for which the fermionic propagator term in ( p. 141) becomes diagonal. 
For these fields, the action reads 

S = J2'^\Q)D^{Q) (3.15) 

Q 

with 

/-/i + 2t(ci + C2) \ 

-/i + 2t(ci - C2) 

-fL- 2t(ci - C2) 

V -fi-2t{ci + C2)J 

(3.16) 



D 



Ci = cos(gj/2). 

We can immediately read off the (cross shaped) Fermi surface for half filling 
from this result (fig. p.3| ). 

The reason why we are so interested in the shape of the Fermi surface is the 
following. Note that the matrix D is nothing else than the propagator matrix 
of four distinct fermion modes (keep in mind that the chemical potential will be 
absorbed in a source term and is not regarded as part of the propagator). Zeroes 
of the propagator are a well known problem in any quantum field theory calcula- 
tion for massless particles, since they lead to divergencies. These zeroes appear 
in our case on the Fermi surface. The usual way to deal with the divergencies 
is to define a regularization scheme. This is not that difficult in simple theories 
where the propagator (in a Euclidean formulation) is proportional to the squared 
momentum — any positive mass like term added to the propagator will cure the 
divergency problem. However, given the complicated (cross shaped) momentum 
structure of the Fermi surface we face in our formalism, momentum cutoffs be- 
come tedious to define and to work with. We will pursue a different way by noting 
that if we do not neglect iu^ in the action above, we end up with a propagator 
of the form iu^ + D, which does not vanish for all T > 0. We exploit this fact by 
using a regularization scheme that uses temperature as a flowing cutoff function. 
By defining some (unphysical) temperature as a function of a parameter k, 
lim^^oo^A: = 00, limfc^o^fe = T, we can lower starting with some large k and 
letting /c — i> in a controlled way until we reach the physical temperature T of the 
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Figure 3.3: The Fermi surface for half fiUing in the Hubbard (dashed hne) and 
the colored Hubbard model (full line). 
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system. The flow of the system with k will be described by exact renormalization 
group equations. We will come to this later. 

Another aspect concerns the motivation of approximations. For the bosonic 
propagators (introduced in the next chapter), we will be able to expand trigono- 
metric functions to quadratic order. This is very convenient, since one achieves 
formal agreement with known theories (e.g. the propagator of a boson 2(2 — 
cos(gi) — cos(g2)) becomes ^ in quadratic order). However, it is not possible, 
even in principle, to expand the trigonometric functions in the fermionic case for 
low temperature without loosing significant information, since the dynamics is 
dominated by modes with energy close to the Fermi surface, not just by modes 
with zero momentum. The main complications of calculations we are about to at- 
tack are that we are forced to keep these trigonometric functions in the fermionic 
sector. 



3.2.3 Symmetries 

In this section we discuss the various symmetries of the colored Hubbard model. 
As already discussed in the introduction, we have the U {!)- symmetry 

i^ix) ex.p(ie)'ijj(x), ^^(a:) ^ '0^(^) exp(-ie) (3.17) 

and the SU (2) -symmetry 

i:{X) ex.p{iae)ilj{X), iP^X) i/j^X) exp(-ia^). (3.18) 

Note that now we only consider global [/(l)-transformations, since no gauge 
bosons are present in our theory (the reason being that these are "integrated 
out" under the assumption of negligible interaction of electrons at different lattice 
sites, giving rise to the Hubbard model as an effective purely fermionic model). 

Furthermore, the model possesses the symmetries of the underlying lattice. 
These may be composed from translations, rotations and reflections. We restrict 
ourselves to the translation by one lattice site in the positive 1-direction, 
the counterclockwise rotation it! by 90° around the origin (at the center of a 
plaquette!) and the reflection I at the the 2-axis containing the origin. All 
other lattice symmetries can be built up by products of these three symmetries 
(for example, the translation Ty along the 2-direction can be composed by a 
translation and a rotation R). 
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In position space, these symmetries act as 

ni,n2)ln 
^(ni ,n2)2n 
'^(ni,n2)3n 



T 

-J- T! 



R: ( 



^ V'(ni,n2)4n 
ni,n2)ln 
ni,n2)2n 
'^(ni,n2)3n 
^ V'(ni,n2)4n 
V'(ni,n2)ln 
ni,n2)2n 
ni,n2)3n 
>. V'(ni,n2)4n 



^(ni,n2)2n 

rai+l,ra2)4n 
rai,n2)3n) 
— n2,ni)4n 
— n2,ni)ln 
i'{-n2,ni)2n 
— n2,ni)3ra; 
^(-ni,n2)2n 
— rti ,n2)ln 
— ni,n2)4n 
— ni,n2)3ra5 



The same apphes for i/)*. Fourier transforming yields 



R : <^ 



^(g) ^ Ai^(Q)exp(zgi/2) 

^t(g)^^t(g)^^exp(-zgi/2) 

/O 1\ 

10 

10 

\0 1 0/ 



'4'niqi,q2) 



'ipn{-q2,qi] 



V^n(9i,92) i^i{-q2,qi) 



/O 1 0\ 

10 

1 

\i oy 



/ : 



^n(gi,g2) 

^n(?l>?2) 



^iV^n(-gi,g2) 
■^i{-qi,q2)Ai 



Again Ai is one of the matrices defined in the appendix [A.2 
It is interesting to consider a transformation like 



Ml 



(3.19) 



(3.20) 



(3.21) 



which is also a symmetry of the action. This transformation does not possess 
an interpretation in position space by means of lattice symmetries as one would 
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expect. Instead, it reflects our freedom in choosing how to label the sites. For 
example, formally translating this transformation back to position space, we end 
up with e.g. ip(ni,n2)in ~^ '4'{ni-i/2,n2)2n- We may read this as the identity trans- 
formation, written by merely choosing another origin of the coarse lattice, shifted 
one lattice site of the original lattice to the left. This class of symmetries is there- 
fore new in the colored formulation and corresponds to the unity transformation 
in the original theory. We call this kind of symmetry relabeling symmetry. In the 
same way we can define the relabeling symmetry transformations La^, Lbq and 
Lb^i which are derived from La^ by replacing the matrix Ai by Aq, Bq and Bi 
respectively. These four relabeling symmetries correspond to the four different 
ways to assign color labels by choosing a different origin of the coarse lattice. As 
in the case of Lai we see that by multiplying with appropriate momentum phase 
factors these additional relabeling symmetries correspond to a translation in the 
2-direction for L^^, a translation along the diagonal connecting color sites 1 and 
3 for Lbq and no translation at all for Laq- 

The last symmetries we mention are reminiscent of time reversal symmetries. 
They are realized by 



where Mi e {A2, B2, B^} and are denoted by Ta^, Tb^ and Tsg. 

3.3 Partial bosonization 

3.3.1 Definitions of fermion bilinears 

Using the color notation, we are now able to define fermion bilinears correspond- 
ing to interesting degrees of freedom in a simple way. In particular, we want to 
include bilinears describing the charge density, antiferromagnetic order as well 
as s- and (i-wave superconductivity. After defining these bilinears, we try to de- 
compose the four fermion action of the colored Hubbard model with respect to 
these bihnears. We see that this decomposition is not possible until a whole set 
of additional bilinears is added. We define this set in this section and discuss the 
decomposition of the four fermion interaction of the Hubbard model with respect 
to this set in the next section. 
We start by defining 



(3.22) 



XabiX) 



-^j(x)m2^:(x). 



(3.23) 
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The operators aab{X) are uncharged spin singlet operators, (/^^^(X) uncharged 
spin triplet operators and Xab{X), Xabi-^) charged spin singlet operators. Note 
that XabiX) = XbaiX) and xlbi^) = xLi^)- ^^r sake of simphcity, we do not 
take into account charged operators in the spin triplet. 

Suppressing the X-dependence, we now define the composite bilinears 

P = 5-11 + 5-22 + Cr33 + 5-44 

p — an — 5-22 + ^33 — 5"44 
% = 5-11 + 0-22 - 5-33 - 0-44 
Qx — — 5-22 — 5-33 + 5-44 

S = Xll + X22 + X33 + X44 

C = Xll - X22 + X33 - X44 
ty = Xll + X22 - X33 - X44 
ix = Xll - X22 - X33 + X44 

s* = Xn + X22 + X33 + Xli 
c* = Xn - X22 + X33 - Xu 

— Xll + X22 ~ X33 ~ X44 
^x — Xll ~ X22 ~ X33 + X44 

The bilinears that we wanted to include in our formalism are the charge 
density p, the antiferromagnetic spin density a, the superconducting s-wave s 
and (ia,2_y2-wave d. In order to make the before mentioned decomposition of the 
four fermion interaction into these bilinears possible, the rest of the bilinears has 
to be additionally included. This set is minimal and cannot be further reduced. 
However, by concentrating on p, a and d (that is, dropping s from the list of 
bilinears we want to include), we may reduce this set by dropping s, c and ix/y. 

Many of the bilinears we considered only in order to be able to perform the 
decomposition in the next section have a simple physical interpretation. For 

— * 

example m describes the ferromagnetic spin density, and the various charged 
composite bilinears correspond to d^^j^-waves (c), extended s-wave (e) or p- waves 

i^^x/y)- 



^ = <^11 + <^22 + ^33 + ^44 
« = <^11 - <^22 + <^33 - <^44 
9y = <fll + <^22 - <^33 - <^44 
9x = <^ll - <^22 - <^33 + 'fu 

e = X12 + X23 + X34 + X41 
^ = X12 - X23 + X34 - X41 

Vy = X12 - X34 

Vx = X23 - X41 

e* = Xt2 + X23 + X34 + xll 
d* = Xt2 - X23 + X34 - xll 

Vy = X12 - X34 

!}C + 

'^x = X23 - X41- 
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3.3.2 Decomposition of the four fermion interaction 

The four fermion interaction in the colored Hubbard model reads 

X a ^ 
X a 

Scoup may be decomposed into our fermion bilinears by use of the identities 
4^^- = /^'+^' + ^' + ^ 

a 

-12 E ^aa = a +m + gl + fy 



(3.26) 

a 

= ^{-p^ +p'^ + m - t) + d* d + e*e + 2{vlv^ + vlvy). 
To prove this, note that 

X*abXcd = ^ca^db + CFcbC^da 

s n _ v3-2(j 

'^ab'^cd = -^ab^cd - '^aadO'cb- 

The clue of ( p.26| ) is that this list of possibilities to write down combinations of 
fermion bilinears to either give a multiple of Yla ^aa is exhaustive. No other 
independent combinations of fermion bilinears can be found to give Yla ^aa 0- 

3.3.3 The partial bosonization 

Consider the fermionic partition function of the colored Hubbard model 

Z = exp l^- ^^^y^ ^ J V^*V^ exp(-S'p' - S^oup - Sj) 

Sf = Y^^KQ) {^^n - 2t (cos(gi/2)Ai + cos(g2/2)Bi)) V'(Q) 
Q 

'^-p = ^EE4(^)4w4w4w = ^EE^-(^)' (3.28) 

X a X a V • / 



40 



Chapter 3. Partial bosonization 



We already know the kinetic and the coupling term. Additionally, now we have 
specified an explicit form of the source term. In this term, we included sources 

G {Ip = lp + fx, Ip, Iq^^y, Irn, ^a; ^g^/y} 
jf3 £ {is^ ic^ it^/y^ ie^ id^ iv^/y] (3.29) 
j (3 ^ {js 5 jc 5 ^t^/y 5 je 5 id 5 3v^/y } 

for all the bilinears introduced in ( p. 241) with 



^ {p,P,Qx/y,rh,a,g^/y} 



(3.30) 



so that denote the uncharged and up the charged bilinears. We also added 
terms quadratic in the sources. Since physical properties are not affected by these 
quadratic terms (which only give rise to a field independent factor to the partition 
function), we have the freedom to do so. Note that the chemical potential is now 
part of one of the sources, as we discussed at the beginning of section |3]^. In 
the case of vanishing sources except for the chemical potential /x, we demand the 
only term to survive to be the one linear in /i. This demand gives rise to the 
quadratic term in /i we added (and wrote as an exponential factor in front of 
the partition function) to cancel the contribution quadratic in /i from the term 
^/^(X)^. The quantities hp, are arbitrary at the moment. We will come to 

them soon. Additionally, 5*^ contains a source term Sj for the fermions. 
We define the partially bosonized partition function by 

Z = exp (^- ^"^J^ ^ J V^*Vi>Vu*VuVw expi-Sp - Sb - Sy - Sj) 

Sf = Y.'^KQ) i^^n - 2t (cos(gi/2)Ai + cos(g2/2)5i)) i,{Q) 
Q 

X \ p 7 / 

Sy = -Y. [Y.^p{ul{X)up{X) + up{X)u}{X)) + h,w,{X)w,{X) 
X \ p 7 

= - E f E ^mX)^^) + 3piX)uliX)) + Y ^l,{X)w,{X)\ + sf 

X \ p 7 ^ / 

(3.31) 



In this partition function, we have a fermionic kinetic term Sp which coincides 
with the corresponding term in the Hubbard model. In the remaining terms we 
introduced bosonic fields ii, u* and w^, one for each fermionic bilinear -u, u* and 
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7- 



Sb is a mass term for these fields. Sy describes a Yukawa like coupling of 
the bosonic fields to the corresponding fermionic bilinears with Yukawa couplings 
hff, h^. The source term now provides sources for the bosonic fields. 

The next step is to prove that this partially bosonized partition function is 
equivalent to ( ^.28 ) for appropriate values of the Yukawa couplings. We realize 
that since the action is quadratic in the bosons, the bosonic functional integral can 
be performed as a simple Gaussian integral. Recall that Gaussian integrals can 
be evaluated by evaluating the exponent at its stationary value. Since the bosonic 
propagators are mass like, we can neglect the constant pre-factor altogether and 
only have to insert the stationary values 



up{X) 
u%X) 



4v]-2 



h 



h 



u*JX) + 



hf3 



h 



(3.32) 



/3 



^ w^{X) + ^^^^^ 



47r2^' 



h. 



'7 



into the exponential of the partially bosonized partition function in order to 
perform the bosonic functional integrals. We see that the source term reduces to 
the source term of {J?2 



Sb + Sy gives a quartic fermionic coupling term of the 



form 



X \ f3 7 



(3.33) 



If (|3.28| ) and (|3.31|) are equivalent, we must have 



X 



»UXf^ (3.34) 



X 



This equation should be read as a condition on the Yukawa couplings ha, h^. At 



this point the identities (|3.26| ) come in handy. We can use them to parameterize 
the solutions of ( |3.34| ). The general solution is further restricted by the fact 
that we demand the couplings to be real. Then the general solution is (with 



K - h^ - \H^) 



H„ 



H. 



Qx/y 



3(A2 - A3) 
3(A2 + A3) 
3Ao 



H, 



9x/y 



2Ai + A2 + 3A3 + 1 
2Ai + A2 - 3A3 + 1 
2A1 + A2 + I, 



(3.35) 



Ht 



/y 



2 ' 



2He = 2Hd = 



•■Iv 



6 A3. 



(3.36) 
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The parameters Aj obey 

A, > Vi 

As > A3 (3.37) 
2Ai + A2 + 1 > 3A3 

to guarantee that the condition hf3, /i^ G R is fulfilled. 

For any choice of the parameters Aj meeting these conditions, the partially 
bosonized partition function is equivalent to the fermionic partition function of 
the Hubbard model we started from. This means that the choice of the Yukawa 
couplings contains a lot of arbitrariness — which does not matter in the exact 
transcription we used here, since all choices are equivalent to the original Hubbard 
model. However, if we use approximations in calculations, the results can and 
will depend on the initial choices of the couplings. This can serve as a test for 
approximation schemes — an approximation is regarded as well justified, if the 
results do not depend on the initial choice of the couplings. However, this problem 
will remain a disturbing one and is the weakness of our theory that we traded in 
for the possibihty to investigate the properties of a system by direct calculation 
of expectation values of bosonic fields. 

Note that is not possible to bosonize the theory without taking into account 
the spin triplet bilinears (due to the signs in ( |3.34] ) and ( |3.26|) ). 



3.3.4 Symmetries 



Most of the symmetries discussed in sec. |3.2.3| can be easily implemented in the 
partially bosonized version of the colored Hubbard model. As we know how the 
fermionic fields transform under the symmetry transformations, we can derive the 
transformation behavior of the fermion bilinears involved in the Yukawa coupling 
terms in ( p.31j ). Since we know that ( |3.31| ) and ( |3.28| ) are equivalent, they should 



have the same symmetries and we define the behavior of the bosonic fields un- 
der symmetry transformations such that the partially bosonized action becomes 
invariant. 

In particular, we find for the f/(l)-symmetry ( |3.17| ) 



w^, Uj3 — > exp{2i9)ui3, u*p exp(— 2z6')u^. (3.38) 

Similarly, for the SU (2)-symmetry we obviously have invariance of all spin singlet 
bosons. The spin triplet bosons transform as three dimensional vectors under 
S'0(3)-rotations around the ^-axis with rotation angle 2\6\. 

The rotation R, refiection I and translation are also implemented in an 
obvious way. However, simple translations cannot be defined in our present 
formulation, since these translations correspond to shifting the field by half a 
lattice site of the coarse lattice and the bosons are defined on the coarse lattice 
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only. To preserve invariance under translations T^, we define a color label a for 
the bosons and bilinears and set 



w^3{X) = Wy{X), Wy4{X) = T~^Wy{X), 



w^i{X) = TyT^ ^Wy{X),w^2{X) = TyWy{X), 
w^3{X) = Wj{X), Wj4{X) = T~^Wj{X), 



(3.39) 



(3.40) 



and similarly for the bosons up.,u*p and the boson and bilinear sources. The trans- 
lations Tj. are now simply implemented by e.g. T^t&^i(X) = w^2{.X), Tj.w^2(-^) = 
w^i{X+ei) etc. Note that with these definitions, the lattice symmetry operations 
may be written in a way completely analogous to ( |3.19| ), since the bosonic fields 
also live on a coarse lattice with the same lattice spacing as the fermions, with 
four color labeled fields attached to each coarse lattice site. One consequence 
is that for the same reason as in the fermionic case we introduce an additional 
relabeling symmetry that acts in the same way on the bosonic color space as 
it did on the fermionic color space. Again, these relabeling transformations all 
correspond to the identity transformation in the uncolored formulation. 

To implement bosonic color into our partially bosonized partition function, 
we write 

Z = exp (^- ^^J^ ^ J ViP*VipVu*VuVw exp(-5^ - Sb - Sy - Sj) 

Sf = 5^ V^^(Q) {^UJ^ - 2t (cos(gi/2)Ai + cos(g2/2)5i)) ^{Q) 
Q 

a X \ (3 7 / 

= -i E E [J2^^(^UX^^P-(X^ + UpaiX)u;,iX)) + J2 h,W,a{X)w,a{X) 
a X \ (3 7 , 

= -z EE f E ^^(i;a(^)^/3a(x) + j,,(x)«;,(x)) + ^UX)W,^{X) 

^ a X \ (3 '^'^ 7 ^ 

(3.41) 



instead of (|3.31|) and realize that the same calculation that we performed to show 
that ( |3.31|) and ( p.28|) are equivalent also goes through here. 
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3.3.5 The colored partition function 

We finally summarize our results for the partition function. In position space, we 
have 

Z = exp (^- ^^^^^ ^ j ViP*VipVu*VuVw exp(-5^ - Sb - Sy - Sj) 

Sf = J2 ^^(^) - 2^ (cos(gi/2)Ai + cos(g2/2)fii)) ^{Q)) 
Q 

= -i E E ( E /^/3(«;a(^)^/3a(X) + ^^,(X)t2;,(X)) + h,W,a{X)w,a{X) 
a X \ (3 7 

= - EE [Y.^JU^)npa{X) + J^,(X)«;,(X)) + ^L,,(X)u;,,(X) ) + S'^ 

a X \ 13 7 / 



(3.42) 

We have redefined the sources by setting 

J^,(X) = ^j^,(X), J*,(X) = 0;,(X), L,„(X) = ^/,„(X). (3.43) 

The partition function can be used as a generating functional for the bosonic 
ra-point functions. In particular, we have 

U(3a = {Upa) = In Z 



6 



a) = T^^^Z (3.44) 



''-^70 



for the expectation values of the bosonic fields^. If we rewrite ( p.28| ) using colored 
bilinears and insert the resulting partition function in the right hand side of 
( |3.44| ), we find the relation between expectation values of bilinears and bosons 



h} . 1 

U/3a — {Ul3a) H nJ/Ba 



1 



> + i-A (3.45) 
47r2 ■ 7r2 



^Note that L* (X) = Lya{X), but L* (Q) = L^a{—Q) in momentum space. If we write 



( p. 44 ) as we did here, the equations hold true both in position and momentum space 
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These relations show that for vanishing sources the expectation values of the 
bilinears and the corresponding bosons are equal up to a factor. This is the reason 
why our formalism makes sense: We know (by construction) that the bilinears 
describe interesting properties of the fermionic system. The expectation value of 
the bilinears tell us whether the system exhibits e.g. antiferromagnetic behavior. 
This is the case if (fia) is non vanishing, which in turn means that also the bosonic 
expectation value does not vanish. We therefore can analyze the properties of 
the system by merely calculating expectation values of bosonic fields in a Yukawa 
like theory. In principle, all calculations in this work are dedicated to do exactly 
this, and to interprete the results by means of the underlying fermionic theory. 

Again we stress the special role of the charge density pa- The source of the 
charge density contains the chemical potential and will in general not vanish. 
The bosonic expectation value of pa is 



Using the Fourier transforms ( p.ll ) for the fermions and defining up and to 
have the same Fourier transform as if), and u*^ to have the same Fourier transform 
as if)*, we can write down the partition function in momentum space 

Z = exp l^- ^^^y^^ ^ j Vi>*VijVu*VuVw expi-Sp - Sb - Sy - Sj) 

Sf = 5^ V^^(Q) {too^ - 2t (cos(gi/2)Ai + cos(g2/2)5i)) ^{Q) 
Q 

a Q \ (3 7 / 

Sy = -J2J2 kq-q'-q") 

abc QQ'Q" 

E {^UQ)^a{Q')vZ{Q\ Q")UQ") + Q"W,{Q": 



-rY,w,mi^l{Q')vZ{Q' ^ -Q")M-Q": 



a Q \ f3 7 / 

(3.47) 



The vertex factors Vab,c in the coupling term are given in appendix |B[ 
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Chapter 4 

A mean field calculation 



Before continuing to improve our formalism, we want to give a first impression of 
the power of our formalism even in a very simple mean field like approximation. 
The results of this chapter have been published in [jl^] . The main ingredients of 
this approximation are the following: 



• The bosonic fields act as constant background fields, so that the bosonic 
functional integrals can be trivially performed by simply setting all bosonic 
fields to these constant background values. 



• The bosonic background fields are homogeneous in the sense that they do 
not possess any spatial dependence. 



• All background fields except the charge density p, the antiferromagnetic 
spin density a and the superconducting (i-wave d vanish. 



We proceed by deriving an expression for the effective action in this approx- 
imation, which depends on the values of the three non vanishing background 
fields. For given temperature and charge density, we then look for minima of 
the effective potential with respect to the expectation values of d and d. If the 
minimum occurs at non vanishing expectation values of the antiferromagnetic or 
superconducting fields, we conclude that the system exhibits antiferromagnetic 
or superconducting behavior at the given temperature and charge density. 
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4.1 The effective potential 

Using ( p.44| ), we define the effective action as the Legendre transform of the log 
of the partition function 



T[U/3a,u}^,W^a] = -\n Z[J pa, J^a^ L^a] 

a Q \ (3 7 / 

(4.1) 

where the sources on the right hand side are functionals of the expectation values 
Uj3a, u*pa^ ^71' ^he usual properties of a Legendre transform tell us 

that 

- r - T - T* (A')\ 

bUpa OU*pa ^^-ya 

For vanishing sources these equations are formally nothing else than the classical 
action principle — therefore the name effective action. They may be regarded 
as the equations of motion for the field expectation values, taking quantum cor- 
rections into account. Assuming that the fermions have been integrated out (we 
will do this explicitly below) and that the bosonic field expectation values do not 
have any spatial dependence, the effective action may be written as 

r[M/3a, M^a) ^^7^] = ^U[u,3a, M^a> ^7^], (4-3) 

where U serves as the effective potential for our theory. U is finite (apart from 
a T-dependent additive constant) and position independent. ( [4.2| ) tells us that 
— for vanishing sources — the system favors the state for which the effective 
potential as a function of the expectation values becomes stationary, and it can 



be shown that not only stationary, but even minimal |]T0|. 



To realize our mean field conditions, we introduce a factor 

n Siu^aiQ) - U^aSm S{UUQ) - ^k^iQ)) Siw.aiQ) " W,a5{Q)) (4.4) 
/37a 



under the functional integral of ( p. 47] ). This sets the bosonic fields to con- 
stant background fields, neglecting all bosonic fluctuations. The momentum 6- 
functions implement our condition that these constant background fields should 
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be homogeneous. Our partition function then becomes 

Zmf = exp (^-—^^^ J T^VT^i^ exp(-S'ir - 83- Sy - Sj) 
Sf = Y.'^KQ) i^^n - 2t (cos(gi/2)Ai + cos(g2/2)5i)) ^(Q) 



2 

7 



7a "^70 



abc Q 7 



^J = -ll[ J2^-JpaiO)upa + Jpa{0)u*f,J + J] L,,(0)«;,, + . (4.5) 

a \ 13 7 / 

Using the same approximations in ( ^?T| ) and inserting (|4.5|), by using ( [4.3| ) we 
arrive at 

^ = 5Z M ' 5Z + Y 5Z ^7a^7a + -^/^' + Af/ (4.6) 

with 

A?7 = In y r'^*r'V^exp(-^A) 

= 5^ ^^(Q) ii^l - It (cos(gi/2)Ai + cos(g2/2)5i)) ^{Q) 



Q 



- E E (E ^7cV^I(Q)Cc(Q' (4-7) 

afec Q 7 

+ [ulM{Q)vtXQ^ -Q)M-Q) + npc^^liQKZiQ, -QW,{-Q) 

As it must be, U can be written as the classical potential terms (which are 
pure mass terms in our theory) and a correction At/ describing the influence of 
fermionic fluctuations. We must now calculate the functional integral in At/. 
The easiest way to do so is to define the vector 
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Note that this vector has 16 components (2 explicitly, 4 in color and 2 in spinor 
space). By aid of this vector, we may rewrite S'a in the form 

^A = ^5^^'^(-g)p(gMg) (4.9) 

Q 

and perform the Gaussian integration to obtain 

A[/ = -^^lndetP(g). (4.10) 

Q 

In general, P{Q) contains contributions from all bosonic fields. From now on, we 
will set all fields except p, a and d equal to zero. For the remaining fields, we set 



Wpi = Wp2 = Wp3 =Wp4 = p 
Wsi = -WS2 = Wss = -WdA = a 
Udl = Ud2 = Ud3 =UdA = d 

<i = = = <4 = d*- (4-11) 

The explicit expression for P{Q) follows by inserting the vertex factors from 
appendix ^ in ([4.7|) and collecting terms. One obtains 

^(^)=W-2tT )-^^P\A, J-^"H^3®a 

- /id(cos(gi/2)Ai - cos(g2/2)fii) ® (^*^ ® la^ (4.12) 

with T = cos(gi/2)Ai + cos(g2/2)-Bi. By using the identity 

lndetA = ilndet(^A(^J j) (j j)) (4-13) 

it is possible to diagonalize P{Q) in the two dimensional space of the two com- 
ponents of ijj, so that the determinant reduces to one over a matrix of dimension 
8: 

AU = 5^1ndet8((cu^)2 + {2tf + hppf + hy 
Q 

+ 2hphapaA3 ®a + hld*d{cos{qi/2)Ai - cos(g2/2)5i)^) . (4.14) 
The determinant can be further diagonalized in spinor space, yielding 

aa — > |a| (T3. (4-15) 
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With {T, A3} = and defining 6 = d*d, a = (f, we find 

At/ = ^ lndet((cj^)2 + (2tr + V + /lav^As ® aa)^ 
Q 

+ hl5 (cos(gi/2) - cos(g2/2)5i)2) . 



(4.16) 



We proceed by evaluating the remaining determinant by brute force. In the 
resulting expression, the Matsubara sum can be performed using the product 
expansion 



13 



cosh(7r2;) = JJ^ ( 1 + 

neINo 



4^2 



(2n+ 1)^ 



(4.17) 



Up to a temperature dependent divergent constant, the final result for the effec- 
tive potential is 



U = 27rV' + 27r2a + A-kH + + At/ 



hi 



(4.18) 



AU = -2T 



In cosh + 62 V4t2(ci + eiC2)2 + hla^ + hl5{ci - eiCa)^ j 



with Cj = cos(gj/2). 



4.2 Discussion of the effective potential 

We will now discuss our result for the effective potential ( [4.18| ) and calculate the 
phase diagram. 

First note that for lar^e temperature AU vanishes and U is given hj U = 
27r^p^ + 27r^a; + 47r^5 + . In this case, the minimum of U with respect to 
a and 6 occurs at a = 5 = for all p. Therefore for large temperature, no 
symmetry breaking, i.e. no antiferromagnetism or superconducting behavior is 
present. If T is lowered, the minimum may be destabilized by the contribution of 
AU. However, note that (as can be seen by expanding AU for a ^ 1) U grows 
as a for large a, which means that even for low temperature the minimum always 
occurs at finite a. The same argument holds for S. 

A sufficient condition for the minimum to occur at non vanishing a or 5 is 
that the masses 

,.2 ^dU dU 

^" a=5=0 9^ 



(4.19) 
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become negative. If e.g. the mass becomes negative, the effective potential 
possesses a local maximum in a = 0. Since for large a the potential V increases 
oc a, the minimum must occur in a > and the symmetry is spontaneously 
broken. Again, the same argument holds for b. The masses can be calculated 
from ( |4.18 ) and we find 

d^q ^ tanh(2^(V + 2te2(ci + eiC2))) 



Ml = 47r2 - hi r /4 y 



, 2te2(ci + eiC2) 



- (2^)' ,,etl,i} ^(ci + eiC2) (cosh ^ + cosh ^i(£l+il£2) 
tanh f — (V + 2te2(ci + eiC2))J 

^ (^^-^^^^)^ ^ . (4.20) 
/ipP + 2te2(ci + eiC2) 



In the expression for we performed the sum over €2 to show that the right 
hand side is finite for ci + eiC2 and T > 0. We see that for 6, the mass 
correction arising from AU always tends to destabilize the symmetric minimum, 
since its contribution to the mass is always negative. The same holds true for a. 
Whether the mass actually becomes negative or not depends on the choice of the 
Yukawa couplings. By increasing the strength of the couplings, we necessarily find 
negative masses and therefore spontaneous symmetry breaking. It is clear from 
this qualitative point of view that our numerical results for the phase diagrams 
in this approximation will crucially depend on the initial choice of the couplings. 
In the frame of the mean field approximation, there is no way to decide which 
choice is the correct one (remember that without any approximations, any choice 
of the couplings respecting ( p.35| ) will be viable and physically equivalent). 

We stress that the condition of vanishing masses is sufficient, but not necessary 
for a phase transition to take place. To understand this, we consider the change 
of the effective potential with temperature as a function of a for fixed 6 (fig. 4.1). 
For high temperature, we know that the minimum occurs in a = 0. Suppose for 
low temperature the symmetry is broken. The phase transition into the state of 
broken symmetry can take place in two different ways: 



• The minimum of the potential remains at a = 0, until the mass M^ changes 
its sign. Then, as a function of temperature, the minimum moves away 
from a = to some finite a. The order parameter a for this kind of phase 
transition starts at zero in the symmetric phase and continuously moves 
away from zero in the broken phase. The phase transition is therefore 
continuous or of second order. 
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Figure 4.1: For high temperature, the potential minimum occurs in a = 0. 
By lowering the temperature, phase transitions of first or second order can ap- 
pear. Only for phase transitions of second order the temperature where the mass 
changes sign coincides with the temperature of the phase transition. 
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• At some temperature, the potential builds up a local minimum at some 
finite a. The value of the effective potential at this minimum decreases, 
until the minimum becomes global. remains positive during this pro- 
cess. In this case the order parameter a changes discontinuously from zero 
to some finite value, so that we face a discontinuous phase transition (or 
equivalently a phase transition of first order). 

This discussion shows that one has to be careful when analyzing the properties 
of the effective potential. For example, if we have found that say M| is negative 
and M\ is positive, we cannot immediately conclude that the system exhibits 
superconducting behavior, since it is possible that the global minimum occurs in 
a > 0, 5 = 0, if we happen to have the case of a first order phase transition in 
the antiferromagnetic channel. The main difficulties of numerically finding the 
minima of ( [4.18[ ) are rooted in these possibilities of first order phase transitions. 



However, even without explicitly calculating the minima of the effective poten- 
tial, it is possible to gain further insight into the possibility of first order phase 
transitions. First consider the case that the minimum occurs in a = 0, (5 > 0. 
Then this minimum must obey 

This equation may be rewritten in the form 

Ml = lhl£^^j:(^,-e,C2r (4.22) 

'tanh {^-^\/{hpp + 2te2(ci + eiCa))^ + hl5{ci - eiC2Y^ 

^{hpp + 2te2(ci + eiC2))2 + hl5{ci - e^c^f 
tanh {^{hpp + 2te2(ci + €102))) 

hpP + 2tt2{Ci + tiC2) 

By using 



tanh a tanh \/ a? + x"^ ^, , r. , a r.^\ 
> ,0^2 ^ 0' ^^-23) 

we see that the right hand side is strictly negative. Therefore solutions with 5 > 
are possible only for M| < 0, which in turn means that the phase transition is of 
second order. In other words: The phase transition between the symmetric phase 
and the superconducting phase is always of second order. The same argument 
holds for a, if p is not too large. For large p, we may (and will) encounter first 
order phase transitions into the antiferromagnetic phase. 
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Figure 4.2: The p-T phase diagram for /i^ = /i^ = /i^ = 10 with symmetric 
(SYM), antiferromagnetic (AF) and superconducting (SC) phase. In the region 
marked by the bold line the phase transition into the antiferromagnetic phase is 
of first order; all other phase transitions are of second order. 



The parameters of our theory are T, p, t and the Yukawa couplings. The 
chemical potential has been removed from the list of free parameters by the 
Legendre transform ([4.1|) and can be inferred from ( [4.2|) to be 



Particularly, we explicitly see by inserting ( [4.1^ ) that p = gives p = 0, so that 
p = corresponds to the case of half filling as it must. 



4.3 Numerical results 

In what follows, we set t = 1, fix the value of the Yukawa couplings and minimize 
U numerically with respect to a and 6 for a large number of (p, T)-pairs, which 
yields a complete picture of the phase diagram in this approximation. We present 
results for four different choices of the Yukawa couplings. 

Some remarks concerning the numerics: To evaluate the momentum integrals 
in the effective potential numerically, we used a very simple self written Rie- 
mannian sum like method on a grid. This works well for the kind of functions 
encountered here and is faster than more sophisticated methods. The search for 
the minimum is performed by using the Broyden-Fletcher-Goldfarb-Shanno vari- 
ant of the Davidon-Fletcher-Powell minimization. This method is implemented 
in the Numerical Recipes routine dfpmin (cf. 0]). To be sure that the method 
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0.1 0.2 0.3 0.4 0.5 

P 

Figure 4.3: The p-T phase diagram for h'^ — h"^ = h"^ = 40 with symmetric 
(SYM), antiferromagnetic (AF) and superconducting (SC) phase. In the region 
marlied by the bold line the phase transition into the antiferromagnetic phase is 
of first order; all other phase transitions are of second order. 
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Figure 4.4: The p-T phase diagram for /i^ = /i^ = 40, h"^ = 10 with symmetric 
(SYM) and superconducting (SC) phase. All phase transitions are of second 
order. 
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Figure 4.5: The p-T phase diagram for = = 40, — 10 with symmetric 
(SYM) and antiferromagnetic (AF) phase. In the region marked by the bold hne 
the phase transition into the antiferromagnetic phase is of first order; all other 
phase transitions are of second order. 




Figure 4.6: The p-T phase diagram for hp = hl = 40, hl = W. 
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Figure 4. 7: The p-T phase diagram for h"^ = h 



40. 



/)2 



10. 



does not return a local minimum, we repeat the calculation several times for each 
point of the phase diagram with random initial values of a and 6. To identify 
phase transitions of first order, we additionally calculate in every step and 
check whether for non vanishing minimum is positive or negative. 

The results for different choices of the Yukawa couplings are presented in the 
figures. Note that the choice of hp does not change the qualitative shape of the 
phases, since hp only enters as a factor for p and therefore does nothing else 
than to rescale the p-axis. For equal values of h^ and ha (figs. and |17^ ), 
we find phase diagrams that already resemble the phase diagram fig. 1.1 of a 
high temperature superconductor. It is interesting that no region exists where 
the minimum of the effective actions occurs for both a > and 6 > 0. An 
expectation value of the antiferromagnetic field tends to suppress spontaneous 
symmetry breaking in the superconducting channel and vice versa. Additionally, 
we find a first order phase transition into the antiferromagnetic phase for large 
p, which is in agreement with our discussion on analytical grounds in the last 
section. 

If we increase either the antiferromagnetic coupling ha or the superconducting 
coupling hd, the antiferromagnetic or superconducting phase dominates respec- 
tively. This feature is illustrated in figs. and ^75| . In fig. we plotted the 
value of S where the effective potential becomes minimal as a function of p and 
T. The choice for the Yukawa couplings is the same as in fig. |4.4| (/i^ = 40, 
hi = 10), so that only the superconducting phase is present. We did the same 
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in fig. |4.7| for a with Yukawa couplings as in fig. |4.5| , in which case symmetry 
breaking always takes place in the antiferromagnetic channel. This plot clearly 
shows the occurrence of a first order phase transition into the antiferromagnetic 
phase for large p. 

We conclude that even the mean field approximation presented here gives a 
rough picture of the phase diagram of high temperature superconductors. How- 
ever, the drawback is also apparent: Different choices of the Yukawa couplings 
lead to different phase diagrams, although the original theory was invariant un- 
der this choice. The strong dependence on the couplings is unphysical, and we 
are not able to remove this dependence in the framework of the mean field ap- 
proximation. To resolve this problem and to build up a more reliable picture of 
the phase diagram, we have to include bosonic fluctuations that we completely 
neglected in our mean field approach. However, the inclusion of bosonic fluctu- 
ations is highly non trivial and can no longer be treated by a simple calculation 
of an effective potential as we did here. The method of choice to deal with them 
is a renormalization group analysis of the effective potential and the couplings. 
In the next chapter we will present the renormalization group formalism that we 
will use in this work. After that, we further transform our partition function of 
the colored Hubbard model into a form more suitable for a renormalization group 
analysis. 
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Chapter 5 



Renormalization group equations 
and the effective average action 



One loop calculations in quantum field theory are usually plagued by divergen- 
cies, which appear as unbounded loop momentum integrals. Two sources of 
these divergencies are possible. For theories defined on a spatial continuum (as 
the usual theories defining the standard model), we face ultraviolet divergencies 
which occur because the integrals over the loop momenta extend to arbitrarily 
large momenta. This kind of divergencies is not present in our case, since the 
spacing of the underlying lattice provides a physical UV-cutoff for momentum 
integrals, constraining them to some finite interval (in our case to the inter- 
val [— 7r,7r]). Another source of divergencies is the presence of massless modes. 
These infrared divergencies do emerge in our theory on the lattice and have to 
be regularized. Whatever regularization scheme we use, this regularized theory 
will depend on some unphysical regularization parameter — an artificially intro- 
duced mass, a momentum cutoff, or whatever. The aim of any regularization 
procedure is to be able to calculate the loop integrals in the regularized theory 
and to remove the regulator afterwards to arrive at physical results. Renormal- 
ization group equations are a well established tool to describe the change of a 
given theory including quantum fluctuations with some flow parameter k, where 
k parameterizes the regulator. We will assume that k = corresponds to a van- 
ishing regulator. The ultimate goal then is to solve the flow equation, which gives 
the regularized theory as a function of k and to go to the limit /c — > yielding 
the physical theory including quantum fluctuations. 

In this chapter we describe the renormalization group formalism that we will 
use to further investigate the properties of the colored Hubbard model. How- 
ever, the treatment of the formalism presented here is rather general and applies 
to any system with fermionic and bosonic degrees of freedom. The concept of 
investigating the flow of the effective average action has been introduced in WE 



and has been used to treat a wide range of problems in quantum field theory and 
statistical physics (for a review, see 
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5.1 Generalized fields and regularization 



We consider a theory with a fermionic field ip*, a real bosonic field w and a 
complex bosonic field u, u* . The generalization to several fields is straightforward. 
We use generalized matrix notation by regarding the fields as vectors with a 
discrete label a (e.g. color, spin) and a position label X or momentum label Q. 
In the same way, we introduce generalized matrices. We then use generalized 
matrix multiplication e.g. in the form 



^t^^ = ^^^*(Q)A,^(Q,Q')^^(Q') 



(A5),^(gi, Q^) = Y^Y^ A„,(gi, Q)B^p{Q, Q,) 



(5.1) 



Similarly, Tr denotes a generalized trace. 
Now define the generalized fields 



u* 

X= ^ 
and the generahzed sources 



= {u\ iF, w\ i>\ i>^) 



K 



J* 
L 

V 



^1 



(5.2) 



(5.3) 



rj and rj* are fermionic sources (these sources enter in (|3.47| ), if we had bothered 
to explicitly write it out). Note that in general w ^ w*, although we introduced 
to as a real boson. But since our notation applies to both position and momentum 
space, we have to take care of the fact that w*{X) = w{X), but w*{Q) = w{—Q)\ 
The same applies to the real source L. 

The general partition function we want to consider is 



Z = 7V(T) / Vi)*Vi)Vu*VuDw exp(-S - S, 



0/1 



(5.4) 



where M{T) is some temperature dependent constant, S is the action without 
sources and 5*^ is a source term, which can be written in our matrix notation as 



Sj = —j'^u — J^u* — L'^w — r]'^il) — rfip* 



(5.5) 
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We regularize the theory by introducing a cutoff function AS'a;[x] so that 

S[x] ^ S[x] + ASkix]. (5.6) 

In the simple case of a theory with IR-divergencies only, the cutoff function could 
be a correction to the propagator terms in S, giving a fc-dependent mass to the 
troublesome massless mode. For k —>■ 0, we would then demand that this k- 
dependent mass vanishes, so that in this limit we recover the physical theory 
described by S and a possible source term only. Specifically, we define the cutoff 
function by 

ASk[x] = h:riR,xx^). (5.7) 

Rk is a cutoff matrix, and again we stress that this equation has to be read in the 
sense of generalized matrix notation. Note that the right hand side is not equal 
to |Tr (x^RkX) due to the fermionic fields contained in x. We assume that Rk is 
diagonal in the space of fields, so that we can set 

Rk = diag (i?^, (i?^)^, R^ Rl {Rlf). (5.8) 

R^ serves as a cutoff matrix for the complex bosonic field, R^ for the real bosonic 
field and R^ for the fermionic field. With this simplification, we can write ( |5.7| ) 
as 

ASk[x] = Tr (u^R'iu + Uj^R'^w + iJ^RU\ . (5.9) 



5.2 The effective average action 

The regularized partition function now becomes 

= M{T) j Vi;*V'tpVu*VuVw exp{-S[x] - ASk[x] + K^x) (5.10) 
and we define 

Wk = \nZk. (5.11) 
The (fc-dependent) expectation values of the fields are then given by 

XHX)^^§. X'HX*)^^. (5.12) 

We can now define 



fk[x] = -Wk + K\ 

Vk[x] = -Wk + K^X-^Su[x]. 



(5.13) 
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Note that these definitions are completely analogous to those we presented in 



(|4.1| ) except the last for rfc[x]. The reason why we subtracted AS'a;[x] the 
last of these definitions will become clear in a moment. For rfc[x], which is the 
Legendre transform of Wk, we have 

^=irtM, ^-^ = MK (5.14) 

OX ox^ 

with M = diag (1, 1, 1, -1, -1). 

Some remarks are in order concerning the interpretation of rfc[x]. This quan- 
tity is called the effective average action. We will now assume that 

lim Rk = 0, lim = oo, (5.15) 

k^O k^A 

where A is either a natural UV-cutoff of the theory (in our case something pro- 
portional to the inverse lattice spacing) or oo for theories with UV-divergencies. 
It is clear that then 

limTfeM =r[x], (5.16) 

where T[x] is the full effective action, since by letting R^ — >■ 0, we remove the 
cutoff function that causes to differ from T. In ( |4.1|) we already discussed the 
interpretation of the full effective action: The effective action yields — by use 
of the classical action principle — the equations of motion for the expectation 
values of the fields of the theory with all quantum fluctuations included. 

We now turn our attention to the limit k —>■ A. We start by rewriting ( |5.10 ) 
in the form 



exp(-rfc[x]) =Af{T) J Vil)*V'4)Vu*VuDw 

exp {^-K^x + A5,[x] - S[x] + K^x - ^Tr {RuXX^)^ , (5.17) 

where on the right hand side the sources are understood to be functions of x- 
Explicitly, we find by starting from (|5.14|) that 



= + x^RkM. (5.18) 

OX 



Inserting this in (|5.17|) and rearranging terms, we find 



exp(-rfc[x]) = A/'(T) J Vij*V^pVu*VuVw 

exp (^-S[x] + ^M{x - x)^ exp (^-i(x^ - x^)RkM{x - x)^ 



(5.19) 
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The second exponential vanishes for Rk ^ oo unless x = so that it effectively 
acts as a 5-f unction and we obtain 

lim exp{-Tk[x\) = const. ■ ^^{T) exp(-^[x]). (5.20) 

fc— >A 

Up to some irrelevant constants, we find that Tk[x] approaches the classical action 
S[x] with the fields replaced by their expectation values as A; — > A. Note that 
to achieve this nice property of Tk[x] it was necessary to subtract AS'a;[x] i^i the 
definition of rfc[x]. 

In conclusion, the effective average action interpolates between the known 
classical action and the unknown effective action. Neither the initial value F^^a 
nor the final value r^^o depends on how we chose to define Rk — it suffices that 
Rk meets the constraints ( p.l5| ). However, it is clear that the interpretation of 
Ffc for some finite k depends on the specific form of Rk- 

One very illuminating choice of the cutoff function is a simple sharp momen- 
tum cutoff Rk ~ k'^9{k'^ —q^). This ansatz meets the conditions ( |5.15|) for A oo. 
This kind of cutoff does not influence the momentum modes with > k"^, but 
gives a mass k"^ to the momentum modes with < fc^. The propagation of these 
low momentum modes is suppressed by the mass, so that in Tk quantum fluctua- 
tions of the high momentum modes only are integrated out, yielding an effective 
theory at scale k for the propagation of the low momentum modes. In position 
space, this means that by lowering k we average over larger and larger regions, 
integrating out the short range fluctuations and building up a an effective theory 
for long range fluctuations only. 

It is a well known feature of the effective action that it preserves the symme- 
tries of the original action in the sense that the same symmetry operations acting 
on the fields in the original action applied to the field expectation values leave 
the effective action invariant, if these symmetry transformations are linear. This 
is the case for all the symmetries we consider here, so that we expect the effective 
action to respect the same symmetries we discussed before for the action of the 
Hubbard model. However, whether Tk[x] respects these symmetries depends on 
the choice of AS'fc[x]- It is convenient to choose AS'^fx] such that it also respects 
the symmetries of the theory. Then Tk[x\ is invariant under the symmetry trans- 
formations for all k. The advantage of Tk being invariant under the symmetry 
transformations is that we will have to write down an ansatz for Tk to solve the 
flow equation calculated in the next section. During the flow an infinite number 
of terms contributing to Tk will be generated, so that we need some guidance 
in selecting terms that we include in our truncation. The fact that only terms 
respecting the symmetries of the original theory are allowed in constrains the 
number of terms we may or may not include in our ansatz, which makes the 
motivation of truncation schemes much easier. Another reason for demanding 
invariance of Tk is that we are often interested in symmetry breaking properties 
of the theory at some scale k and not only in the limit k ^ 0. Of course, if 
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Tk breaks the symmetries at finite k explicitly, spontaneous symmetry breaking 
properties are completely obscured. 



5.3 The flow equation for the eff^ective average 
action 

In this section we derive an exact flow equation for Tk\x] that can be used to 
calculate the effective average action. 
We start by noting that 

^=(«'>-<x)(xt>, (5.21) 
which is the connected 2-point function, and 



Recall that in our notation, expressions like jj^^ are matrices. The derivation 
of the flow equation is now straightforward: 



so that 

-|r. = iTv((WM(|^ + (5.24) 

This is our master equation for deriving all the flow equations we need. The 
ultimate goal is to solve this equation with an initial condition given by the 
classical action. In general it is not possible to find an exact solution to the flow 
equation. However, the fact that we know that F/^. equals the classical action at 
the beginning of the flow k ^ A allows to motivate sensible truncation schemes 
for possible solutions. 

We have to stress that up to this point no approximations entered the calcu- 
lation, so that the result is exact. This is particularly interesting if one compares 
it to the perturbative one loop calculation of the effective action, which yields 

|r, = iTv((W/(gM + (5.26) 
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This is completely the same, except that the exact propagator is replaced by the 
classical propagator here. The simple replacement of the classical propagator by 
the exact propagator in the one loop correction to the effective average action 
renders this equation exact, including arbitrary high loop orders and genuinely 
non perturbative effects. The formal similarity to perturbative one loop expres- 
sions allows to use well known calculation techniques and leads to results that 
can be simply interpreted in a diagrammatic language. For further discussion. 



applications and references, see [16 
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Chapter 6 

Diagonalization of the 
propagator matrix 



At this point, we are in principle done with the preparations to extract the prop- 
erties of the Hubbard modeL We have derived a partition function (p.47|) , which 



contains the interesting degrees of freedom in an exphcit way, and we have pre- 
sented the formahsm of renormahzation group equations for the effective average 
action, which we can use to derive the effective action. However, as we aheady 



stated, ( |5.24| ) cannot be solved exactly. Instead, we have to invent some trun- 
cation for the effective average action, which consists of terms respecting the 
symmetries of our theory containing /c-dependent variables (masses, couplings, 
wave function renormahzation constants etc.). Since our theory is non renor- 
malizable, it is clear that one cannot avoid approximations at this point. The 
question is: Which terms should be included in our truncation without making 
the truncation error too large? As we mentioned in the course of discussing sym- 
metries of Ffc, it is clear that the smaller the number of terms allowed by exact 
symmetries, the smaller the truncation error that we make by discarding terms. 

Unfortunately, in the present formulation of our theory (|3.47|) , a huge number 
of terms in the effective action are allowed, all of which are of the same order of 
magnitude. Particularly, consider the propagator terms of the form 

J2'>^UQ)PmAQ,Q>MQ'), (6-1) 

QQ' 

Y,Wia{Q)P^-y'AQ^Q>i'b{Q')- (6.2) 
QQ' 

We know that — due to f/(l)-invariance — propagator terms mixing complex 
and real bosons will not occur in the effective average action. Similarly, the spin 
singlet bosons and spin triplet bosons will not mix due to SU (2)-invariance. To 
keep things simple, we will set Ai 0, which because of ( |3.35|) is equivalent to 
discarding the bosons s, c and tx/y from our theory. Then we have three sets of 
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bosons (real bosons in the spin singlet, real bosons in the spin triplet and complex 
bosons) with four boson species each. Every boson species occurs in 4 different 
colors. This means that we face 3(4-4)^ = 768 boson propagator matrix entries, 
non of which vanish due to symmetries. Of course, it is completely hopeless to 
include so many terms in any useful truncation scheme, but on the other hand, 
by considering only some of the terms and neglecting all others, we introduce 
large truncation errors. 

The problem that we face here is rooted in the fact that we are not able to 
make use of the additional lattice symmetries (translations, rotations and reflec- 
tions) to narrow down the number of allowed terms, since the bosons mix under 
these transformations. If we were given a set of bosons belonging to different 
inequivalent representations of these lattice symmetries, it would be clear that 
mixing between these different bosons could not happen. In this case, most of 
the terms of the propagator matrix vanish due to exact symmetries. The idea 
of this chapter is to find linear combinations of the existing bosons to make the 
corresponding states eigenstates of the symmetry transformations. Truncation 
schemes are then proposed for a theory which is written in terms of these new 
bosons. 

The simplest way we found to attack this problem is to explicitly calculate the 
oneloop corrections to the bosonic propagators perturbatively and to diagonalize 
the resulting one loop improved propagator matrix. As a byproduct, we gain 
some information about the propagation of bosonic modes that will be useful 
when we define truncation schemes. Although we diagonalized the perturbative 
expressions only, it will turn out that the propagator matrix remains diagonal 
to all orders in perturbation theory and that in fact this transcription leads to 
bosonic states which are eigenstates of symmetries of our theory. Finally, we 
will be able to write down a partition function of our model with a new set of 
bosons that all belong to different inequivalent representations of translational 
symmetries and therefore do not mix. Much of the material in this chapter has 



been published in [12 



6.1 The diagonalization procedure 



6.1.1 The oneloop calculation 



Our starting point is ( 3.47| ). We want to calculate the one loop corrections to the 
bosonic propagators. It is clear that these one loop corrections correspond to the 
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diagrams 

K + Q 




Q 




where the sohd hnes denote fermions and the dashed hnes bosons. Thus only 
fermionic fluctuations enter these one loop corrections. To calculate these, we 
expand the fields around their expectation values 

+ r{Q)^r{Q)+5r{Q) (6.3) 

and find for the action S = Sp + Sb + Sy + Sj 

S\u, u*, w, -i/^*] = S[u, u*, w, ip, + S[u, u*, w, 5ip, Stp*] . (6.4) 

^ V ' ^ V ' 

=5o =52 

The effective action (with bosonic fiuctuations neglected) is then given by 

T[u,u*,w,ip,ij*] = + So - Sj - \n J VSij*VSijexp{-S2). (6.5) 

Note the similarity to the meanfield calculation of chapter ^. The only difference 
is that our bosonic background fields are not assumed to be homogeneous, and 
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we will only consider corrections to the propagator (which, in the mean field case, 
would correspond to expanding the effective potential up to quadratic order in 
the fields). A similar transcription as in chapter ^ 

m - (^) (6.6) 

allows to perform the functional integral in (|6.5|) and we get 

T[u,u*,w,ip,ip*] = + 5*0 - S J -^Indet P[u,u*,w] (6.7) 

P ^ V ' 

=Ar 



with 



P[u, u\ w] (Q, Q') = ^^^^^^ - (6.8) 

^ v ' ^ V ' 

=Po =AP 

P^iQ, Q') = i^oon - 2t(cos(gi/2)Ai + cos{q2/2)Bi))6{Q - Q') (6.9) 

A = Y,Y.'^:'^Q,Q>,c{Q-Q') 

7 c 

B = '2Y.Y. Q'>Pc{Q + Q') (610) 

/3 c 
f3 c 



Recall that the vertex factors V^c are matrices in color and spinor space with 
components Vah,c in color space. 

We can now expand AF in numbers of bosonic fields 

Ar = -^Tr InP = -^Tr ln(Po(l - Pq^^P)) 

= -i ^Tr InPo - MPo'^P) " ^Tr (PcT^^P)' + ■■•)■ (6-11) 



Only the third term of this expansion contributes to the propagator corrections 
we want to calculate. Note that the trace involves summation over color, spin 
and generalized momentum. 
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More explicitly, we have for AFa = l/4Tr(Po"^AP)2 

= $^[$^$^^^,.(-K)Tr{(p^)-i(g)vr(Q,ir + g)(p'^)-i(ir + g) 

KQ 77' cc' 

v7(K + g,g)}w7y,,(ir) 

-2 E E^^cW Tr{(p^)-^(g)Kr''(g, - Q){{p'')-'f{K - g) 

/3/3' cc' 

vf{K-Q,Q)yp,,{K)\. (6.12) 



Note that we have written out the momentum sums, so that the trace only sums 
over color and spin. 

From ( |5.22| ) we know that the inverse propagator is simply given by the sec- 
ond derivative of F with respect to the field expectation values. The one loop 
corrections to the inverse propagators are therefore given by 



af(2) (K) 

= J] Tr { iP^r\Q)Vp{Q, K + Q){P^r\K + Q)V;r' {K + Q,Q)} 
Q 

^^^lu.AK) (6.13) 

= -2 5^Tr {(P'^)'^(g)v:c^(g, k - q){{p^)-Y{k - Q)v^}'{K - g, g) 

Q 

where ATw!^^wy^,{K) denotes the second derivative of AF2 with respect to the 
field expectation values. 

The one loop calculations will be simplified if we make a slight change to our 
bosonization prescription. We replace the fermion bilinears e and d by 



ey = l{i + d), ~e. = \(e-d) (6.14) 



The bosonization procedure gives the same result as before with the following 
modifications: Replace e, h^, by ty, h^.^ and d, by e^, h^^ everywhere. In ( |3.35| ) 



we have to replace 2He = 2Hd = H^^^^ = 6A3 by H^^^^ = H^^^^ = 6A3. The 
vertex factors are given by V^y = 1/2{V^ + V^) and = 1/2{V^ - V^). 

In principle one could now insert the vertex factors and fermionic propagators 
and proceed by calculating the traces for all possible 768 propagator matrix 
entries (we actually did that, but it is not very illuminating to present it here). 
However, a simple transformation renders the propagator matrix diagonal in color 
space. We will discuss this transformation in the following section and present 
the results for the one loop corrections afterwards. 
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6.1.2 Diagonalization in color space 

The idea to diagonahze the propagator matrix in color space is to note that for 
homogeneous fields, translations can be easily applied by color transformations. 
Particularly, we have 

Txba = ^{Ai)abbb, Tyba = ^iBi)abbb, (6.15) 
b b 

where ba is any of our bosons with color a. It is clear that Tj^^ = 1 in this case and 
Ai and Bi commute. The group of transformations consisting of 1, Tj., T^^ and 
T^Ty is therefore isomorphic to G = Z2X Z2, where Z2 is the cyclic group of order 
2. G is Abelian and possesses therefore 4 classes. The number of inequivalent 
irreducible representations of any group is equal to the number of classes of the 
group, so that we have four inequivalent irreducible representations, which are 
necessarily one dimensional. A basis for these irreducible representations is easy 
to define. Simply take 

1 - 1 

bi = -{bi + 62 + 63 + 64), &2 = -^{bi - 62 + 63 - 64), 

63 = 1(61 + 62-63-64), h = \{b,-b2-bs + b,). (6.16) 

The representations of the elements of G = {l,Tx,Ty,TrcTy} in this basis are 
{1, 1, 1, 1} for 61, {1, -1, -1, 1} for 62, {1, 1, -1, -1} for 63 and {1, -1, 1, -1} for 
64. Since we demand all terms in the effective action to respect the translation 
symmetries, it follows that all propagator matrix entries mixing different colors 
have to vanish. For example, 

blP~^fH2 = (Tj*)P'^i'^2(T^62) = -blP^*H2 = 0, (6.17) 

if pKb2 ig ^Jig propagator matrix element coupling the bosons bl and 62. 

If we write ( |6.16|) in the form ba = \ Matbb, we have the inversion 6^ = 
'Y^i^Mabbb- Inserting this in ( 3.47|) , we find that as functionals of the new bosons 



6, the bosonic terms become 

a Q \ p 7 

Sy = -Y.Y. KQ-Q'-Q") 

abc QQ'Q" 



W-ya{~Q)w^a{Q) 



E fe(Q)V^r(Q')Kt(Q', Q")UQ") + S/3c(Q)^I(Q')Cc(^'' Q")rb{Q"] 







+ 5^^,,(Q)V^t(g')v;-;^(g',_g'/)4(_g")^, (6.18) 
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where the new vertex factors are given by V'^ = MahV\. 

The above argument only apphes to homogeneous fields. In general, the 
propagator matrix elements for different colors will not vanish for non vanishing 
momentum. However, by choosing the Fourier transforms for the original bosonic 
fields as in ( |A.2|) , we can show — by explicitly calculating the propagator matrix 
entries — that indeed the transition to the fields h renders the propagator matrix 
diagonal in color space for arbitrary bosonic momentum. 



6.1.3 Color diagonal one loop results 

We present our results for the one loop corrections to the bosonic propagators of 
the new bosons h which can be derived from ( |6.13|) . 

Consider the 12 sets {p,p,qy,qx}a, {a,m,gy,g^}a and {ey,ex,Vy,Vx}a, where 
a = 1, . . . , 4 is the color index. All propagator matrix elements connecting bosons 
belonging to different sets vanish (due to U{1)- and 5'[/(2)-invariance and diago- 
nalization in color space). The full propagator matrix is then block diagonal with 
the blocks given by 12 4 x 4 matrices. We call the one loop corrections to the 
4 X 4-propagator matrix blocks corresponding to the four colors of {p,p, %, qx}a 
(Ar^)a, (Ar^)„ for {a,m,gy,g^}a and (Ar^)^ for {ey,ex,Vy,Vx}a- 

For the sets of real bosons in the spin singlet, we find for {Ar^)^a\K) , K = 

(Ar^)f (ir) 



= ^ r T^^'^^i' +) + 9Li- -)), 

= -T 1^ + 2vr, h + 2vr)(^?;(+, +) + gl{-, -)), 

= T r ^,S^{h, h + 2n){gL{+, -) + gL{-, +)), 
= -T f ^,S''\h + 271, k,){g\{+, -) + gl{-, +)). (6.19) 



>S"'(/ci, ^2) is a 4 X 4- matrix given by 



hphpSiS2 
2 



hphpSiS2 
4/ipSnfsn2 



2ihphgyCsfs2 2ihphgySisnl 
\-2ihphgJiCsl -2ihphg^snjs2 



-2ihphqyCsls2 2ihphqjicsl\ 

-2ihphqy sisnl 2ihphq^sn\s2 

4^g,CS?Snf -hqyhqJlS2 

-hqyKJi^2 4:hl^snlcsl ) 



(6.20) 



with csj = cos(A;j/4), suj = sin(A;j/4) and Sj 



= sin(A;j/2). The functions g are 
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Matsubara sums {m,n e Z, e {1, —1}): 

rc( ) ^^^^ ^^^^^ ^^^2) + ^S^^' 

9e's [^u e2) - 2^ ^^^^ ^ ^ ^ ^ 

which depend on m, k and q. The upper sign apphes to real bosons (r) and the 
lower sign for complex bosons (c) (wc will use it later when we write down the 
one loop results for the complex bosons). The sums Si and 5*2 are 

Si (m, a, h) := g [(2^ + 1)2 + a^] [{2{n + m) + If + b^] 



TT 6(4m2 - + b^) tanh(f ) + a{Am'^ + - b^) tanh(f ) 
2 a6[4m2 + (a + b)^][Am^ + (a - bf] ' 

(2n+ l)(2(n + m) + 1) 



(6.22) 



,^2 (m, a, b) := ^ 



[(2n + 1)2 + a2] [(2(n + m) + 1)2 + b^] 



(6.23) 



_ TT a(4m2 + - fe^) tanh(f ) + 6(4m2 - + b"^) tanh( 
~ 2 [4m2 + (a + 6)2] [4m2 + (a - 6)2] 

The frequencies a;, a;' appearing in the definition of g read 

uj^{2n+ DttT a;' = I (2("^ + "-) + '^)t^T for real bosons ^4) 
^ ' [ (2(m — n) — l)7rT for complex bosons. ^ ' ' 

For the arguments of 5'i_2 we use the abbreviations (e^ e {1, —1}) 

ae, = (ci + eiC2)/(7rT) 

6,^ = (ci + e,4)/(7rT), (6.25) 

where q and c[ are given by 

, , f 2i cos((/ci + Oi)/2) for real bosons 

Cj = 2tcos(gi/2), c = <^ );/ ^*(',^( , , , (6.26) 

' " ' \ 2tcos{{ki- qi)/2) for complex bosons. ^ ' 

For the real bosons in the spin triplet, we get the same results with hp ^ hm, 
hp^ha and V/, ^ V/,- 

The result for the complex bosons is similar. Define 

/ -2hl^cl_cl^ \heyhe^CiC2 iheyKyC\_S2+ -\ih^yKjiC2\ 

S^{q q') = i^e,/'e,SiC2 -2hl^cl^4_ -\ih^^^K^CiS2 ihe^h^A+ci^ 

' —iheyhvyCi_S2+ 2'^he^hyyCiS2 — 2/l^^C^_S2+ 2^V:^hyySiS2 

\ l'iheyhvJiC2 -ihe^K^sl^4_ \K^hyySiS2 -2hljl^cl_ J 

(6.27) 
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Figure 6.1: The second derivative of the one loop correction to the bosonic 
kinetic term in the effective action of 0.2 for high temperature T = 1. We have 
set t = 1; the Yukawa couphngs are at their Hubbard model values ( 6.46|) with 
hi = 10. We have plotted the Matsubara mode m = 0. 

with Ci± = cos((gi ± g-)/4), Si± = sin((gi ± g-)/4), q = cosgj + cosg- and Si = 
siuQi + sing^. Then we have for {AT^)^\k) 

{Ar^)f\K) =tI^ {S%q, k - q)g'Li-, -) + A^S^q, k - q)A,gl{+, +)) 

(Ar^)f (K) = -T 1^ ^2^S\q + 27r(ei + e^), k-q- 27r(ei + e^M+i-, -) 
+ A^S%q + 27r(ei + 62), k-q- 27r(ei + e2))A^gX{+, +)) 

{AT^)t\K) =T r ^,{S%q + 27762, k-q- 2Tre,)g^L{+, -) 
+ AaS^iq + 27re2, k-q- 27re2)A3^^(-, +)) 

{Ar^)f\K) = -T 1^ ^2^S\q + 27rei, k-q- 2T,e,)gl{+, -) 

+ A^S'{q + 27rei, k-q- 2'Kei)A:,gl{-, +)). (6.28) 

6.1.4 Discussion of the one loop expressions 

As a byproduct of the diagonahzation procedure of the propagator matrix we have 
at hand the exphcit one loop corrections to the bosonic kinetic terms. We here 
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Figure 6.2: The second derivative of the one loop correction to the bosonic 
kinetic term in the effective action of 0,2 for low temperature T = 0.01. We have 
set t = 1; the Yukawa couphngs are at their Hubbard model values ( |6.46 ) with 
hi = 10. We have plotted the Matsubara mode m = 0. 

discuss the generahzed momentum dependence of the one loop expressions given 
by ( |6.19| ) and ( |6.28| ). As the generic case, we show plots of these expressions as 



functions of the bosonic momentum for high and low temperature for the boson 
02. The most remarkable features of the one loop expressions are the following: 

• The one loop expressions are periodic in 4n, not in 2tt as one would naively 
expect from the boundaries of the integrations. In fact, even if one takes 
into account the periodicity properties of the bosonic fields, it turns out 
that no single bosonic propagator term in the effective action is periodic in 
2n. This leads to interesting consequences discussed in the last section of 
this chapter. 

• The plots show that the one loop correction becomes minimal at = 0, 
which holds true for any temperature. Since the classical propagators are 
momentum independent, the momentum dependence of the one loop cor- 
rected propagator is given by the one loop correction term and our result 
means that the propagation of momentum modes near vanishing momen- 
tum is facilitated in comparison to the higher momentum modes. The zero 
momentum mode corresponds to a spatially homogeneous field. We con- 
clude that an approximation of the effective action that mainly keeps the 
dependence on spatially homogeneous fields will be justified, as the dynam- 
ics of the system is dominated by these homogeneous fields. 
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• For high temperature the momentum dependence is quite simple. It is 
mainly given by — (cos(fci/4) cos(fc2/4))^ in the case of 02. For low tem- 
perature, the momentum dependence becomes more complicated. Note the 
appearance of the cross shaped region in fig. ( |6.2| ), where the momen- 
tum dependence becomes non analytical for low temperature. This cross 
exactly corresponds to the Fermi surface in the colored Hubbard model, 
which makes perfect sense. 

The high temperature limit of the one loop expression can be easily calculated 
analytically. To do this, we can expand the one loop expressions with respect to 
1/T. Using 

9Uei,e2) = j^Smo + 0{l/T'), g^Je^e^) = -j^.S^o + 0{1/T') (6.29) 

we see that in this limit the momentum dependence of the propagators is com- 
pletely dominated by the m = 0-Matsubara mode. For low temperature, where 
we expand the sums Si and 5*2 with respect to T, we note that m enters only 
via the product {mT)'^, so that the first m-dependent term in the expansion with 
respect to T is oc m^. 



6.1.5 The final diagonalization step 

Although we have greatly reduced the number of non vanishing propagator matrix 
entries, we can do even better by diagonalizing the remaining 4 x4-matrices ( |6.20| ) 
and ( |6.27|) . Define the transformation matrices 



/ipSnisn2 
z/igyCSisn2 
y-2/;.q^snicn2 



-/ipsnisn2 ihg^sn2 -i/ig^snA 



/;,pCSiCS2 





and 



Then the transformations 



e,CS2 





\ 










hpCS2 





iKySn2 


heyCS2 








hpCSi ) 



\ 

z/itj^sni 


/ie.CSi / 



(6.30) 



(6.31) 



Vl{ki,k2) 



U^kiM) 

W{ki,k2) 



U;{ki,k2)=U%ki,k2) 



Ul{h,k2) 



U'{ki,k2) 



|csi— »— isni 
CS2— >jsn2 
srii— >i csi 
sn2--*-ics2 

|cs2^i sn2 
sn2^— ics2 

|csi — ^— zsni 
sni— >j csi 



U'Aki,k2) 
U^ih,k2) 

Ul{ki,k2) 



U%ki,k2) 
U\ki,k2) 

U%ki,k2) 
U\ki,k2) 



-i cSi 



|CS2— ►«sn2 
sn2— >— icS2 



I csi— »i sni 
sni— >— i csi 



(6.32) 
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render the one loop propagator matrices diagonal for the real bosons and block 
diagonal for the complex bosons. More explicitly, we have 

S'ih, k2) = Uli-h, -k2fS'ih, k2)Ul{h, k2) = 

diag (^4 (cslihlcsl + hlsnl) + snl{hlsnl + /ij.csa))' , 0, 0, 0^ , 

and 

U^{-ki, -k2fS'{ki + 27r, k2 + 2n)U^{ki, k2) = S^k^ + 27r, k2 + 27r) 

Ul{-ki, -k2fS''{k,, k2 + 2n)U;{k,, k2) = S'iki, k2 + 27r) 

Uli-k,, -k2fS'{k^ + 27r, k2)Ul{k^, k2) = S'{ki + 27r, k2) (6.33) 

for the real bosons and 

S'^iq, k-q)^ U^iki, k2)^S%q, k - q^k,, A^s) = 

^ -2cl ( hi c4 + hi snlY 2c,^C2-{hlcsl + hlsni) 

^ «cs? + <snf) 

2c^^C2-{hlcsl + hlsni] ^^2 ^^2 ,^2 ^^2.2 q 

/'l,2 2,7 2 2\ 2- y Cx L ' Vx 1/ 



V y 

^72(^1, hyS'^iq + 27r(ei + 62), fc - g - 27r(ei + e2))U;iki, ^2) = 

+ 27r(ei + 62), fc - g - 27r(ei + 62)) 

C/|(A;i, A;2)^5^(qr + 27re2, k-q- 2T:e2)Ul{k^, ^2) = 

S\q + 2T:e2,k - q - 2T:e2) 

U^iki, k2yS%q + 27rei, k-q- 27rei)Ul{ki, k2) = 

^% + 27rei,fe-qf-27rei) (6.34) 

for the complex bosons. 

The next step is to define the set of bosons corresponding to the transforma- 
tions we performed in order to diagonalize the propagator matrix and to express 
the partition function with respect to these new bosons. Define the new bosons 
Ra-, Sa and Xa by 



(p, P, qy, q.)jK)^U:iK)R,{K) 
(e, d, vy, %)l{K) = Ul{K)xa{K). (6.35) 
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Here the new bosons Ra, Sa and Xa are vectors with four components for each 
color a. Wc have written these transformation rules for the expectation values, 
but the same applies for the fields Ra, Sa and Xa- Then the one loop correction 
to the propagator matrix for these new bosons has the form calculated above. In 
analogy to the color transformation in the last section, the vertex factors of the 
theory with the new bosons are given by linear combinations of the old vertex 
factors 

T/^" ^ VP 

^2 _ (TTr\T 

yRa - Wa) ygy 

and similarly for Sa and Xa- Again, be aware of the notation: For any fixed color 

a, the old vertex factors V^, etc. are 4 x 4-matrices in color space. The different 
vertex factors for p, p and qy/^ are linearly combined by the matrix {Uiy to give 
the new vertex factors V^" , which are again 4 x 4-matrices. The vertex factors 
can be calculated explicitly, and it turns out that 

V^-^ = ^3^" = » = y/" = Vi'^ = Vt = V^a^" = = Va. (6.37) 

This leads to a great simplification of our formalism, since all the bosons with 
vanishing vertex factors decouple from the fermionic sector of the theory and can 
be integrated out! The theory depends only on the remaining set of bosons, i.e. 
on four real bosons in the spin singlet, four real bosons in the spin triplet and 
eight complex bosons. For each vector Ra with fixed a, one boson remains that 
we again call Ra, but Ra now understood to represent exactly one boson. In the 
same way, Sa will denote one boson for each a, and Xa denotes a vector with 
two bosons for each a. We will present the complete expression for the partition 
function with the new bosons in the following section. 

6.2 The final form of the partition function 

In the last section we diagonalized the propagator matrix in two steps, first in 
color space, and then in the spaces of boson species blocks for each color. We 
now present our final form of the partition function that we will actually use for 
renormalization group calculations. Since the stepwise diagonalization procedure 
is very error prone, we show that by starting with our final form of the partition 
function, we can reproduce the Hubbard model by integrating out the bosons 
and inserting the Hubbard model values for the Yukawa couplings. The rest of 
this section will discuss the partition function and the one loop expressions. 



(6.36) 
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6.2.1 The partition function 

Our final version of the partition function is 



Z = N{T) j DrDi;DrDxDRDse^v{-{S^in + SL + Sy + S,)] . 



(6.38) 



The real fields R, s as well as the complex fields 



X 



X 

* I ) 



X 



Xx 

Xy 



(6.39) 



are understood to carry a color- index a — 1 ... 4. The terms of the action read 



Q ah 



K ab 



-U-K)Pi{K)s,{K) 



+ xi{K)Pl{K)UK) 



+ Sc 



.ml{Q)vUK)MQ')) 
5{K-Q- Q') (xl{ml{Q)vX.{Q. Q')MQ') 

+ximi{Q)viMQ'mQ') 

-^^^ = - EE [Lf{-Q)Rc{Q) + Ll{-Q)l{Q) + Lf\Q)xc{Q) + Lf{Q)xl{Q) 



KQQ' abc 



+ 



Q c 



+ v:{Q)MQ)+Vc{Q)rc{Q) 



The propagator matrices are {B stands for i?, s or x) 

Pi{Q) = [lUn - 2t (cos(gi/2)^i + cos(g2/2)5i)]„, , 

PabiK) = {27:fP!{K)5ab 



(6.40) 



(6.41) 



with 



P^{K) = hlcos^{ki/A) cos2(A;2/4) + hlsiv?{ki/A) sm^{k2/A) 

+ cos^(A;i/4) sin^(A;2/4) + h]^ sin2(A;i/4) cos^(A;2/4) 

P({K) = hi cos\ki/A) cos\k2/A) + hi sm\ki/A) sin^k^/A) 

+ hi cos^(A;i/4) sin2(A;2/4) + hi sin^iki/A) cos^{k2/4) 



(6.42) 



6.2. The final form of the partition function 



83 



2 _.2n. iA\ , u2 „,^2,u iA\ (6-43) 



for the real bosons and 

P^{K) = hi co^\h/A) + sin^(fci/4) 
= hi coslh/A) + hi s\n\h/A) 

for the complex bosons. P^{K), P^{K) and P-f{K) for a = 2 ... 4 are given by 
the expressions for a = 1 with the replacements sin(A;j/4) ^ cos(A;i/4) for a = 2, 
sin(A;2/4) ^ cos(fc2/4) for a = 3 and sin(A;i/4) cos(A;i/4) for a = 4. 
The vertices for the real bosons read 

V^K) = P^{K){D,U ® o-r", Vl,{K) = P^{K){D,U ® ^r", (6.44) 

For the complex bosons, the vertices are 

^a6.iW,yj - 2 (^cos((g2-g^)/4)Pr(g + Q')(i?i)aJ ^ 

_ 1 ^sin((gi - gl)/4)P|(g + Q'){A,),,\ 
K,,2{Q.Q)- 2 ^sin((g2-g^)/4)P2^(g + g')(i?2)aJ ® 

yx* (n n') - ^ - ?0/4)P|(g + g')Pi)aA „ ..^ ^ (6.45) 

VsW^gJ - 2 1^ sin((g2-g^)/4)P3ng + g')(^^i)a6 J ® ^'^2^' 

yx* (o m - 1 ^-sm((gi - gO/4)Pi'(g + Q'mu\ 



2 V- COs((g2 - g2)/4)Pr(g + Q'){C2)ab^ 

Kyg,g') = -K€(-g,-g')- 

The matrices A^, P^, and P^ are defined in the appendix ( |A.2| ). 



6.2.2 Equivalence to the Hubbard model 



In this section we will show that our new partition function ( 6.38|) is equivalent to 



the one of the Hubbard model, if we set the Yukawa couplings to their Hubbard 
model values according to 



hi = '^H,U 

Hp = 3(A2 - A3) Prn = A2 + 3A3 + 1 

Hp = 3(A2 + A3) H^ = A2 - 3A3 + 1 

= Hgy = 3A2 Pg, = P^^ = A2 + 1 

He:, = Hey = Hy^ = Hy^ = 6 A3. 



(6.46) 
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In this case, the propagators simphfy to 

P/^(i^) = /^(ir, 1, 1) Pf = f{K, 1, 1) 

p,''{K) = f\K, -1,1) Pi =nK, -1,1) 

P^{K) = f iK, 1, -1) Pi = fiK, 1, -1) (6.47) 
P^iK) = f {K, -1, -1)P( = f{K, -1, -1) 

with 



(6.48) 



/^(ir, ei, es) = TT^ (^A2 - ei Y ( cos(A;i/2) + cos(fc2/2) 
fiK, e„ e^) = TT^ Q + ^ + ( cos(fci/2) + cos{h/2) 
The solution of the field equation for Rc{K) is 

RciK) = ^-^^^ E E ^(^ + ^ - Q')^mylc{K)MQ') + i^^^- 

(6.49) 

Inserting this in the action, we see that only the source independent part of this 
solution contributes to the four fermion term in the purely fermionic theory. We 
will therefore set the sources equal to zero from now on (it is a simple task to 
check that the source dependent terms produce the correct source dependent 
terms of (t3.28|)). The contribution from the i?-bosons to the four fermion term 



in the action is 

Sa=\ E ^(Q - Q' + - Q') ( - A2([1111] + [2222] + [3333] + [4444]) 

QQ'QQ' 

+ As cos((gi - g0/2)([1122] + [3344]) + A3 cos((g2 - g0/2)([1144] + [2233])) 

(6.50) 

where we use the short hand notation 

[ahcd] = ^Pi{Q)MQ')^l{Q)MQ')- (6.51) 

We proceed similarly for the s-bosons. To bring 5*1 into a form which can easily 
be compared with ( |6.5CI| ), we use the identity 

^i{Q)aMQ')^liQWMQ') 

= -^iiQ)MQ')^liQ)MQ') - 2i;i{Q)MQ')^liQ)MQ') (6-52) 

= —[abcd\ — 2 [a(ic6] 24. 
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Here [abed] is defined as in ( |6.51| ) and [abcd]ij is equal to [abed] with the momenta 
of the i-th and j-th field exchanged. We find 

S!=^ J2 5(g-g' + g-g')((l + A2)([llll] + [2222] + [3333] + [4444]) 

QQ'QQ' 

+ A3Cos((gi - g;)/2)([1122] + [3344] + 2([1221]24 + [3443]24)) 



+ A3 cos((g2 - g2)/2)([1144] + [2233] + 2([1441]24 + [2332]24)) 



(6.53) 



(6.54) 



For the x-bosons, we use the identity 

= -^;(Q')v^6(g')^KQ)V'a(Q) -^l(g)4(Q')^;(Q')4(Q) 

to calculate 

^4^ = ^ E m-Q' + Q-Q') 

QQ'QQ' 

( - 2A3 cos((gi - g;)/2)([1122] + [3344] + [1221]24 + [3443]24) ^^'^^^ 
- 2A3 cos((g2 - g2)/2)([1144] + [2233] + [1441]24 + [2332]24 
From this we find 

^4^ + 5'!+^^ = ^ J2 5(g-g' + g-g')([llll] + [2222] + [3333] + [4444]) 

QQ'QQ' 

(6.56) 

which is exactly the four fermion Coulomb term of the Hubbard model in the 
form of ( p:28D . 



6.3 Discussion 



The last two sections were very formal and it is time to get some physical insight 
into the properties of our new formulation. In this section we will address the 
following issues: 

• When first partially bosonizing the Hubbard model, we introduced the orig- 
inal bosons to be able to express physical degrees of freedom by expectation 
values of bosonic fields. Although we simplified our formalism analytically, 
we seem to have lost the intuitive grip on the physical significance of the 
bosons. However, we will show that in the case of homogeneous fields our 
new bosons really describe the physical degrees of freedom we want to in- 
vestigate. 
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• The new formahsm is based on a perturbative one loop calculation. It is 
not yet clear that the full propagator matrix with all quantum fluctuations 
included remains diagonal beyond one loop. We will show that this is 
indeed the case and what we actually have done is to switch to a set of 
bosons that belong to different irreducible representations of translations, 
so that mixing of the bosons is prevented by translational invariance. 

• The one loop expressions offer a first insight into the way quantum fluc- 
tuations are included into bosonic propagators. They give some valuable 
information for the motivation of truncation schemes, e.g.: Which momen- 
tum modes are favored in propagation? How does the bosonic Matsubara 
frequency enter the propagators? How do the propagators behave as func- 
tions of temperature? 

• The last topic of this section will be the periodicity properties of the new 
fields, which will also turn out to play an important role in the definition 
of truncation schemes. 



6.3.1 Homogeneous fields 

The renormalization group analysis we present in the following chapter will be fo- 
cused on the investigation of the properties of the effective potential. The effective 
potential is a function of constant bosonic fields, b{X) — b — const, in position 
space or respectively b{Q) = S{Q)b in momentum space. The interpretation of 
our new bosonic fields in this limit is therefore of great interest to understand 
what the effective potential tells us about the physical degrees of freedom of the 
theory. 

For homogeneous fields, the bosonic kinetic term reads 

= 2n\hlRl + hlRl + hlkl + hlRl) 

+ 2ix\hlf-, + hit, + hl/l + hlX) 
+ 47r2(xldiag(/i2^, hljxi + X2diag(/i'^, 

+ X^diag(/i2^, /i^Jxa + xldiag(/i2^, hljxi). (6.57) 

The Yukawa coupling term becomes 

Sy = - Y,i^\Q){HKA^) + HhlA^) + HhlDo) + R,{hlD^)i:{Q) 
Q 

+ i,\Q)(s,{hlA,a) + UKA.a) + UhlD^a) + UhlD,a))ij{Q) 
_ 1 V (^l^H-Q) (y' r^'^ + (hi sin(g,/2)AA 
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V ^lsm(g2/2)Ci y V"^! cos(g2/2)C2y y"^^^^ 

+ (^/j ■?/;*, sin — > - sin, ■icr2 ^ -i(T2, X^)^ • (6.58) 

We see that in this hmit the bosons R and s couple to the same fermionic bihnears 
as the original bosons p, p etc. This means that for homogeneous fields, we can 
interprete Ri as the charge density, si as the magnetic spin density, S2 as the 
antiferromagnetic spin density etc. 

The interpretation of the complex bosons is not so simple. To understand the 
coupling terms of the complex bosons, we translate them to position space. For 
example, this yields for the first component of Xi 

- l^lxl J2^MX)ta2MX) + MX)ta2MX - ei) 

X 

+ i,,iX)ta2MX) + MX - e^)ta2MX)). 

We represent this result as a pictorial expression in fig. The plaquette 

with position label X is indicated by a dashed line. The results for the other 
complex bosons are given in the same way by the remaining diagrams in fig. 
[H7^ . Particularly we see that Xii ~ X12 describes a d-wave. Similarly, bosons 
with different spatial symmetry properties can be built up by simple inspection 
of these diagrams. 

In conclusion, we have found that despite the transformations we performed 
to diagonalize the propagator matrix, the physical interpretation of the boson is 
as simple as for the original bosons in the most interesting case of homogeneous 
fields. 



6.3.2 Symmetries 



Since our final form of the partition function is equivalent to the original Hubbard 
model for the Hubbard values of the Yukawa couplings, it is clear that it should 
respect the lattice symmetries as well as U{1)- and SU (2)-symmetry of the origi- 
nal model. In principle we could infer the symmetry transformation properties of 
the new bosons from the known transformation behavior of the bilinears (under 
the assumption that the action should be invariant). However, there is no need to 
do this with one important exception: The behavior with respect to translations. 
As we have mentioned in the discussion following ( p.21j ), the translations Tj., Ty 

cLlld Tr^Ty cLCij cLS 



/ 1 \ 

T 

\TxTyJ 



( 



\ 



A, 

Ai exp(igi/2) 
5iexp(ig2/2) 
VBoexp(i(gi + g2)/2)/ 



i,{Q). 



(6.59) 
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Figure 6.3: The complex bosons couple to fermion bilinears that have a position 
space structure as shown. 
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For completeness, we also wrote down the identity transformation 1. Then, to 
guarantee translational invariance, we have 



/ 1 \ 




( ^ ^ 






Ty 


Ki[Q) = 


exp(%/2) 
exp (2^2/2) 


RiiQ) 


\T.Ty) 




^exp(^(gi + q2)/2)J 






( 1 ^ 




( ^ ^ 






T 

Ty 




-exp(2gi/2) 
-exp (2^2/2) 




\T.Ty) 




Vexp(^(gi + g2)/2)y 






1 1 ^ 




/ 1 


\ 




T 


MQ) = 


exp(igi/2) 
-exp (2^2/2) 




MQ) 


\T.Ty) 




V-exp(z(gi + g2)/2)y 




1 1 ^ 




/ 1 


\ 




T 


MQ) = 


-exp(igi/2) 
exp(zg2/2) 




HQ) 






V-exp(z(gi + g2)/2)y 





(6.60) 



The same transformations apply to Sa and Xa- This tells us that the bosons Ra 
belong to inequivalent irreducible one dimensional representations of the transla- 
tion group] The same reasoning holds for Sa and Xa- In other words, to preserve 
translational invariance the full propagator matrix (not only the one loop cor- 
rected one) has to be diagonal in the same sense in which the one loop corrected 
propagator matrix is diagonal (which means diagonal except for the 2 x 2-blocks 
for the complex bosons). In conclusion, f/(l)-invariance tells us that the real and 
complex bosons cannot mix, SU (2)-invariance forbids mixing of bosons from the 
spin singlet and triplet, and finally translational invariance excludes mixing of 
bosons with different a. Note that these invariance arguments do not only hold 
for the propagator matrix, but can also be used to narrow down the possible 
form of arbitrary n-point functions. We will exploit this when writing down a 
truncation of the effective action in the following chapter. 

Translations leave each term of the action separately invariant. This is not 
the case for rotations. Although we know that the action as a whole is invariant 
under rotations, this is not the case for single terms in the action. For example, 
if we only consider the terms in the original action containing the boson qy and 
translate these terms back into the purely fermionic theory (as we did in ( p.32| )), 
we find that these terms give rise to a four fermion interaction term ~ {ip\ii)i + 
'^2'^2 — i^ypi — '^\i^AY- This term breaks rotational invariance. However, the 
conditions for the Yukawa couplings guarantee that these symmetry breaking 
terms cancel each other, but only, if all terms are taken into account. To preserve 
rotational invariance during the flow, we must keep in mind that during the 
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renormahzation group flow the individual terms we wrote down in our action 
wiU not develop wave function renormahzation constants or couplings that are 
independent from each other. For example, the condition that hq^ = hq^ at 
the beginning of the flow remains true during the whole flow due to rotational 
invariance. In the same way, the ratio of wave function renormahzation constants 
for say Ri and Ri, i G {2, 3, 4} will always be constant, if we truncate the wave 
function renormalizations to be momentum independent. However, symmetry 
considerations do not keep mass terms or couplings in the effective potential from 
flowing independently, since the effective potential is a function of homogeneous 
flelds. For homogeneous flelds, the new bosons coincide with the original ones, 
and symmetry transformations do no longer mix different terms. 

The last symmetry we want to discuss is the "time reversal" ( |3.22[ ). In the 
limit of spatially homogeneous bosonic flelds Ran{q = 0) etc. the corresponding 
symmetry transformations for the bosons are simply 

TaMO) = diag(-l, 1, -1, l)R.niO) 
TbMO) = diag(-l, 1, 1, -l)/2_„(0) 

TbMO) = diag(-l, -1, 1, l)i?-„(0). (6.61) 



Particularly, 



TA,TB,TB,Ran{0) = -Rai-n){0) Va. (6.62) 



The last equation also holds if R is replaced by s or x. We can use this sym- 
metry to argue that in any term of the effective action that only depends on 
homogeneous flelds Ra{n=o){Q = 0), etc. the number of bosonic flelds must be 
even. 



6.3.3 Discussion of the momentum dependence 

Before our flnal diagonalization step carried out in sec. |6.1.5| , the bosonic prop- 
agator terms in the action (as functionals of the original bosonic flelds or of the 
flelds bi) were pure mass terms. The momentum dependence emerged at one loop 



level as discussed in sec. |6.1.4| . However, in our flnal form of the theory, already 
the classical propagators ( |6.42| ), ( |6.43| ) are momentum dependent. For high tem- 
perature, this momentum dependence of the classical propagators dominates over 
the momentum dependence induced by quantum corrections. 

In the truncation for the effective action that we will use in this work, we 
will keep the Yukawa couplings constant at their initial values. If we insert 
the initial Hubbard values for the Yukawa couplings, the propagators reduce to 



(|6.47| ). Then the momentum dependence of the propagators for Ri and S2 goes 
as — cos(A;i/2) — cos(/c2/2), which takes its minimum at k = 0. The propagator 
matrix for Xi becomes diagonal and momentum independent. To see how the 
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Figure 6.4: The second derivative of the one loop correction to the bosonic 
kinetic term in the effective action of xii for high temperature T = 1. We have 
set t = 1; the Yukawa couphngs are at their Hubbard model values (|6.46|) with 
= =10. We have plotted the Matsubara mode m = 0. 




Figure 6.5: The second derivative of the one loop correction to the bosonic 
kinetic term in the effective action of X22 ^or high temperature T = 1. We have 
set t = 1; the Yukawa couplings are at their Hubbard model values ( |6.46| ) with 
hi = = 10. We have plotted the Matsubara mode m = 0. 
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momentum dependence emerges for xi, we investigate the one loop correction 
to the diagonal elements of the propagator matrix of xi the case that the 
Yukawa couplings are set to their Hubbard model values (we will neglect the 



off diagonal elements in our truncation). The results are shown in figs. 



JTSj. In both cases, the minimum of the propagator correction again occurs in 



A; = 0. Thus we find that the dominating momentum dependence of the bosonic 
propagators of our final theory exhibits minima at k = for the bosons describing 
the charge density, antiferromagnetism and d-wave superconductivity in the limit 
of homogeneous fields. We will use this fact to argue that the propagation of 
momentum modes close to = is facilitated, so that we can approximate these 
fields homogeneously by their k = 0-mode. 

6.3.4 Periodicities 

The bosonic integrals extend over the range [— 7r,7r]. This range originated in 
the periodicity of the integrand in our definitions of the Fourier transforms. It 
is interesting to analyze the periodicity properties of our new bosons and to see 
how the periodicity of the integrand is maintained. 

The periodicity behavior of the new bosons can be inferred from the known 
behavior of the original bosons, which follows directly from the form of the Fourier 
transforms. By applying all the transformations we performed to arrive at our 
new bosons, we find after a tedious but straightforward calculation 

R{K + 2Txnei) = {-1)''/^R{K) 

x{K + 2nne,) = {-ir/'x{K) (6.63) 

for n even, i G {1,2} and 

R{K + 27mei) = i{-l)^''+^^/'^BiR{K), 
R{K + 27me2) = BoR{K), 

b 

XaAK + 2ixne,) = z(-l)("+i)/2 5^(5o)aa..(i^) 

b 

XaM + 2^ne,)=z{-lY-^'^/'J2^B,),,XbM) 

b 

XaM + 27^62) = J2(Bo)abXbM) (6-64) 
b 

for n odd. The equations for R also hold respectively for s. Apart from a possible 
phase factor, the fields are periodic in Air. They are not periodic in uneven 
multiples of 27r, but transform into each other. Thus it is clear that no single 
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term in the bosonic part of the action is periodic in 2ti. However, periodicity of 
the integrand is restored when considering sums of terms. For example, whereas 

]^Rr{-K)P^MK) (6.65) 

is not periodic in 27r, the sum 

\Y.R^i-K)PaMK) (6.66) 

ab 

is periodic in 2tt. 

The interesting thing about this behavior is that it is possible to write terms 
as higher momentum modes of other terms. For example, we can write 

(6.67) 

for any c G {1,2,3,4}. We have shown that the same property holds for the 
propagator terms of the complex bosons, as well as for the coupling terms. One 
possible equivalent transcription of the bosonic terms in (|6.40 ) that we will use 
is 

K 

Sy = -J2J2 [s{K-Q + Q')(^MK)^i{Q)V,l,{K)MQ') 

KQQ' ab 

+ UK)i^l{Q)VX^{K)UQ'] 
+ 6{K-Q- Q') [x\{K)i,l{Q)vX,{Q, Q')MQ') 

+ xl{K)^i{Q)VX,{Q,Q')rb{Q'] 



(6.68) 



where it is understood that bosonic momentum integrals extend over the interval 
[— 27r, 27r] (that is, = T J2n /-L (1^' ^ denotes the momentum of a boson) 
and fermionic integrals over [— vr, tt]. The 5-function is assumed to be periodic in 
27r for Q, Q' and in Air for K. Only the physically interesting bosons Ri, S2 and 
Xi enter our formulation in this transcription. The other bosons are included as 
higher momentum modes of these three bosons. 



94 Chapter 6. Diagonalization of the propagator matrix 



Chapter 7 

Renormalization group analysis 



The history of renormahzation group approaches to the Hubbard model is a short 
one [g 



Up to now all these investigations were performed in the framework of 
the purely fermionic model. In this section we will see how to apply the renormal- 
ization group formalism developed in chapter ^ to our final partially bosonized 
form of the partition function (|6.38| ), where the action terms are given by (|6.4CI|) . 
If convenient, we may rewrite the action terms by using ( 6.6^ ). The first task is 
then to write down a suitable truncation for the effective action. This truncation 
ansatz is then to be inserted in ( [5.24| ), from which we can derive the flow equa- 
tions for masses, couplings and wave function renormalization constants. Our 
renormalization group analysis will be focused on the properties of the bosonic 
effective potential, and the truncations proposed in the first section will be ad- 
justed to this aim. 



7.1 The truncation 



Since ( p.24|) cannot be solved exactly, we propose an ansatz for its solution, which 



is a truncated version of the effective average action. The guidelines for doing so 
are the following: 

• Due to the property lim^^A = S the effective average action will resemble 
the classical action as a function of expectation values of the fields at the 
beginning of the flow. The truncation will therefore include terms that 
look like the corresponding terms of the classical action and systematic 
generalizations of these. 

• The generalizations are limited by the fact that they should respect the 
symmetries of the theory. 

• We only keep terms that seem to be absolutely necessary to describe the 
behavior of quantities we want to calculate (in our case bosonic masses and 
quartic couplings). 
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It is clear that by truncating the effective average action we introduce errors 
due to the approximation. The most difficult part in a successful renormalization 
group analysis is to decide which terms to include in the truncation. This requires 
physical intuition, some systematic expansion of the terms which possibly appear 
in the effective average action and a trial and error procedure — at last, the 
only reliable way to estimate truncation errors is to include more terms in the 
systematic expansion and to see how they change the results. We want to stress 
that for theories with strong couplings (so that no perturbative expansion with 
respect to some small quantity is possible) there is no way to circumvent these 
approximation problems. One large advantage of the method of the effective 
average action we present here is that the property lim^^A Tfc = 5* enhances our 
intuitive grip on possible truncation schemes. 

In general, the effective average action can be written in the form 



where contains only bosonic fields, only fermionic fields and T^^ coupling 
terms between bosonic and fermionic fields. Since we are mostly interested in the 
properties of the effective potential, which is part of Ff , we propose the following 
simple ansatz for Ff and Ff-^: 



The reasoning behind this approximation is that we are mainly interested in the 
flow of the effective potential. Terms with more than two fermionic fields do not 
contribute to the flow of the effective potential. However, they do contribute to 
the flow of the Yukawa couplings — that we keep constant in our truncation (see 
below). For constant couplings we can therefore ignore all terms with more than 
two fermionic fields. This leaves one fermionic propagator term which is part of 
Fj^ and coupling terms of two fermionic fields to an arbitrary number of bosonic 
fields in T^^ . The propagator term can be written as the classical propagator 
term times a momentum dependent wave function renormalization constant. To 
arrive at (|7.2| ), we additionally make the following approximations: 

• The fermionic wave function renormalization constant Z[ is kept constant 
at its initial value Z[ = 1. 

• All terms with two fermions and more than one boson are neglected. 

• The Yukawa couplings are kept constant. This approximation would not be 
too good for most simpler theories, where the initial values of the couplings 
are known. It is even worse in our case, since the many different choices of 
the initial values of the couplings (although equivalent if everything is exact) 
lead to different results if approximations are made, as we already saw in 
the mean field case. However, to get a first impression of the properties of 



F, = Ff + Ff + F: 



k ; 




pF _ nF -pBF _ Q 

^ k — '-'kin^ k — '^y- 



(7.2) 
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the effective potential, we will nevertheless take the Yukawa couplings as 
parameters and discuss the results as functions of these parameters. The 
inclusion of the flow of the couplings into the renormalization group analysis 
is subject to current work ||r?| . 

Again recall that with all these approximations the effective average action terms 
involving fermionic fields necessarily coincide with the corresponding terms of the 
classical action due to the property lim^^A Tfc = S. 

We now turn our attention to Ff . Ff is a functional of the bosonic fields 
bi G {Ra,Sa,Xa,xl} with a G {1,2,3,4}. A systematic general expansion of Ff 
in powers of the fields reads 

K1K2 1112 

+ ■■■ (7.3) 

with /c-dependent coefficients {k is the flow parameter from chapter Each 
coefficient contains a momentum conserving ^-function. From ( |6.63|) we know 
that 

h{K) = -h{k + Aney2) (7.4) 

which implies 

Cn...^„,k{Ku ... ,Kj+ 47rei/2, . . . , K^) = (-l)"Q,...i„,fc(iri, . . . , ir„) (7.5) 

to preserve the periodicity of the integrands. This tells us that we can expand each 
coefficient Cj^...j„^fc(-^i5 • • • ,Kn) with respect to the functions {Kj = {{u!^)j, kj), 
kj = ((%)i,(%)2)) 

{1, cos(mj(A;_,)i/2) — 1, cos(mj(A;j)2/2) — 1, sin(mj(A;j)i/2), sin(mj(A;_,)2/2)} 

(7.6) 

for n even and 

{cos((2m, + l)(A;j)i/4),cos((2mj + l)(%)2/4), sin((2mj + l)(A;,)i/4), sin((2m, + l){kj)2/4:)} 

(7.7) 

for n odd and with m G N (we have chosen to expand with respect to cos(mA;i/2/2) — 
1 instead of cos(mA;i/2/2) to achieve that all momentum dependent terms for even 
n vanish for k = 0). The coefficients Cm,k of this new expansion with respect to 
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trigonometric functions only depend on the Matsubara frequency and the flow 
parameter k. 

In our truncation of we include a momentum dependent propagator term 
(more on this later) and a second term containing all other bosonic terms 



rf = rL,. + u,- (7.8) 

We simplify the coefficients Cm,k by neglecting their Matsubara mode dependence. 
Furthermore, for homogeneous fields we have Ki = OVi. Then (|6.62| ) tells us that 



all terms of the effective action with an odd number of fields have to vanish 
to preserve "time reversal" -invarianceQ. For the remaining terms with an even 
number of fields, only the terms with momentum independent coefficients survive, 
so that we truncate Uk in the form 

UkM = UoV+ $^Q,,,,,6,,(iri)6,,(i^2)^(...) 

K1K2 1112 
K1K2KSK4 jlj2«3«4 

+ ••• (7.9) 

with homogeneous fields bi{K) = bi5{K), coefficients that only depend on the 
fiow parameter k and appropriate momentum conserving 5-functions 6{. . .). 

Due to U{1)- and 5'f/(2)-symmetry the effective potential Uk can only depend 
on the invariants 

Pab{KuK2) = ^Ra{Ki)R,{K2) 
aab{Ki,K2) = ^Sa{Ki)Sb{K2) 

Ki,bAKi,K2) = xUKi)XbAK2). (7.10) 

By use of ( |6.68| ) we can write the effective potential as a function of pn, 022 and 
Aii_ii, Aii_i2, Ai2,ii, Ai2,i2 only. To simplify the notation, we write Aab := Ai^^i;,. 
Recall that all other bosons are included as higher momentum modes of the 
bosons explicitly present in our truncation as discussed after ( |6.68| ). This means 
that although we will only consider homogeneous modes of the invariants pn, 
022 and All, ^12, ^21, ^22 as external lines, the exchange of all virtual bosons 



^Note the tricky part of this argument: The symmetry transformation ( |6.62D hved in the 
space of bosonic fields, which are all independent. If we switch to the effective action, the 
(linear) symmetry transformation carries over to the expectation values of the fields, except 
for the charge density i?i(0) oc p. The charge density expectation value p is regarded as a 
parameter controlled by the source p that is no longer explicitly present in the effective action. 
The transformation p — s- — p does not follo w from the transformation behavior of p, since it is 



no longer a free field, but from p — > —p in ( 3.22 ) and ( 4.24 ) 
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is included as we also integrate over the higher virtual boson momentum modes 
representing the bosons not explicitly present in the boson set we selected. 

The momentum dependent term in ^ is truncated to a propagator term 
only: 

KK' 

+ X\{K)P^,^,{K,K')UK')\ (7.11) 
For the real bosons we set 

Plk{K, K') = {27:fZ^{P,%K) - P^{K = 0) + (mvrT)2)5(K + K') 
Pi,^,{K, K') = {27i fzi{Pl{K) - Pl{K = 0) + {m'nTf )5{K + K'). (7.12) 

This has to be compared to the classical action propagators 

P^^{K, K') = {2T,fP^{K)5{K + K') 

P^K, K') = {2TxfPi{K)5{K + K'). (7.13) 

First note that the terms P^{K = 0) and Pi{K = 0) are momentum independent 
and therefore can be compensated by adjusting the coefficients in the effective 
potential Uk, so that we are free to add them. If we note that the propagators 
Pf{K) and P^iK) (for the Hubbard model values of the Yukawa couplings) 
become minimal in i^' = 0, we see that (for m = 0) the advantage of this addition 
is to make the propagator vanish at zero momentum and positive otherwise. This 
allows to define simple truncation schemes for the bosonic propagators. 

As we will see below, we use a temperature like cutoff in the fermionic sector. 
This means that during the beginning of the flow the system behaves as in the high 
temperature limit. The interesting physics emerges gradually as the temperature 
cutoff is lowered. 

The last term we added involving the Matsubara frequency m is needed to 
make loop Matsubara sums finite (since we will keep the Yukawa couplings con- 
stant, the complete m-dependence of loops will be provided by the propagators). 
By adding a term ~ m^, we mimic the low temperature behavior of the one 



loop result as discussed in section |6.1.4|. We fit our truncation to the original 



action term by comparing the m = 0-Matsubara mode (in the high temperature 
limit, this is the only one contributing in the one loop calculation). We therefore 
set = 1 at the beginning of the flow, where k is large. Note that we have 
truncated Z^, Z| to be momentum independent. 

In principle there is nothing that prevents us from replacing (mvrT)^ by 
Cfc(m'7rT)^, where Ck is a /c-dependent quantity. We have checked numerically 
that such a factor only has a small effect on our results, so that we set it equal 
to unity. 



100 



Chapter 7. Renormalization group analysis 



For the complex bosons, there is no classical propagator that can be used as 
the main ingredient of the truncation, since for the Hubbard model values of the 
couplings the classical propagator of the complex bosons is momentum indepen- 
dent. We therefore propose a simple ansatz taking into account the first Fourier 
term in a general expansion of the effective action with respect to momentum 



iA3(2 - cos(A;i/2) - cos(A;2/2)) + {mixTf ) 5{K - K'). 



(7.14) 

Since we know that the first momentum dependent contributions to the propa- 
gators of the complex bosons emerge at one loop level in 0{1/T^) (cf. ( |6.29| )), 
and our regularization scheme in the fermionic sector (as discussed below) will 
replace mvrT by mirTk = m7r(T + /c^), the momentum dependence develops as 
0{l/k^). We therefore approximately set Z^^j^ = 1/A^ as the initial value of the 
wave function renormalization constant. 

In conclusion, our truncation for the effective average action reads 

-^BF 



kin.k 



Y,i'\Q)P'^{Q.Q')m) 



QQ' 



KK' 



Uk[bi 



Uk[pn, 022, All, Ai2, A21, A22] 

KQQ' ab 

+ 6{K-Q- Q') {x\{KWa{Q)vX,{Q, Q')MQ') 

+ xI{K)iji{Q)VX,{Q,Q')rb{Q') 



(7.15) 



with 

pHQ,Q' 

,R 

11, fc 



2t(cos(gi/2)Ai + cos{q2/2) Bi))6{Q - Q') 



Pf\,{K, K') 



Pkk^K.K') 



(2TifZ^{P^{K) - PP{K = 0) + {mnTy)6{K + K') 
{27ryzi{Pi{K) - PiiK = 0) + {m7TTf)6{K + K') 



P^M.K') = {2TifZl -A3(2 - cos(fci/2) - cos(A;2/2)) + (rnvrT) 



5{K - K') 
(7.16) 



and the vertices given by ( |6.44| ) and (|6.45|) . 
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7.2 The flow of the eff^ective potential 
7.2.1 The second derivative of Tk 

We will now start to derive the flow equations for the effective potential. In 
general, the flow equations for any interesting quantity can be derived from our 
master equation ( |5.24| ) for some particular truncation (in our case ( |7.15| )). As we 
see, we need the second derivative of with respect to the fields. Explicitly, we 
need 

/ 



Tf\K,K') 













SuHK)5u{K') 


5u-\{K)5u''{K') 


5u^K)5w{K') 


SuHK)5-ij){K') 


Su'f (K)5ip'- {K') 
&^Tk 


5uT(K)Su(K') 


SuT(K)5u* (K') 


5uT{K)5w{K') 


5uT(K)S^(K') 


SuT(K)5iP*{K') 
S^Vk 


5wT{K)5u{K') 


5wT(K)5u* (K') 


Sw'^{K)Sw{K') 


SwT (K)5il)(K') 


5wT{K)5iP*{K') 

5 


5ipt (K)Su{K') 


5ip'r{K]Su* (K') 


Sil)"! {K)5w{K') 


S^I)^K)5il)(K') 


S^Tk 


S-il)T (K)5u{K') 


5tl>T(K)Su*{K') 


{K)5w{K') 


5ipT{K)S^{K') 





(7.17) 

{'2) 

where T\. {K, K') denotes the matrix of second derivatives occurring in ( |5.24| ) 
(we use this notation to avoid confusion between the complex bosonic fields x we 
use here and the generalized fields x i^i ( |5.24[ )). 
In our truncation, the kinetic terms yield 



/O \ 





P "^ 

\0 P"^) 



{K, K') 



kin,kJ 



l(2) 



{K, K') 



11, fc 







V 





pX 

11,' 













diag(Pj^.,P; 



11, fc' 22,k) 






0\ 



0/ 



{K,K'), (7.18) 



where we have used that the propagator matrices are symmetric. 
The contributions of the coupling terms are 



64{K)6xm 

5x^{K)5ram 

5^l{K)5x*{K') 



Q b 
Q b 

J2J2S{K-K' + Q)VX,{K, Q)MQ) 

b 
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2-rBF 



6'T 



with 



x2-nBF 

-^-^^^^ = - E E ^(^ + ^' - 

- E E '^(^ - ^' - Q)vi'{Q>m 

Q 7 

- ^ + X' - Q)xkQ)V;^,,i(i^, K') 
Q 

J2S{K + K'-Q)x^{Q)VlAK,K') 
Q 

- E E '^(^ - ^' + Q)n?(Q)^7 (7.19) 



V{K, K') = X') - V^{K', K) (7.20) 

and ty-y e S2}. 

The effective potential Uk depends only on the invariants pn, q;22, Ajj. It 
is therefore convenient to express the derivatives with respect to the fields by 
derivatives with respect to the invariants. Let 



2 \5pii{K,K') Spu{K',K) 



2 \5a22{K,K') 5a22{K',K) 
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Tlien we can write (generalized matrix notation!) 

f/f = diag(DA, Dl, Dp, D^, 0, 0)Uk (7.22) 



+ 



/Da(xX^)^a 


da{xx^: 


Dl 


Da{xRi)Dp 


DA{x^2)Da 





o\ 


Dl{x*X^)D^ 


{Da{xX^I 


DaV 


Dl{x*Ri)Dp 


Da{x*^)D^ 








Dp{R,x^)DA 


Dp{RiX^ 


)Dl 


Dp{R,Ri)Dp 


Dp{Ri^)D^ 








Dc.{s2X^)Da 


Dais2X^ 


)Dl 


D^{s2Ri)Dp 


Da{s2^)Da 




























V 
















0/ 



7.2.2 The flow of the effective potential 



In this section we will derive the flow of the effective potential from (|5.24 ) or 



equivalently, the flow of the coefficients Ci^,,,i„ in ( |7.9| ). In general, there is an 
infinite number of coefficients. To be able to calculate anything, there are two 
possible ways to render the number of flow equations finite: 

• Truncate the effective potential at some finite power of the fields, e.g. at 
quartic order in the fields. This yields a finite set of flow equations for the 
remaining coefficients. 

• Instead of considering the flow of the coefficients, consider the flow of the 
full potential as a function of discretized homogeneous fields. For example, 
if U[(j)] is an effective potential depending on the (continuous valued) field 
0, one could analyze the flow of the potential at given points ?7[0 = 0o]) 

= = 02] etc. A finite number of points yields a finite number 

of flow equations. 

Both methods have their advantages and drawbacks. The first method is easier to 
apply, more stable and faster in the numerical treatment, but has the drawback 
that — since it is nothing else than a polynomial series expansion around some 
given point to a given order — the possible solutions for the effective potential are 
reliable only close to the point around which we expand. Usually, one expands 
around the minimum of the effective potential. Physical properties that can be 
inferred by looking at the flow of the potential close to this minimum can be well 
understood by this approach. However, phase transitions of first order typically 
can not be derived from properties of the flow of the potential near the minimum 
(cf. fig. O)- 

The second method does not have this kind of problem and any possible 
shape of the potential can be described by it — not only those corresponding to 
some finite power expansion in the fields. But on the other hand the numerical 
implementation of this second method is far from trivial. The reason for this is 
that we need the second derivatives of the effective potential on the right hand 
side of the flow equation. If we discretize the region the effective potential is 
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defined on, these derivatives also have to be discretized. As is already known 
from more simple systems, it is a formidable task to achieve this discretization 
in a stable way. We will therefore use the first method, keeping in mind that we 
possibly face difficulties with first order phase transitions and postpone the more 
general treatment by the second method to future work. The implementation 
of the first method requires an explicit truncation of the effective potential. We 
will investigate different truncation schemes for the effective potential later on. 
In this section, we will derive the flow equations as far as possible without fixing 
some specific truncation. 

The idea to extract the flow of the effective action from ( 5.24|) is to note that 

(7.23) 



lim lim = Uk 

ip,'<p'^Ob,{Q)~*h,S{Q) 



SO that 



dk 



- lim lim Tr (dkRk)M 

2 ip,ip*^0h{Q)-^h6{Q) * 



From the results of the last section, we have 



lim F 



(2) 



B 






F 



(7.24) 



(7.25) 



with 



B 



diag(Pi\,, + D^Uk, iP^J'^ + DlUk. Pi, -r -p.., . 22,. 

f D^{xx'')DAUk D^{xx^)DlUk D^{xRi)DpUk D^{x^)D^Uk\ 

DI{x*X')DaUu {DA(xx')DAVUk Dl(x'Ri)DpUk DA{x*^)D,Mk 

Dp{Rix'<)DAUk Dp{Rix^)DlUk Dp(RiRi)DpUk Dp{Ris^)DcUk 

\D^{s2X^)DAUt, Do,{s2X^)DlUk Dc(s2Ri)DpUk D„(s2sf )Dcl/fc/ 



+ DaUk 



+ 



F{K, K') = diag(P'^(ir, K'), P^{K, K')) 



+ 



-EqHQ-k-k') 



Eq 6iQ + K- K') 



(7.26) 



(2) 

We see that lim^^^.^orfc is block diagonal in the bosonic and fermionic sector. 
With 

Pf = diag(i?^, {Rlf, ^r), Rk = diag(i?r, (R^f) (7-27) 
the flow equation now reads 

Tr ((Mf) + - Tr (^(Mf) (f + P^)"' 

(7.28 



—Uk = - lim 

dk 2 bi(Q)^bi5{Q) 
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We separate the flow equation by Uk = Uf: + , so that 

4rU^ = I lim Tr ( (dkR^) (b + R^Y'] 

4rU[ = lim Tr ((dkR^) (f + R^Y'] . (7.29) 

The right hand side of the flow equation for U[ does not depend on Uk, so that 
we can further calculate it without specifying a truncation for Uk- 
It is very instructive to rewrite the flow equation for as 

^U^ = hm 9,lndet(F + /?f), (7.30) 

dk 2 bi(Q)^biS{Q) 

where 

If the right hand side depended on k only via R^ , we could replace dk by ^ and 
immediately integrate the equation to get Uk as a function of k. For A; — 0, we 
know that i?f vanishes. In this case we have 

UCo = U^-^j, - I hm In det F (7.32) 

which is nothing else than the one loop corrected potential we calculated in our 
mean field analysis (cf. ( [4.10| )) — besides the fact that we used a different set of 
bosons there. 

The calculation mainly goes through as in the mean field case. We set 

R^{K) = i{2n + l)'Kk'^ (7.33) 

and finally arrive at 

In cosh (|^\/ {\I'^KPii + ^2^2h^22 + 4t2(ci + eic"^)' + A(ei) j 

(7.34) 

with Cj = cos(gj/2), Tk = T + k"^ and 

A(ei) = /^^^An cos^ | + /^t A22 cos^ | + </^L(Ai2 + A^^e^ cos | cos |. 

(7.35) 
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For Aij = 0, this reproduces the mean field result ( ^.181 ) with 5 = 0, if we set 
Ri = p/hp, S2 = a/ ha and assume that the right hand side depends on k only via 
Tfc. This was to be expected from comparing the original coupling term in ( |3.47[ ) 
with (|6.58|) in the limit of homogeneous fields. The term involving the complex 



bosons looks somewhat different, which is due to our redefinition of the complex 
fields in ( p.l4| ). The original field e corresponds to Xi + X2; whereas the rf-wave 
corresponds to Xi ~ X2- The corresponding invariants are 

e:{x*i + xDiXi + X2) = An + A12 + A21 + A22 

d:ix*i- xDiXi - X2) = An - A12 - A21 + A22. (7.36) 

If we assume that the original boson e always has a vanishing expectation value, 
we can set A12 + A21 = —An — A22. If we further assume that rotational 
invariance is not broken in the sense that An 7^ A22, we can set 

A = 2An = 2A22. (7.37) 



Since we keep all couplings at their Hubbard values, we additionally have h 



hi = 2h'^, so that 



A(ei) =2/i^A(ci-eiC2). (7.38) 



With these assumptions, we therefore are in complete agreement with the mean 
field result, if we additionally set 6 = 2/i^A. 

7.3 The flow of the bosonic wave function renor- 
mahzation constant 

In the following section we will also need the fiow of the wave function renor- 
malization constant Z^. We will derive the corresponding fiow equation in this 
section and show that it can be analytically solved in our truncation. 

First note that can be extracted from our truncation of the effective action 

by 

4 d 52 

Zt = —, lim lim ^, lim lim — ; — - — — — — - lim Fj. = 0'3'™™Fi, (7.39) 

in the symmetric phase and 

Zf = — i— lim hm lim lim lim F^ = O^J^^Vk (7.40) 

in the broken phase, where denotes the minimum of the effective potential, K = 
{uj^, k) and by e we denote some small momentum |e| -C 1. Note that the method 
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of extracting Zf. is not uniquely determined by our truncation. Particularly, the 
limits s — >■ and s* — *■ are introduced only to yield more consistent equations 
when deriving the flow of from ( [5.241 ). By substituting the momentum k by e, 
expansions of the trigonometric functions with respect to e become possible. This 
is very convenient, since the quadratic expansion of limfc^e(ci + C2) is rotationally 
invariant and a function of |e| only. 
From ( |5.24| ) we immediately have 



d,Zl = 0^'"/^=^°'^ ( iTr((9fci?.)M(rf + R,)-') ) . (7.41) 



7.3.1 The symmetric phase 

As we already discussed, the flow equation for the effective average action has 
the form of a one loop equation. The one loop S2-propagator term corrections 
have the diagrammatical form 




for fluctuations in the fermionic sector and 



1; 

K \ / K 



for fluctuations in the bosonic sector. Acting with O^™™ on these expressions, 
the external legs are amputated. Since we have to differentiate with respect to |e| 
(which is nothing else than a small external momentum) and the loop momentum 
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does not depend on K in the second diagram, we conclude that in the symmetric 
phase we do not have any contribution to the flow of Z| from the bosonic sector. 
The flow equation therefore reads 



(7.42) 



where we have already performed the limit ip, ip* ^ 0. Inserting F on the right 
hand side, performing the trace and all limits and derivatives, a lengthy but 
straightforward calculation yields 



with 



9 , tanh 



te. = ^^^^^7^^, 3. = 2tsin(g,/2). (7.44) 



Since in our truncation the right hand side of the flow equation depends on k via 
Tfc only, we have dk = dk, so that we can directly integrate the flow equation to 
get 



vr^Aa (27r) 



2 X] rp2^e^,k{Q) 

£1 



(7.45) 



From (|7.43|) and ( [7.45|) , we immediately find the anomalous dimension 



?7fc = -kdklTiZl 



kdkZj 



(7.46) 
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7.3.2 The broken phase 

In the broken phase, an additional diagram contributes to the one loop correction 
of the bosonic propagator term. 

4X Q 



\ 



/ \ 

\ 



/ \ 



^ / K-Q 



\ 

\ 



Here we have the possibility of a coupling to a non vanishing external field s^, so 
that we have a contribution from the bosonic sector in contrast to the symmetric 
case. The broken phase appears for small values of k. Since the flow in the 
fermionic sector goes as kd^Z^ oc fc^, it is a good approximation to neglect the 
contribution from the fermionic sector in the broken phase. Furthermore, as we 
will see when discussing the flow equations for the effective potential below, only 
the Matsubara mode m = and momenta q close to zero contribute to the flow 
at small k. This means that — with regard to the flow of — we effectively 
face a two dimensional bosonic theory, which can be mapped to a simple 0(3)- 
model in two dimensions. But for the 0(3)-model the flow of the wave function 
renormalization constant has already been calculated to be 



Vk = ^«°2,.(AD^mi,(2A^a°2,„ 0) (7.47) 

where V2 = l/(87r) and Tnl 2{'^^k^22,ky^) a so called threshold function whose 
exact form depends on the choice of the cutoff function (again, for more details 
we refer to pO[). If we choose Rf.{q) = Z^k"^ — q^)6{k'^ — q"^), the threshold 



function simply reads 



m|2(2A^<., 0) = ^2 • (7-48) 

With some minor adjustments reflecting different constant factors in our formu- 



lation as compared to the treatment of the 0(3)-model in |2^, the anomalous 
dimension finally reads 

AT 1 
Zl can then be derived from dkZ^ = —Zf,rik by numerical integration. 
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7.4 Results for different truncations of the ef- 
fective potential 

A lot of interesting questions can be addressed by the formalism we developed so 
far. We will restrict our attention to two aspects: 



The Mermin- Wagner theorem tells us that in two dimensions at non 
vanishing temperature no long range order is present. However, in the mean 
field case we observed phase transitions for non vanishing temperature. 
Furthermore, in the effectively two dimensional cuprates antiferromagnetic 
and superconducting behavior is actually being observed. In a very simple 
setting (half filling, only antiferromagnetic degrees of freedom taken into 
account) we will see that this puzzle is solved by the renormalization group 
analysis. 



• Mean field calculations tend to overemphasize symmetry breaking. It is 
interesting to study the influence of bosonic fluctuations on the phase di- 
agram to better understand the basic mechanisms of the interplay of the 
different degrees of freedom. 



7.4.1 Antiferromagnetic behavior at half fiUing 

At half filling, we have p = pu = 0. Furthermore, we set A = 0, since we are 
only interested in the antiferromagnetic behavior. 

Fermionic contribution to the flow 

In this truncation, the flow equation for reads 



with 



k 



SkTV ^ ^^(«22) tanh(/,,(a22)) (7.50) 



/.,(a22) = ^\/2/^a«22 + 4t2(ci + e.c.f . (7.51) 



We have replaced ha by ha to indicate that this coupling is not renormalized. 



7.4. Results for different truncations of the effective potential 



111 



Truncation of the potential 

We work with a simple quartic truncation of the potential. We set 



K1K2 



\ Ka22{KuK2)a22{Ks,KMKi + K2 + K3 + K^) (7.52) 



2 

K1K2K3K4 



in the symmetric phase and 



2 

(a22(i^i,i^2) -a^2,fc<^(^i)<^(^2))(«22(i^3,i^4) -a^2,fc<^(^3)5(ir4)) (7.53) 
in the broken phase. For homogeneous fields, these equations become 

Uk = V (Uo + «)'a22 + \\>l2) (7.54) 



and 



U, = ^V\Ua22 - al2,k? (7.55) 



plays the role of a fc-dependent mass, k is the minimum of the potential 
if the symmetry is spontaneously broken and is the quartic bosonic coupling. 

Bosonic contribution to the flow 

For our truncation, the matrix B in ( [r.26| ) becomes 

B + R^ = dmgiP^,^, + Rl (P^, + Rlf, P,l, + Rl P^, + Rl + (D„ + D^is2^)D, 

(7.56) 

If we neglect all terms in 

-jrUk = Idk lim Trln (b + R^] (7.57) 
that do not depend on fields, we get 

71^^ = V^'^ . J'^^.n. ^'^"^ (^^Ifc + Ri + {D^ + D^{s2sl)Do.)Uk) . (7.58) 
If we now insert our truncation for Uk and perform the trace, we arrive at 

^^u^ = E ^ E r HPkkiQ) + RiiQ) + 4), (7.59) 

A ryy~, '2,11 V / 
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where Q = {u^, q) and 



{{fhlY + 3A^a22, {iTT^kT + ^k^22, {'rrikY + for the symmetric phase 

A^(3a;22 — «22 fc; '^22 — «22 A;; "^22 " «22 fc); ^^^^ broken phase. 

(7.60) 



Choice of the cutoff function 

We define the cutoff function Rl{Q) by 

RiiQ) = i2n)\'zi (k' - iA3(2 - c, - c^) - m^T' 

e (k^ - iA3(2 - ci - C2) - tu'tA (7.61) 



where 0{x) is the usual step function^. This cutoff function respects (|5.15|) and is 
therefore a viable choice. The step function cuts off large Matsubara frequencies 
and momenta. By lowering fc, we average over larger and larger regions in position 
space. For this cutoff, we therefore can relate properties of the effective action at 
some given value of k with long range order properties of the theory in position 
space at scales of order 1/k. This will become important when interpreting the 
results of the renormalization group analysis with respect to the Mermin- Wagner 
theorem puzzle. Note however that from a purely numerical standpoint our 
choice of the cutoff function is not optimal, since the cutoff for the Matsubara 
sum is momentum dependent, so that the Matsubara sum cannot be carried out 
analytically. In the next section, where we consider a different truncation and 
are not dependent on the position space interpretation of the flow, we will use a 
different cutoff for which the Matsubara sum can be evaluated analytically. 
Inserting the cutoff function into ( [7.59|) , we flnally arrive at 



dk^ k ^ ^ J (27r)2 

i m=-M ^ ' 

-Vk{k' - iA3(2 - cos(gi/2) - cos(g2/2)) - (mT)^) + 2^ 
-i , (7-63) 

1.2 _|_ ^fc 

^ (27r)2,r2Z^ 

where M = max{m G N|m < k/T} and the momentum integration runs over 
{{qi,q2}h e [0,2n],e ~ {X^/2){2 - ci - c^) - m^T^ > 0}. % = -kd^lnZiis 
the anomalous dimension. 



^Explicitly, it is given by 
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Extraction of the coefficients 



The flow equations for the coefficients follow from the flow equations of the ef- 
fective potential by 



1 , d 

— hm 

y 022^0 aa22 

1 , <P 

— lim 

V 022^0 da22 



dk 



d_ 

dk 



(7.64) 



for the symmetric phase and 



di.aZ 



— lim 



lim 

d^ 



d 



at!. da22 



dk 



a22^a22,fe dal2 



d_ 

dk 



(7.65) 



for the broken phase. Note that if we set 0:22 = in the symmetric phase 
expression of Cfc, we see that Ck — > {ifnlY, {'^k)'^i if'^kY) describes three modes 
with equal mass. In contrast to this, if we set 022 = «22 k broken phase 

expression of Cfc, we find Ck = (2X^0:22 fc; 0; 0); ^^^^ massive mode and 

two massless modes. Of course, the massless modes are nothing else than the 
Goldstone modes which have to appear when breaking SU{2) down to U{1). 



Introduction of rescaled, renormalized quantities 

It is very convenient to introduce rescaled, renormalized quantities defined by 



[m 



a\2 

k) 



a\2 

k) 



a 



22,k 



^fc"22,fc' 



"'a 



(7.66) 



The advantage of rewriting the flow equations in terms of these variables is that 
they (at least for small k) no longer depend on and k explicitly. 



The fiovi^ equations 



Using ( [7.64 ) and ( 7.65| ), we can derive the flow equations for the rescaled and 
renormalized quantities ( [7.66|) (except for the coupling ha that as we said will be 
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kept constant) from ( [7.50[ ) and ( |7.63| ). They read 

2Thl ^ r _fq_ / tanh^(O) 1 

T3 Z^j_^(2^f\ /^^(O) +cosh^ 



lOTA^ 

(27r)%2 



E 

m=—M 



dPq 1 2 - r^fc (1 - ^(2 - Ci - C2) - (mT/A;)^) 
(27r)2 A;2 



1 + 



(m-)2 



(27r)27r2 

r dS /^(0)-sinh^(0)cosh^(0) 



2T5 



(27r) 



^(0)3coshV..(0) 



2tanh^(0) 
^(0) cosh (0) 



+ 



44T(Ag)^ 
(27r)%4 



M 



m=— M 

-2A^(l-r/fc) 
in the symmetric phase and 



d\ 1 2-r;fc(l-^(2-ci-C2)-(mT/fc)2) 
(27r)2 F 



1 ^ ("^g)'. 



(7.67) 



2Tht ^ _fq_ I tanh/,,(a°,,) ^ 



cosh^ (aS 



+ 



2T 



(27r)%2 Z 



2 + 



^(27r)2|^ /,,«,) 



(27r)27r2 



„0 



2T| 



IE 



.(27r) 



/ej fc)~'''"h/e-^ (q^j fe) ^^1 ('^M.fc' 2 tanh /ej ("22, fe' 

/^l("22,fc)^'=°*2'/,^{<.^2,fe) 



n^"22,fc^ 



+ 4T 



(27r)%^ 



2 Af 
m=—M 



2 + 



1 + 



^^fc"22,fc 
(27r)27r2 



in the broken phase with dt = kdk (that is, t = ln(A;/A)), 
max{m G N|m < k/T}, q = cos(gj/2), r^jt = —kdklia Zl, 

/.i(a22) = l/(2Tfe) v/2/i4a22 + ^t\ci + eiCa)^, 



(7.68) 
T + k\ M = 

(7.69) 
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and the integration runs over 

{{qi, q2}h e [O, 27r], e - (A3/2)(2 - ci - C2) - m'T^ > 0}. (7.70) 

The initial conditions are (wi^)^ ~ ((27r)^/ia)/A^, where A is the initial value of 
/c, = and = 1. We did not write down a flow equation for Z^, since Z^ 
and the anomalous dimension r]k can be calculated analytically in our truncation 



as we saw in sec 7.3. 



Remarks 

For small k, we have Tj. tuT and the integral in the bosonic sector 



d'q _ )^2 - r^, (^1 - ^(2 - ci - c^) - {mT/kYj j (7.71) 

can be evaluated analytically. First note that for small k we have M = 0, so that 
only the m = 0-mode contributes to this integral. The integration region is then 
given by 

{{qi,q2}\qi e [0,27r], fc^ - (A3/2)(2 - ci - c^) > 0}. (7.72) 

Thus for small k, only momenta close to g = contribute, and we can expand 
the trigonometric functions to quadratic order around q = 0. It is then simple 
to perform the integration and to see that all explicit fc-dependences cancel. For 
small k, the right hand sides of the flow equations therefore indeed do not depend 
explicitly on k. This is a very nice property, since it allows to extend the numerical 
investigation to very small values of k without the complications usually arising 
when dealing with very small numbers in numerical calculations. 

The possible zeroes in the denominators of the terms from the fermionic sector 
do not pose any problems, since they have the form 

tanhx , X — sinhxcoshx 2 
hm = 1, hni = --. (7.73) 

x^O X x^O x^ 

Note that the emergence of the two different terms in the bosonic sector 
in the broken phase enclosed by the square brackets, one with denominator 1 
and one with denominator (1 + 2A^a22 fc/((27r)^vr^))", is a consequence of the 
Goldstone bosons appearing in the broken phase. The first term stems from the 
two Goldstone modes and the second from the massive mode. 

It is illuminating to discuss qualitatively the flow in the broken phase. In the 
broken phase k is very small. As can be numerically confirmed, the anomalous 
dimension rik is also small. In the flow equation for A^ in the broken phase, we 
can therefore approximate (1 — r/^) ?a 1 and 2 — r]k{...)^2. The quartic coupling 
becomes large, so that (jqi^J^ — jja ~ 0. If we keep ha fixed, then (since Z^ becomes 
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large) ha becomes very small, so that we can neglect the contributions from the 
fermionic sector. In this case, the flow equation for is 

dAt = {K)\---)-'^K- (7.74) 

This means that for decreasing k and initially increasing A^, we eventually reach 
some value of k where the right hand side vanishes, which means that we have 
reached a fixed point of A^. Assuming that this fixed point has been reached, we 
can regard A^ as constant with regard to the other flow equations. For small k 
(so that the momentum integral can be performed), small 7]k and large )\'k(^22k 
(A^ at its fixed point), the flow equation for the potential minimum reads 

T 

dt(y22,k = 2 ■ - a^22,kVk (7.75) 

We emphasize the factor 2 to indicate its origin from the presence of 2 Goldstone 
modes. Inserting r/^ from ( [7.49[ ), we find 



dta'22,k = 2 • - = T^- (7.76) 

AsTT" AaTT^ AsTT^ 

The anomalous dimension correction exactly cancels the contribution of one of 
the two Goldstone modes! This is a well known feature which has already been 
discussed in the context of 0(A^)-models in The right hand side of the flow 



equation in this simple approximation is positive. Thus the flow tends to lower 
•^22 k decreasing k. If 0:22 k becomes too small, our approximation A^a22 k ^ 
breaks down. However, if the minimum goes to zero, we can expect the system 
to return to the symmetric phase for small k. We will see that the numerical 
results confirm this expectation. 

The above discussion was based on fixing the unrenormalized coupling ha- 
Instead, we can fix the renormalized coupling ha- In this case, the fermionic 
fluctuations will not be suppressed in the broken phase. As in the mean field 
case, the fermionic contributions will stabilize the symmetry breaking. We expect 
that in this case we will have true symmetry breaking and do not return to the 
symmetric phase even for very small k. Again, this expectation is confirmed by 
our numerical results. 

Numerical results 



We solve the flow equations ( |7.67|) and ( [7.68| ) numerically. We set h^ = 10 and 



A3 = /i-a/vr^ at the beginning of the flow k = A. We choose A large, so that the 
results no longer depend on the actual choice of A in the limit A; — > 0. One finds 
that it actually suffices to set A = 10. The Yukawa couplings are kept fixed. 
There are two ways to do so: Either we fix the unrenormalized coupling ha or 
the renormalized coupling ha- We give the results for both cases. We find that 



7.4. Results for different truncations of the effective potential 



117 



3 - 



1 - 



- 




-1 I 1 1 > 1 1 > > 1 1 1 

0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

-t 

10-3(m^)^ lO-^A^-- Zk Vk 

Figure 7.1: The Row of the mass (m^)^, the quartic coupUng A^, the wave 
function renormahzation and the anomalous dimension rjk in the symmetric 
regime for the antiferromagnctic boson S2 at half filling and temperature T = 
0.15. The Yukawa coupling is h\ = 10. In this plot we keep the unrenormalized 
Yukawa coupling hxed. 
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10-\n4)^ 10-6 A^-- Zfc r)k 

Figure 7.2: The How of the mass {m'^)^, the quartic coupling A^, the wave 
function renormahzation Z^ and the anomalous dimension r]k in the symmetric 
regime for the antiferromagnetic boson S2 at half filling and temperature T — 
0.15. The Yukawa coupling is — 10. In this plot we keep the renormalized 
Yukawa coupling fixed. 
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ao,k— 10-^A^-- 10-^Zfc 

Figure 7.3: The Sow of the minimum aok, the quartic coupling and the 
wave function renormahzation in the regime of broken symmetry for the 
antiferromagnetic boson S2 at halffiUing and temperature T — 0.15. The Yukawa 
couphng is — 10. In this plot we keep the unrenormalized Yukawa coupling 
fixed. 
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Figure 7.4: The How of the mass ctofe, the quartic coupling and the wave 
function renormalization in the regime of broken symmetry for the antiferro- 
magnetic boson S2 at half filling and temperature T = 0.15. The Yukawa coupling 
is hi = 10. In this plot we keep the renormalized Yukawa coupling fixed. 
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Figure 7.5: The Bow of the minimum «ofc in the regime of broken symmetry 
for the antiferromagnetic boson S2 at half fiUing for different temperatures. The 
Yukawa coupling is — 10. In this plot we keep the unrenormalized Yukawa 
coupling fixed. 
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Figure 7.6: The flow of the minimum aok in the regime of broken symmetry 
for the antiferromagnetic boson S2 at half filling for different temperatures. The 
Yukawa coupling is /i^ = 10. In this plot we keep the renormalized Yukawa 
coupling hxed. 
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the critical temperature that describes the onset of symmetry breaking for 
T < Tc is slightly decreased in comparison to the mean field case (Tc ^0.2 here). 
In figs. [7.1H7.4I we plot the results for the flowing variables at fixed temperature 
T = 0.15. For this temperature, we find spontaneous symmetry breaking if k is 
sufficiently small. In figs. [7.1| and [7.2| we plot the flow of the mass (m^)^, the 
quartic coupling A^, the wave function renormalization Z| and the anomalous 
dimension rjk. We see that the mass becomes large if k is lowered, reaching 
some maximum value and then drops to zero, indicating the phase transition. 
Note that we plotted 10~^(m^)^, so that the initial value of the mass (which is 
(m^)^ = (27r)^/i^/A^) is not distinguishable from zero in the plot. The quartic 
coupling also reaches a maximum during the flow in the symmetric phase and 
begins to decay if k is lowered towards the value where the phase transition occurs. 
Since the bosonic fluctuations (which enter the flow equation for the mass as a 
term oc A^) tend to prevent the phase transition, this behavior was to be expected. 
The wave function renormalization becomes large and remains large at the phase 
transition. This behavior will be important to explain the apparent contradiction 
to the Mermin- Wager theorem. In the symmetric phase the qualitative behavior 
of all quantities is independent of whether we fix the unrenormalized or the 



renormalized coupling. This is quite different in the broken phase. In fig. 
and fig. |7|^ we show the flow of the minimum a22,fc; quartic coupling A^ and 
the wave function renormalization in the two cases. In both cases, the quartic 
coupling reaches some fixed point and diverges if k is lowered. If we keep 
ha fixed, the minimum reaches some maximum value and returns to zero. If we 
keep ha fixed, the maximum converges to some finite value. Both results are in 
agreement with our expectations based on analytical reasoning. 
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In the case of fixing the unrenormalized coupling, no contradiction to the 
Mermin- Wagner theorem appears, since for — > the symmetry becomes again 
unbroken. The symmetry breaking in some finite range of k can then be in- 
terpreted as an antiferromagnetic order on large clusters that disappears if we 
average over even larger scales. Note the scale at which the broken symmetry 
becomes again unbroken, t = —180 corresponds to k = Aexp(— 180), which is 
extremely small, so that the symmetry becomes unbroken only when averaging 
over extremely large scales. For any probe of realistic size to be examined ex- 
perimentally, we will find antiferromagnetic properties. In this interpretation, 
antiferromagnetic properties of superconductors do not contradict the Mermin- 
Wagner theorem, since they are finite size effects that would disappear if the 
probes were made large enough. 

However, if we fix the unrenormalized coupling, the symmetry remains bro- 
ken for k —>■ 0. To see how this can be reconciled with the Mermin- Wagner 
theorem, we have to be careful in distinguishing the renormalized and unrenor- 
malized minimum of the potential. The Mermin- Wagner theorem states that the 
unrenormalized minimum has to vanish for A; — > 0. But we have analyzed the 
flow of the renormalized minimum, and since the wave function renormalization 
diverges for k ^ 0, k = '^22 k/^k actually vanishes in complete agreement with 
Mermin and Wagner. In fact, this is the mechanism how phase transitions can 
actually appear even in the case where they are forbidden by Mermin- Wagner 



(Kosterlitz-Thouless type phase transitions |22 



In figs. ^]5|, |7.6| we show how the flow of the minimum changes as a function of 
temperature. For fixed unrenormalized coupling, we see that the strength of the 
symmetry breaking and the scale at which the broken symmetry becomes unbro- 
ken increases if the temperature is lowered. This is what we intuitively expect, 
since the antiferromagnetic clusters should be larger for small temperature. For 
fixed renormalized coupling, the maximum value of the minimum also becomes 
larger for smaller temperature. Note that all the curves converge, but those for 
low T on much larger scales than shown in the plot. Fig. |7[^ enlarges the region 
of small —if: from 7^6, which is actually accessible by experiments. We see that 
the temperature dependence is weak in this region. 

In the framework of our truncation, we cannot decide which of the two possi- 
bilities of reconciling the occurrence of spontaneous symmetry breaking with the 
Mermin- Wagner theorem is realized. To do this, we have to include the flow of 
the Yukawa couplings. 



7.4.2 Charge density fluctuations and superconductivity 

In this section we will investigate the influence of charge density fluctuations on 
the superconducting properties of the theory. We set ^22 = and keep all wave 
function renormalizations constant. 
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Fermionic contribution to the flow 

The flow equation for reads 



In cosh (^/2^^ + 2ie2(ci + 6102))' + 2/i4A(ci - eic^)^ 



€162 



with 



^,,,(p,A) = — ^(^2^+2ie2(ci + eiC2)) + 2/i4A(ci - £102)2. (7.78) 

Although we keep all wave function renormalization constants flxed, we write 
instead of h^i since the initial value of 7^ 1. 



Truncation of the potential 

Our truncation for the potential reads 

K1K2 
K1K2K3K4 
K1K2K3K4 

{pu{Ks, K,) - p1,S{Ks)S{K,)){pn{K,, K^) - pl,S{K,)S{Ke)) (7.79) 
in the symmetric phase. For homogeneous flelds this becomes 

Uk^V (^{mifA + ^XtA' + ^Kipn - pl,f + \k,A{p,, - pl,f^ . (7.80) 

We included terms up to quadratic order in A just as in the antiferromagnetic 
case. Additionally, we expanded pu around the minimum pJi- Recall that this 
minimum is not /c-dependent, since the charge density is an external parameter 
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controlled by the doping. The term oc kk induces an interaction between the 
charge density and A. 

We will restrict our attention to the flow in the symmetric regime, so that we 
do not write down a truncation for the potential in the broken phase. 



Bosonic contribution to the flow 

In much the same way as in the antiferromagnetic case we find 

^P^iQ) + 3X^A + h,{pu - p'n?^ (P.^(g) + (A^ + «:.A)(3pn - P?i)) 
- AklApuipu - p'u)'] , (7.81) 



where 



P^ = Plk + Rl + i^kf, P^ = Plu + Rk- (7.82) 



Choice of the cutoff functions 

We define the cutoff functions to be 



^Rl{Q) = R§{Q) = {2n)\\e - {mTme - {mTf). (7.83) 



Z 



In contrast to the momentum dependent cutoff we chose in the last section, 
this cutoff allows to perform the Matsubara sum analytically, which speeds up 
the numerical calculation significantly. The drawback is that we can no longer 
interprete the fiow as an averaging process over larger and larger clusters, but 
since in this section we will be mainly interested in whether a phase transition 
takes place for some temperature and p^i? we do not need this interpretation. 
Instead of inserting the cutoff function into ( [7.81 ) immediately, it is more 



convenient to first extract the fiow equations for the couplings and masses. 
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Extraction of the coefficients 

Prom the flow of the potential we can obtain the flow equations 



V A^o.pii^p?, dA \dk 



V A^o.pii^p?! dA^ \dk 

V A-vo,pii-»p?i dpi^ \dk 
dki^k = ^ . lim „ -rfj^ (i-U, ) . (7.84) 



V A-^o,pii-*p?i dpl^dA \dk 



Introduction of rescaled quantities 

The rescaled quantities read 



) - — ' ^k - \- -j^^ i^k - -j^- (7.85j 



The flow equations 

The flow equations are 



-(2.)T(2M + l)y^ — (^-^ + ^J-2K) 

'scchV6i62 (1 + 2/,2,,2sech2/.i62) tanh/,,,2 



/■2 fi 



_Th^ r ^( (l + 2/,^^,^sechV.,.Jtanh^,; 



+ (2.)T(2M+l)y^ 1^1W"WFJ"''^ 
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+ {2TifT{2M + 1 



2 cosh /. 







V (^;?)' (n^)^ (Pi^r {p^yp^ p^{p^y 

2nk (7.86) 



with M = max{m e N|m < /c/T}, 



i^.^ = 2n% ^-2~"' + 47r^ + 2A^p?„ (7.87) 



where p^^ = p^/(2/i^). The initial conditions are (m^)^ = Sn^hl/A^ and A^ = 
K,A = 0. = 1/A^ is kept fixed. The initial condition for A^ is set to Af^ = 
since the one loop correction to this coupling is of this order (we do not set A^ = 
in order to have a potential minimum in 



Remcirks 

As in the antiferromagnetic case, the zeroes of the denominators of the fermionic 
contributions to the flow do not cause any problems, since they can be canceled 
against numerator zeroes. 

Note that the contributions in the bosonic sector have a simple diagrammatic 
interpretation. To see this, note that in our flow equations the derivatives dk 
have been carried out on the right hand side (which is of course necessary for 
the numerical investigation). For the diagrammatical interpretation, it is more 
useful to consider the expressions before taking the derivative. For example, the 
bosonic contribution to the flow of {m^Y is 
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Figure 7.8: Dependence of the superconducting region (including A- and p- 
Huctuations) on the Yukawa couphng hp. We take h\ — 20. 

and similarly for the other variables. The diagrammatic interpretation is then 

mt: — — + - 



+ — 




K 




+ 





(7.89) 



where the dashed line represents the propagation of Xi the solid line the prop- 
agation of Ri and the double line the coupling to the external charge density 



Numerical results 



We have solved the flow equations numerically for different temperatures, charge 
densities and Yukawa couplings hp. We set /i^ = 20. We follow the flow of {m^y 
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until it either reaches zero or diverges. The first case indicates an instability in 
the (i-wave superconducting channel (this should be taken with a grain of salt, 
since we know from our mean field discussion that phase transitions of first order 
render this criterion untrustworthy. However, in the superconducting regime we 
expect no phase transitions of first order to occur if no antiferromagnetic order is 
present). In fig. |7.8| we show the borders of the regions where the mass happens to 
vanish during the flow. The comparison between different values of hp shows that 
a strong charge density coupling tends to enlarge the region of superconductivity. 
Note that we use the parameter hpp on the horizontal axis instead of p in the 
mean field case. Recall that in the mean field case the results depended on hp 
only via hpp. This means that if we ignored the bosonic fluctuations (which 
would reduce our flow equations to the mean field case), no difference between 
the phase borders for different hp would appear in our plot (a comparison of the 
mean field results and the results for the inclusion of various bosonic fluctuations 
can be found in the next section). The differences are consequences only of the 
inclusion of bosonic fluctuations. The choice of the couplings in our plot is rather 
extreme; in the mean field case changing /i^ or h\ only by a factor 2 had a strong 
effect on the phase diagram. By changing /i^ by a factor of 2 (not 20 as in the 
plot), we see that the effect of charge density coupling is relatively small. We 
will see in the next section that — in contrast to the charge density — the phase 
diagram is very sensitive to the strength of antiferromagnetic fluctuations. 



7.4.3 Antiferromagnetic fluctuations and superconductiv- 
ity 

In this section, we will analyze the effect of antiferromagnetic fluctuations on the 
phase diagram. The derivation of the flow equations is very similar to the case 
of including charge density fluctuations in the flow. We set pn = p^i (that is, we 
do not consider charge density fluctuations in this section). 

Fermionic contribution to the flow 

The flow equation for U[ reads 

^U^ = ^^J2[ 7Si^,.2(«22,A)tanh^,,,(a22,A) (7.90) 

with 



/(«22,A) = 




(7.91) 
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Truncation of the potential 

We truncate the potential as 

K1K2 



K1K2K3K4, 

+ J2 «)'«22(i^l, K2)6{Ki + K2) 



K1K2 



+ 1 Yl K<^22{K,, K2)a22{K3, K4)6{K, + K2 + + K^) 
K1K2K3K4 

+ Yl %Oi22{Ki,K2)A{K3,K^)S{K, + K2-Ks + K4). (7.92) 
K1K2K3K4 

This is an expansion up to quadratic order in 0:22 and A. Again, we restrict our 
attention to the symmetric phase. 

Bosonic contribution to the flow 

The flow of the potential in the bosonic sector is given by 

+ "^^.T Y.j ^^J^2 + A^22 + 7. A) 



m 



2tt 



(27r)2 



where 



In {{{P^ + 3A^A + 7fca22)(P|'+ 3A^a22 + 7feA)) - 47|Aa22) 

(7.93) 



Choice of the cutoff' function 

We choose the same cutoff function as in the last section 
1 



i^RHQ) = RliQ) = {27:)\\e - {mTf)e{e - {mTf). (7.95) 
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Extraction of the coefficients 

Prom the flow of the effective potential we obtain 

a,(m^)2 = l hm 4^(^U, 

V A,a22^0dA \dk 

OkK = — hm — r —Ui 



= ^ lim {-^Uk 
dklk = :J: lim ^ {-JT^k 



(7.96) 



Introduction of rescaled quantities 

The rescaled quantities read 

ft ft ft ft ft 

The flow equations 

We find the flow equations 



£162 



SechVeiea (1 + '^fle^^^^^ U^e^) tauh /, 



/•2 f3 



eie2 J €162 

21-^2^ /0/17X/'\A\2 e^,2 
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dtik 



5 ^ / ^ (27r)2 
hpp 



4T| 



1 +62 



2t(ci + eiC2); V fl 



[ sechV.iea (1 + S/.^i.^sechVciea) tanh 



£2 



+ (27r)^T(2M+ 1) 



2- A /22(A^)2 



fc2 



{Pi? 



T5 



4T| 



(Ci - eiCs)^ 1 + 62 



/ sechVei 



hpp 



^ (27r)2^^^ ^^^^^ ' ^'2t(ci + eiC2) 
(1 + 2/2 ,^sech2/.ie2) tanh/,,, 



^2 



V P 

+ (27r)2T(2Af + 1) 



27r 72 



A;2 



^12Z,^7.A^ ^ 107fcAg ^ 4Z,^7fc^^ ^ 



47| 



(p^r 



k > fc 

with M = max{m G N|m < fc/T}, 

1 



2Tfc 



(/ipP + 2te2(ci + eiC2)), 



-I I, 



27r^A, 



A;2 
ci - C2 



,A\2 
k I 1 



fc2 



o\2 



(7.98) 



(7.99) 



The initial conditions are {m^Y = Sn'^hl/A'^, (ml)^ = in^hl/A'^, = A 
ka = 0. 



Remarks 

Again the flow equations seem to be ill deflned for /,,,, ~^ and again all these 
denominator zeroes are canceled by corresponding numerator zeroes. However, 
we additionally face possible singularities in the flow equations which involve 
derivatives with respect to 0:22 in ( |7.96| ) for ci + eiC2 0. These singularities 
cancel if the sums over ej are performed. For example, 



£2 



hpp 



2t(ci + eiC2) / V 



-4sech^^^ ( tanh 

V _ 1^ (7.100) 
2Tfc \2Tk 2Tfc / 
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Figure 7.9: Boundaries of tiie superconducting region for tlie mean field case, 
the inclusion of bosonic A- fluctuations, A- and a-fluctuations as well as A- and 
p- fluctuations. The Yukawa couplings are = 20, = 10 and h"^ = 20. 

for Ci + eiC2 —>■ 0, which is perfectly finite. In the same way, the sum over €2 in 
the flow equations for and 7fc can be carried out in the limit ci + eiC2 0. 

Note that as in the last section, the contributions in the bosonic sector have 
simple diagrammatic interpretations: 




The dashed line represents the propagation of whereas the solid line the prop- 
agation of s. 

Numerical results 

In fig. |7.9| we compare the effect of including different kinds of bosonic fluctuations 
into the flow with the mean field result. The plot for A- and a-fluctuations has 
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Figure 7.10: Dependence of the superconducting region (including A- and a- 
Huctuations) on the Yukawa couphng ha- We take h"^ = 15. 



been generated by using the flow equations of this section, the plot for A- and 
p-fluctuations by using (|7.86|) . By setting (m^)^ = (m^)^, = in (|7.98|) we 
get the plot for the inclusion of A-fluctuations only. We set /i^ = 20, = 10 
and h"^ = 20. For this choice of the Yukawa couplings only the superconducting 
phase is present in the mean field result. The same is true if bosonic fluctuations 
are included, but compared to the mean field case, the superconducting region 
becomes smaller. This was to be expected, since the fermionic fluctuations favor 
the phase transition, whereas bosonic fluctuations have the opposite effect. In 
the mean field approximation we ignore the bosonic fluctuations, so that we 
overestimate the symmetry breaking behavior, which leads to a larger region 
of broken symmetry in the phase diagram. For this special choice of Yukawa 
couplings, the boundaries of the superconducting region approximately coincide. 
As we have seen in the last section, the phase boundary is not very sensitive 
to changes in hp. However, fig. [7.10| shows that we have a strong dependence 
on ha — as we already had in the mean field case. In so far, the inclusion of 
bosonic couplings do not qualitatively alter the mean field results. The interesting 
feature of the plot is that the superconducting region is shifted to the right 
if ha is increased. Of course, this is a feature not present in the mean field 
approximation, since there ha entered only via the combination h'^a, so that 
outside the antiferromagnetic phase (where a = 0) varying ha could not have 
any effect on the boundary of the superconducting region. This means that 
strong antiferromagnetic coupling tends to enlarge the superconducting region 
by means of antiferromagnetic fluctuations, even if there is no antiferromagnetic 
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order present! In fact, it is suspected that antiferromagnetic fluctuations are 
crucial for the understanding why cuprates remain in the superconducting state 
even for large temperature. 

Also note that — as in the last section — the phase boundaries have been 
inferred from an analysis of the masses only. Therefore our results have to be 
taken with some care — possible phase transitions of first order will shift the 
phase boundaries. This is particularly true for the left boundaries of the super- 
conducting regions in fig. [7.10|, where the superconducting region is bordered 
by a region of antiferromagnetic behavior, since in our mean field calculation we 
found phase transitions of first order exactly at this boundary. However, at the 
right boundaries no phase transitions of first order appeared in the mean field 
approximation, so that hopefully our interpretation remains intact even for the 
more general case considered here. 



7.4.4 Final remark 

Although many questions remain and a lot of work has to be invested to free our- 
selves from the limits of the problem of the coupling ambiguity and the possibility 
of first order phase transitions, the point we want to make is that we see that our 
formalism is in principle suitable for analyzing the properties of the model both 
in the symmetric and in the broken phase. Renormalization group approaches so 
far have been limited to the investigation of the fiow in the symmetric phase and 
were not able to describe features of the fiow beyond the point of spontaneous 
symmetry breaking. The problems we face in our analysis so far are not intrinsic 
to our formalism and can be cured by more refined truncation schemes (partic- 
ularly including the fiow of the couplings and some more general truncation for 
the effective potential). A systematical enhancement of the truncation scheme is 
straightforward and will be attacked in the future: 

• In our approach so far, we approximated the effective potential by a polyno- 
mial in the fields, which raises the problem of dealing with phase transitions 
of first order. Calculations including the flow of the full potential are possi- 
ble and have already been carried out successfully for a number of systems 
(cf. e. g. 0). 

• The phase diagrams we calculated in this chapter depend on the choice of 
the initial values of the Yukawa couplings that we kept constant. To get 
rid of this dependence, we have to include the flow of the Yukawa coupling 



into our sets of flow equations. This task is already worked on |17 . 



One qualitative feature of the phase diagram of high temperature super- 



conductors flg. 1.1 is the separation of the antiferromagnetic and super- 



conducting region at low temperature and intermediate doping. In our 
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calculations, the superconducting region is always bounded by the anti- 
ferromagnetic region towards small p. The reason for this shortcoming is 
our oversimplified homogeneous treatment of the charge density More 
reahstically, the charge density caused by doping is not homogeneously dis- 
tributed, but forms stripe shaped regions of alternating high and low charge 
density. The width of the stripes depend on doping, and it turns out that 
for some intermediate doping the alternating charge density induces the 
formation of parallel spin ladders that decouple from each other. The effect 
is that long range order is lost, and antiferromagnetic order occurs only in 
the ladders. The material then becomes paramagnetically over large scales. 
This behavior should be reproducible with our formalism if we give up 
the assumption of homogeneous charge density and generalize it to stripe 
structures. 

• Although we introduced a large set of bosons in the formalism, we concen- 
trated on only a small subset of them in our calculations. The choice of 
this subset was motivated by experiment: We know that high temperature 
superconductors exhibit antiferromagnetic and (i-wave superconducting be- 
havior and therefore chose exactly those bosons representing these proper- 
ties. It is interesting to see whether a more unbiased choice, that is, taking 
into account more bosons representing e.g. s-wave superconductivity or fer- 
romagnetism, confirms that we have chosen those bosons actually leading 
to spontaneous symmetry breaking. 

Some words are in order comparing our formalism to the renormalization 



group approaches in |I8|. As already mentioned, all these calculations were per- 
formed in the purely fermionic theory. The authors basically investigate the flow 
of quartic fermionic couplings in different channels. The conventional way of 
these approaches was to introduce a regularization scheme that cut off momenta 
near the Fermi surface. In this scheme, following the flow to small k corresponds 
to the inclusion of momentum modes increasingly close to the Fermi surface. 
More recently, temperature cutoffs resembling the regularization scheme used in 
our work have been applied (Honerkamp, Salmhofer and Rice 2002 in [0). In- 
dependent of the regularization scheme, the flow of quartic fermionic couplings 
indicates instabilities in certain channels by the emergence of divergencies of the 
corresponding couplings. This means that in these approaches it is only possible 
to follow the flow until the point of symmetry breaking is reached, and no infor- 
mation about the behavior in the broken phase is available. By the bosonization 
and the investigation of the effective potential in this work, it becomes possible 



to analyze the flow in the broken phase. However, in contrast to , our formal- 
ism is plagued by the problem of the coupling ambiguity, and its success will be 
ultimately measured by our capability to overcome this complication. 

It should have become clear in the course of the explicit calculations presented 
over the last sections that our formalism provides an elegant and suitable starting 
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point for the investigation of the Hubbard model and the properties of high 
temperature superconductors, and is also easily implemented for investigations of 
the behavior inside the broken phase. At this point, we have reproduced the gross 
qualitative features of the phase diagram of high temperature superconductors, 
with correctly placed regions of antiferromagnetic and superconducting behavior. 
We understand how the Mermin- Wagner theorem can be reconciled with the 
existence of antiferromagnetic long range order for non vanishing temperature, 
and we took a first glimpse on how antiferromagnetic fluctuations can enlarge 
regions of superconductivity. All these results let us suspect that we are on the 
right way and we hope to provide more insights into the nature of the Hubbard 
model and high temperature superconductors in the future. 
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Conclusion 



High temperature superconductors have a two dimensional layer structure with 
small interlayer coupling. They can be modeled by the two dimensional Hubbard 
model on a square lattice. This model describes electrons on a quadratic lattice 
that experience a local Coulomb interaction and are able to hop to adjacent lattice 
sites. The partition function for this model depends on the chemical potential 
(or equivalently the charge density), the temperature and the relative strength of 
Coulomb interaction and hopping amplitude. To compare the predictions of the 
model with experimental results for high temperature superconductors, one has 
to calculate the phase diagram of the model in the charge density-temperature 
plane. 

Our way to tackle this task is to identify the most prominent degrees of 
freedom of high temperature superconductors, which are antiferromagnetism and 
d-wave superconductivity, and to define a set of bosonic "particles" , so that every 
particle corresponds to one degree of freedom of the model. We found an exact 
transcription of the partition function, which describes a Yukawa-like theory. The 
expectation values of the bosons in this rewritten theory indicate a possible long 
range order in the channel (antiferromagnetic, ci-wave superconducting, etc.) to 
which the bosons correspond. 

A mean field calculation in this partially bosonized theory, neglecting all 
bosonic fluctuations and integrating out the fermions, already reveals the main 
features of the phase diagram of high temperature superconductors. 

More refined calculations can be performed by means of exact rcnormalization 
group equations. We use the effective average action method. To simplify the 
definition of truncation schemes and to minimize the error induced by approxi- 
mations, we rewrite the partially bosonized theory as a function of bosons which 
are eigenstates of translations on the lattice. Due to the lattice symmetries, these 
bosons are no longer mixed in the full effective action (e.g. the full propagator 
matrix becomes diagonal). 

We use this modified theory as a starting point for a renormalization group 
analysis. This analysis shows how the Mermin- Wagner theorem can be recon- 
ciled with the existence of antiferromagnetic long range order for non vanishing 
temperature and indicates that antiferromagnetic fluctuations tend to favor su- 
perconducting behavior in certain regions of the phase diagram. 
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The drawback of the present analysis is mainly the arbitrariness of the Yukawa 
couplings. Current work is dedicated to this problem. However, the cure is basi- 
cally an improved truncation scheme for the effective action and does not intrinsi- 
cally limit our formalism. Planned future work includes a more general treatment 
of the effective potential with special regard to first order phase transitions and 
the inclusion of more bosons in the truncation scheme, allowing to test whether 
e.g. superconducting channels other than that with spatial d-wave symmetry 
play a role. 

In conclusion, we hope to have provided a formalism which is easy to imple- 
ment for renormalization group studies, which introduces a convenient interpreta- 
tion of non local fermionic order parameters as local expectation values of bosonic 
fields and that will continue to help investigating the rich and beautiful spectrum 
of properties of the Hubbard model and high temperature superconductors. 



Appendix A 
Conventions 



We use units for which h = c = = 1- Bold symbols (n, x, q, etc.) denote two 
dimensional vectors. Symbols with arrow (a, m, etc.) denote three dimensional 
vectors. Generalized momenta are called Q, P and K, whereas X, Y and Z are 
generalized positions (see below for the definition of generalized quantities). By 
a " we indicate fields. The same symbol without " is the expectation value of the 
corresponding field. ~ is used to indicate fermion bilinears (to distinguish them 
from their bosonic counterparts). 



A.l Fourier transforms 

We define 



The Fourier transforms for both fermionic and bosonic fields are given by 



Q = {uJn,q), X^{T,n), 
QX = ujnT + nq, 




(A.l) 



5{Q-Q')^P5nA'^T:f5{q-q'), 
5{X - X') ^ 5{t - T')5r.,^> . 



Xa(^) = ^Xa(Q)exp(z(gX + 2„q)), 



Q 



(A.2) 



xliX) = xliQ) exp(-z(gX + z„g)), 



Q 
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where x stands for ip, or Uis, whereas x* stands for ip* or u*^ and the Za are 
given by 



Zi = (-l/4,l/4) Z2 = (l/4,l/4) 

Z4 = (-1/4, -1/4)^3 = (1/4, -1/4). ^ • ' 



A. 2 Matrices 

{(Ji}, i G {1,2, 3} is the usual set of Pauh matrices. Additionally, we identify ctq 
with the unity matrix. We then define the matrices o"^ <S> (7v by 



-^-i's:y --C^O' ^-(AtO' ^-Co 4 

(A.4) 



where ji G {0, 1, 2, 3}. The matrices and have the properties 

{Ai, Bj} = 2S,jBo, {A, Aj} = {Bi, B^} = 2Sij, 
[A, Bj] = 2ie,,fcBfc, [A, Aj] = [B,, Bj] = 2ie,jkAk, (A.5) 
BqA^ = A^Bo = Bi^, BoB^ = B^Bq = A^ 

with i G {1,2,3}. 



A. 3 Fermion bilinears 



XabiX) = tijl{X)ia2MX) 
X:,{X) = -i,l{X)ia,r,{X). 



(A.6) 



A. 3. Fermion bilinears 
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^11 + ^22 + ^33 + 5'44 


— * 


= 'fu + ^22 + <^33 + <^44 


p = 


C^ll — 5-22 + 5-33 — 0-44 


— * 

a 


= '^ll - '^22 + ^^33 - '^44 


% = 


5^11 + d'22 — 5"33 — a44 


9y 


= "fw + ^^22 - '^33 - <^44 


Qx = 


&U - O22 - 5-33 + Cr44 


9x 


= - <^22 - <^33 + <^44 


s = 


Xll + X22 + X33 + X44 


e 


= X12 + X23 + X34 + X41 


c = 


Xll - X22 + X33 - X44 


d 


= X12 - X23 + X34 - X41 


ty = 


Xll + X22 - X33 - X44 


Vy 


= X12 - X34 




Xll - X22 - X33 + X44 


Vx 


= X23 — X41 


S* = 


Xii + X22 + X33 + 


e* 


= Xt2 + X23 + X34 + Xll 


c* = 


Xii ~ X22 + X33 ~ X44 


d* 


= X12 ~ X23 + X34 ~ X4I 


f* = 


Xii + X22 ~ X33 ~ X44 




— X12 ~ X34 


f* - 

^x 


Xii ~ X22 ~ X33 + X44 




— X23 ~ X41 



(A.7) 
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Appendix B 

Vertex factors for the Hubbard 
model 



The vertices V^{Q', Q") for the bosons p, q^^y depend only on the momentum 
Q = Q' — Q" . With = (1, 0), = (0, 1) and Za, a = 1 ... 4, given in appendix 
|A.1| , they can be written in the form 

KIc(Q', Q") = KZc(Q) = ^e--v-^M- ,(g) ® rr, (B.i) 

The color matrices (with matrix elements ^) read 

MC(Q) = diag{l, e'^^'^, e'^^^'-^y^'^ , e^'^^''}, M^{Q) = diag{l, 1, 6"'^^'^, 6"^'='"'}, 
M^iQ) = diag{l, 1, 1, 1}, M^iQ) = diag{l, e'^^'^, e'^^\ 1}; 

(g) = {-ly-' diag(l, -1, 1, -1) MP{Q)- 
Ml^{Q) = Mf (Q) ■ diag(-l, -1, 1, 1), Ml^iQ) = M^{Q) ■ diag(-l, -1, 1, 1), 
M|^(Q) = M^iQ) ■ diag(l, 1, -1, -1), (Q) = M^(Q) ■ diag(l, 1, -1, -1); 

Mf (g) = Mf(g) ■ diag(-l, 1, 1, -1), (Q) = M2^(Q) ■ diag(l, -1, -1, 1), 

M|^(Q) = MP{Q) ■ diag(l, -1, -1, 1), Mf (Q) = M^{Q) ■ diag(-l, 1, 1, -1). 

(B.2) 

The same can be obtained for the bosons with spin index, m, a,gx y, by substi- 
tuting 1'/'' ^ a'P'''. 

For the bosons s, c, t^^y one finds similarly: 

vXAQ'.Q") = ^e-'"('''+''")e--('''+'^")Mr,,(g',g") ® ^^2, 

Mf{Q',Q") = M^{-Q'-Q"), Mf{Q',Q") = M^{-Q'-Q"), 
Mf (Q', Q") = i-Q' - Q"), M'J (Q', Q") = Mf (-Q' - Q"), 

while (i, e, Wa;^^ are a bit more complicated. Let us define e*-' = e*(^''''+^j''") and 
a *-product C = A * _B by Cij := AijBij (no sum over indices here!). One then 
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Appendix B. Vertex factors for the Hubbard model 



obtains 
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(B.4) 



With the aid of the *-product the other vertices can now be obtained from these 
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m:^(q',q")^ 



f -1 \ 
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/ 1 \ 

10 

-1 
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*M,^ (g',g"); 



*Mf{Q',Q")-X,, A =(-1,1, 1,-1); 



*Mf{Q',Q").X,, A =(-1,-1, 1,1). 



(B.5) 



The transition from the 'U*'0'0- vertices to the 'U'0*'0*-vertices can be 
carried out by 

V""*(<9',Q") ^ ^"(<3',Q") = -(v^"*(<3',Q"))* = -v^"*(-g',-Q")- 
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